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Plan  of  this  Report 

This  report  is  divided  into  two  major  parts.  In  PART  I,  brief 
summaries  of  the  principal  research  results  are  given  together  with 
lists  of  published  papers,  reports,  dissertations,  theses,  and  oral 
presentations  that  resulted  from  the  work.  Also,  the  personnel  who 
worked  on  the  project  are  listed.  In  PART  II,  summary  accounts  of 
three  selected  technical  areas  are  given  and  suggestions  for  fruit¬ 
ful  areas  of  future  work  are  listed.  Many  results  uncovered  during 
the  course  of  the  project  are  not  dealt  with  in  PART  II;  rather, 
three  subjects  representative  of  what  are  felt  to  be  of  principal 
importance  in  the  overall  work  are  discussed.  Further  details  of 
the  technical  results  can  be  found  in. the  papers  and  reports  listed 


in  PART  I. 


1 .  INTRODUCTION 


This  document  represents  the  concluding  report  of  research  results  ob¬ 
tained  on  the  project.  Computational  Methods  in  Nonlinear  Mechanics,  during 
the  period  June  1,  1980  to  September  30,  1982.  The  project  was  supported 
by  the  Air  Force  Office  of  Scientific  Research  under  Contract  F-49620-C-0083 
and  was  monitored  by  Dr.  Anthony  Amos  of  AFOSR. 

All  objectives  originally  listed  in  the  Statement  of  Work  of  the  pro¬ 
posal  for  this  work  have  been  accomplished.  Indeed,  a  number  of  research 
results  which  are  regarded  as  very  important  in  the  study  of  the  nonlinear 
mechanical  behavior  of  complex  structures  by  finite  element  methods  have 
been  established.  These  have  represented  advances  in  specific  areas  of 
numerical  analysis,  computational  algorithms,  computer  codes  for  nonlinear 
structural  problems,  modelling  and  characterization  of  various  nonlinear 
effects,  and  resolving  theoretical  questions  connected  with  nonlinear 
theories  of  structural  behavior.  More  details  on  these  results  are  given 
later. 

During  the  two  year  contract  period,  the  research  effort  led  to  the 
publication  of  33  journal  articles  and  scientific  papers,  11  technical 
reports,  5  Ph.D.  dissertations,  2  M,S.  theses,  and  50  oral  presentations 
at  conferences,  seminars,  and  invited  lecture  series.  The  principal 
investigator  of  the  project  was  Professor  J.  Tinsley  Oden  and  the  co¬ 
principal  investigator  was  Professor  E.  B.  Becker,  both  of  the  University 
of  Texas  at  Austin.  Lists  of  other  personnel  who  worked  on  the  project 
during  the  contract  period  are  given  in  Chapter  3. 


2.  PRINCIPAL  ADVANCES 


Significant  research  results  in  a  number  of  distinct  areas  were  ob¬ 
tained  during  the  course  of  this  project.  However,  it  is  felt  that  espe¬ 
cially  penetrating  and  important  results  in  several  specific  areas  were 
obtained.  In  the  following,  we  list  what  are  considered  some  of  the  prin¬ 
cipal  contributions  made  in  this  research  project. 

1.  Finite  Element  Methods  for  Contact  Problems  in  Elastostatics 

A  large  volume  of  work  was  done  on  the  problem  of  contact  between  deformable 
elastic  bodies.  In  this  general  area,  a  variety  of  results  were  obtained 
on  the  description  of  the  physical  phenomena,  development  of  mathematical 
models  using  the  theory  of  variational  inequalities,  the  study  of  existence, 
uniqueness,  and  regularity  of  solutions  of  the  variational  inequalities, 
development  of  corresponding  approximation  theories,  the  development  of 
numerical  stability  criteria,  error  estimates,  convergence  criteria,  de¬ 
velopment  of  families  of  completely  new  algorithms  for  solving  resulting 
systems  of  linear  and  nonlinear  inequalities,  development  of  programming 
strategies  and  finite  element  codes  for  the  analysis  of  such  problems,  and 
the  actual  analysis  of  numerous  example  problems  drawn  from  elasticity, 
plasticity,  and  elastoplastlcity  in  which  contact  conditions  are  encountered. 
Included  in  this  large  collection  of  results  are  several  completely  new 
variational  principles  for  contact  problems  in  elasticity,  new  finite 
element  methods  which  employ  such  devices  as  Interior  and  exterior  penalty 
methods,  reduced  integration,  mixed  finite  element  formulations,  and  techniques 
drawn  from  linear  and  nonlinear  programming. 

The  theoretical  results  were  not  limited  to  simple  unilateral  contact. 


Indeed,  new  results  were  obtained  for  prc 
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buckling  of  thin  elastic  plates  and  shells  supported  by  unilateral  constraints 
and  studies  were  made  of  a  complex  family  of  bifurcation  problems  characterized 
by  nonlinear  variational  Inequalities. 

Many  of  the  results  obtained  in  this  general  area  are  summarized  in  a 
forthcoming  treatise  entitled,  "Contact  Problems  in  Elasticity,"  by  the  prin¬ 
cipal  investigator  and  Professor  N.  Kikuchi,  who  worked  on  the  project 
periodically  during  the  contract  period.  This  volume  is  to  be  published  by 
SIAM  Publications  (Society  of  Industrial  and  Applied  Mathematics,  Philadel¬ 
phia,  PA)  and  should  appear  in  early  1983.  It  is  believed  that  there  has 
never  before  been  a  more  detailed  and  thorough  analysis  of  this  difficult 
class  of  nonlinear  problems  in  structural  and  solid  mechanics.  A  brief 
summary  of  some  of  these  results  is  given  in  chapter  4  of  this  report. 

2.  Reduced  Integration-Penalty  Finite  Element  Method  for 
Constrained  Problems  in  Elasticity  and  Fluid  Flow 

Some  n”merical  instabilities  encountered  early  in  the  project  led 
the  investigators  into  the  study  of  a  collection  of  difficult  numerical 
and  theoretical  issues  connected  with  the  convergence  and  numerical  sta¬ 
bility  of  a  variety  of  popular  finite  element  methods  frequently  used  to 
study  nonlinear  problems  in  solid  and  fluid  mechanics.  This  particular 
thrust  of  the  research,  which  was  never  anticipated  in  the  original  research 
plan,  proved  to  be  one  of  the  most  fruitful  and  important  areas  on  investi¬ 
gation  of  the  entire  project. 

Particular  attention  was  focused  on  the  idea  of  using  exterior  penalty 
methods  and  selective  reduced  integraiton  to  handle  constrained  problems  in 
incompressible  elasticity,  incompressible  elastoplasticity,  contact  problems 
in  elastostatics,  and  constraints  in  plate  and  shell  theories.  These  methods 
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ve  been  in  wide  use  on  an  international  scale  in  1976.  It  was  discovered 
the  present  project,  however,  that  many  of  the  more  popular  methods  in 
e  can  be  dramatically  unstable  when  the  finite  element  mesh  is  refined  and, 
fact,  lead  to  divergent  approximations  of  stresses  and  pressures.  These 
suits,  it  is  felt,  had  a  significant  impact  on  this  subject  as  a  whole  and 
s  changed  the  thinking  on  the  use  of  penalty  methods  throughout  this  country 
d  abroad.  It  was  discovered  that  a  key  condition  for  the  numerical  stability 
such  methods  rests  in  the  so-called  discrete  LBB  condiditon  (Ladyszhkenskaya, 
buska,  Brezzi)  which  involves  a  stability  parameter,  ,  which  governs 
e  stability  and  convergence  characteristics  of  most  of  these  methods.  As 
ted  earlier,  it  was  discovered  that  many  of  the  more  popular  finite  element 
thods  now  in  use  do  not  satisfy  this  condition,  and,  therefore,  can  be 
stable.  The  question  then  arose  as  to  whether  or  not  elements  exist  which 
e  stable,  numerically  robust,  and  which  converge  at  optimal  rates  of 
ivergence.  In  the  latter  phases  of  the  research  effort,  several  such 
:imal,  stable,  and  convergent  methods  in  this  general  family  have  been 
.covered.  These  have  been  completely  analyzed  from  a  mathematical  point 
view  and  also  by  numerical  experiments,  A  summary  of  some  of  these  results 
given  in  Chapter  6  of  this  report, 

3.  Contact  Problems  with  Friction 

Several  classes  of  contact  problems  with  friction  have  been  analyzed 
ing  the  course  of  the  project.  Numerical  schemes  and  algorithms  have 
n  developed  for  certain  special  cases,  together  with  proofs  of  commer¬ 
ce  and  error  estimates.  In  particular,  special  problems  in  which  the 
mal  contact  pressure  is  prescribed  on  surfaces  on  which  frictional  forces 
be  developed  have  been  analyzed  in  some  detail  and  several  papers  have 
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been  written  on  this  .subject.  It  was  concluded  during  the  last  year  of 
the  project,  however,  that  many  of  the  principal  physical  and  numerical 
difficulties  encountered  in  this  class  of  problems  arose  from  the  inadequacy 
of  the  description  of  friction.  Consequently,  a  completely  new  line  of 
research  was  undertaken  to  study  new  models  of  friction.  This  is  summarized 
briefly  in  the  paragraphs  which  follow. 

4.  Non-Classical  Friction  Laws 

During  the  last  four  to  six  months  of  the  project,  considerable  at¬ 
tention  has  been  given  to  the  study  of  modelling  of  static  friction  be¬ 
tween  metallic  bodies  using  friction  laws  which  deviate  markedly  from  the 
classical  pointwise  law  proposed  by  Coulomb  200  years  ago.  In  particular, 
on  the  basis  of  micro-mechanical  mechanisms,  some  new  friction  laws  have 
been  proposed  which  feature  1)  a  nonlocal  description  of  the  contact  stress 
in  the  criteria  for  sliding  and  2)  a  nonlinear  friction  law  in  which  the 
elastic  and  elastic-plastic  response  of  metallic  junctions  on  the  contact 
surface  are  taken  into  account.  This  has  led  to  some  new  variational  prin¬ 
ciples  for  contact  problems  with  friction.  At  this  writing,  the  full  impli¬ 
cation  of  these  new  theories  has  not  yet  been  understood,  but  it  is  clear 
that  they  are  sufficiently  general  to  lead  to  a  much  better  modelling  of 
friction  effects  in  solid  and  structural  components.  A  summary  of  some  of 
these  ideas  is  given  in  chapter  5  of  this  report. 
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5.  Large  Deformation  Plasticity  and  Metal  Forming 

Considerable  effort  was  spent  on  the  development  of  finite  element 
simulators  of  large  deformations  of  elastic-plastic  materials.  A  computer 
program  was  developed  for  the  study  of  large  deformation  plasticity  prob¬ 
lems  and  metal  forming  which  is  applicable  to  a  broad  class  of  problems  in 
plane  stress,  plane  strain,  and  axisymmetric  deformation  of  bodies  of  revo¬ 
lution.  Some  impressive  numerical  results  have  been  obtained  from  this 
code,  which  apparently  surpass  all  existing  coramerical  codes  in  terms  of 
accuracy,  efficiency,  and  overall  applicability.  Nevertheless,  the  con¬ 
vergence  and  numerical  stability  characteristics  of  these  methods  are  still 
not  well  understood,  and  it  is  clear  from  some  of  the  computed  results  that 
the  choice  of  an  appropriate  friction  law  on  contact  surfaces  has  a  signif¬ 
icant  impact  on  computed  distributions  of  residual  stresses. 

6 .  Other  Areas 

Significant  research  results  in  a  number  of  other  areas  were  produced 
during  the  contract  period.  These  have  been  discussed  in  great  detail  in 
some  of  the  interim  reports  submitted  to  AFOSR  over  the  past  24  months, 
and  therefore  are  not  discussed  in  detail  here.  However,  some  results  in 
this  general  area  are  of  sufficient  importance  that  they  deserve  mention. 

We  first  note  that  new  and  useful  results  were  obtained  on  the  behavior  of 
finite- element  methods  for  singular  problems  in  structural  mechanics,  par¬ 
ticularly  for  problems  with  stress  singularities.  These  mathematical  results 
represent  the  most  general  that  have  yet  been  obtained  in  this  area 
and  provide  concrete  estimates  for  singular  problems.  The  results  focus  on 
the  behavior  of  finite  element  models  which  employ  various  singular  elements, 
a  subject  not  discussed  in  any  mathematical  literature  to  date.  It  was  shown 
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during  this  research,  for  example,  that  the  use  of  singular  elements  can 
produce  high  rates  of  convergence  for  reasonable  finite  element  meshes  even 
though  asympotically  these  methods  may  perform  poorly  or  even  diverge.  This 
type  of  phenomena  has  been  observed  in  numerical  calculations  for  a  number 
of  years,  but  there  has  not  existed  a  complete  mathematical  theory  as  to  why 
this  sort  of  behavior  is  prevalent  in  finite  element  models  which  employ 
so-called  singular  elements.  During  the  project,  this  particular  class  of 
problems  analyzed  in  great  detail;  error  estimates  were  obtained  which  were 
subsequently  confirmed  by  numerical  experiments. 

In  addition  to  these  results,  advances  were  made  in  the  modelling  of 
finite  deformation  of  elastic  bodies  including  finite  strains  of  compressible 
and  incompressible  elastic  materials.  A  very  sophisticated  finite  element 
program  was  developed  for  two-dimensional  problems  and  axisymmetric  problems 
in  finite  elasticity.  This  program  employs  a  large  library  of  finite  elements 
special  features  such  as  the  capability  of  using  an  arbitrary  strain  energy 
function  for  characterizing  materials,  provided  the  material  is  isotropic. 

The  capability  of  handling  very  large  deformations,  multi-valued  solutions, 
limit  points,  and  bifurcations,  was  also  incorporated  in  this  analysis. 

Some  results  on  nonlinear  shell  theory  were  also  recorded  during  the 
course  of  the  project.  The  first  existence  theorem  for  buckling  and  large 
deflections  of  shells  was  derived,  and  this  was  later  extended  to  include 
results  on  contact  problems  with  shells.  In  particular,  bifurcations  and 
nonlinear  buckling  phenomena  of  unilaterally  supported  shells  were  analyzed 
and  several  papers  were  published  on  this  subject.  These  papers  included 
new  theorems  on  nonlinear  eigenvalue  problems  and  bifurcations  and  are  re¬ 
garded  as  a  useful  contribution  to  the  literature  on  these  subjects. 


We  also  note  that  new  results  on  viscous  incompressible  flow  problems 
were  also  produced.  These  results  were  primarily  the  outcome  of  studies  on 
reduced  integration  penalty  methods,  and  on  the  numerical  stability  of 
various  finite  element  techniques  which  were  developed  to  study  problems 
in  finite  elasticity,  plasticity,  and  contact  problems  in  elastostatics. 


3.  SUMMARY  OF  PUBLISHED  RESULTS  AND  PROJECT  PERSONNEL 


In  this  chapter  we  summarize  the  publications  and  oral  presentations 
of  results  obtatined  during  the  contract  period. 

3.1  Journal  Articles  and  Scientific  Papers 

1.  Oden,  J.T.  and  Kikuchi,  N. ,  "Theory  of  Variational  Inequalities 
with  Applications  to  Problems  of  Flow  Through  Porous  Media,"  International 
Journal  of  Engineering  Science,  Vol.  18,  No.  10,  pp.  1173-1284,  1980. 

2.  Oden,  J.T.  and  LeTallec,  P.,  "On  the  Existence  of  Hydrostatic  Pressure 
in  Regular  Finite  Deformations  of  Incompressible  Hyperelastic  Solids," 

Nonlinear  Partial  Differential  Equations  in  Engineering  and  Applied  Science, 
Marcel  Dekker,  Inc.,  pp.  1-8.  1980. 

3.  Oden,  J.T.  and  Bernadou,  M. ,  "Theorerae  d'Existence  pour  une  Classe 
de  Problemes  Nonlineaires  de  Coques  Peu  Profondes,"  Compte-Rendus-Academy  of 
Science,  Paris,  (Analyse  Numerique) ,  pp.  1025-1028,  1980. 

4.  Oden,  J.T.,  "Penalty-Finite  Element  Approximations  of  Unilateral 
Problems  in  Elasticity,"  Approximation  Theory,  Edited  by  W.  Cheney,  Academic 
Press,  New  York,  pp.  693-697,  1980. 

5.  Oden,  J.T.  and  Kikuchi,  N. ,  "Use  of  Variational  Methods  for  the 
Analysis  of  Contact  Problems  in  Solid  Mechanics,"  Variational  Methods  in  the 
Mechanics  of  Solids,  Pergamon  Press  Ltd.,  Oxford,  pp.  259-264,  1980. 

6.  Oden,  J.T.,  Ohtake,  K. ,  and  Kikuchi,  N.,  "Analysis  of  Certain  Uni¬ 
lateral  Problems  in  Von  Karman  Plate  Theory  by  a  Penalty  Method  Part  1.  A 
Varitional  Principle  with  Penalty,"  Computer  Methods  in  Applied  Mechanics  in 
Engineering,  Vol.  24,  pp.  187-213,  1980. 

7.  Oden,  J.T.,  Ohtake,  K.,  and  Kikuchi,  N.,  "Analysis  of  Certain  Uni¬ 
lateral  Problems  in  Von  Karman  Plate  Theory  by  a  Penalty  Method  Part  2. 
Approximation  and  Numerical  Analysis,"  Computer  Methods  in  Applied  Mechanics 
and  Engineering,  Vol.  24,  pp.  317-337,  1980. 

8.  Oden,  J.T.  and  LeTallec,  P.,  "Existence  and  Characterization  of 
Hydrostatic  Pressure  in  Finite  Deformations  of  Incompressible  Elastic  Bodies," 
Journal  of  Elasticity,  Vol.  11,  No.  4,  pp.  341-357,  1981. 

9.  Oden,  J.T.  and  Bernadou,  M. ,  "An  Existence  Theorem  for  a  Class  of 
Nonlinear  Shallow  Shell  Problems,"  Journal  de  Mathematiques  Pures  et  Appl- 
lquees,  Vol.  60,  No.  3,  pp.  285-308,  1981. 

10.  Oden,  J.T.  and  Kubrusly,  R. ,  "Nonlinear  Eigenvalue  Problems  Chara¬ 
cterized  by  Variational  Inequalities  with  Applications  to  the  Postbuckling 
Analysis  of  Unilaterally  Supported  Plates,"  Journal  of  Nonlinear  Analysis, 

Vol.  5,  No.  12,  pp.  1265-1284,  1981. 
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11.  Oden,  J.T.  and  Campos,  L. ,  "Some  New  Results  on  Finite  Element 
Methods  for  Contact  Problems  with  Friction,”  New  Concepts  in  Finite  Element 
Methods,  ed.  by  M.  Spilker,  T.J.R.  Hughes,  and  D.  Gartling,  A.S.M.E.  Mono¬ 
graph,  New  York,  pp.  1-12,  1981. 

12.  Oden,  J.T.  and  Carey,  G.F.,  "Variational  Inequalities  in  Finite 
Element  Analysis"  New  Concepts  in  Finite  Element  Methods,  ed.  by  M.  Spilker, 
T.J.R.  Hughes,  and  D.  Gartling,  A.S.M.E.  Monograph,  New  York,  pp.  133-145, 
1981. 


13.  Oden,  J.T.,  "Exterior  Penalty  Methods  for  Contact  Problems  in 
Elasticity:  Nonlinear  Finite  Element  Analysis  in  Structural  Mechanics," 
Lecture  Notes  in  Mathematics.  Springer-Verlag,  Heidelberg,  pp.  655-665, 
1981. 


14.  Oden,  J.T.  and  Demkowicz,  L. ,  "On  Some  Existence  and  Uniqueness 
Results  on  Contact  Problems  with  Nonlocal  Friction,"  Journal  of  Nonlinear 
Analysis,  1982  (to  appear) 

15.  Oden,  J.T.,  Campos,  L. ,  arid  Kikuchi,  N. ,  "A  Numerical  Analysis  of 

a  Class  of  Contact  Problems  with  Friction  in  Elastostatics,"  Computer  Methods 
in  Applied  Mechanics  and  Engineering,  1982  (to  appear) 

16.  Oden,  J.T.,  "Analysis  of  a  Class  of  Contact  Problems  with  Friction 
by  Finite  Element  Methods,"  The  Mathematics  of  Finite  Elements  with  Appli¬ 
cations,  ed.  by  J.R.  Whiteman,  Academic  Press  LTD.,  London,  1982  (to  appear) 

17.  Oden,  J.T.  and  Kim,  S.J.,  "Interior  Penalty  Methods  for  Finite 
Element  Approximation  of  the  Signorini  Problem  in  Elastostatics,"  Computers 
and  Mathematics  with  Applications,  "Vol.  8,  No.l,  pp.  35-56,  1982. 

18.  Oden,  J.T.,  Song,Y.J.,  and  Kikuchi,  N. ,  "Penalty-Finite  Element 
Methods  for  the  Analysis  of  Stokesian  Flows,"  Computer  Methods  in  Applied 
Mechanics  and  Engineering,  1982,  (to  appear). 

19.  Oden,  J.T.  and  Kikuchi,  N. ,  "Finite  Element  Methods  for  Constrained 
Problems  in  Elasticity,"  International  Journal  for  Numerical  Methods  in 
Engineering,  Vol.  18,  pp.  701-725,  1982. 

20.  Oden,  J.T.,  "Analysis  of  Galerkin  Approximations  of  a  Class  of 
Pseudomonotone  Diffusion  Problems,"  SIAM  Journal  of  Mathematical  Analysis, 
Vol.  12,  No.  6,  pp.  917-930,  1982. 

21.  Oden,  J.T.  and  Pires,  E.,  "Numerical  Analysis  of  Certain  Contact 
Problems  in  Elasticity  with  Non-Classical  Friction  Laws,"  Computers  and 
Structures,  1982,  (to  appear) 

22.  Oden,  J.T.  and  Whiteman,  J.R.,  "Analysis  of  Some  Finite  Element 
Methods  for  a  Class  of  Problems  in  Elasto-Plasticity,"  International 
Journal  of  Engineering  Science,  Vol.  20,  1982. 


23.  Oden,  J.T.,  "Penalty  Methods  for  Constrained  Problems  in  Nonlinear 
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Elasticity,"  IUTAM  Symposium  on  Finite  Elasticity,  Martinum-Nijhoff  Pub., 

The  Hague,  pp.  281-300,  1982. 

24.  Oden,  J.T.,  "RIP  Methods  for  Stokesian  Flows,"  Finite  Elements  in 
Fluids,  Vol.  IV,  ed.  R.H.  Gallagher,  O.C.  Zienkiewicz,  N.  Norrie,  John  Wiley 
&  Sons,  LTD.,  London,  (to  appear) 

25.  Oden,  J.T.,  "Finite  Element  Methods  for  Constrained  Problems  in 
Elasticity  and  Fluid  Mechanics,"  Proceedings  of  Symposium  on  Finite  Element 
Methods,  Hefei,  China,  1981,  Science  Press,  (to  appear) 

26.  Oden,  J.T.,  "Mixed  Finite  Element  Approximations  via  Interior 

and  Exterior  Penalties  for  Contact  Problems  in  Elasticity,"  Hybrid  and  Mixed 
Finite  Element  Methods,  ed.  by  S.  Atluri,  John  Wiley  &  Sons  Ltd.,  London, 

(to  appear) 


27.  Oden,  J.T.  and  Pires,  E.B.,  "Algorithms  and  Numerical  Results  for 
Finite  Element  Approximations  of  Contact  Problems  with  Non-Classical  Friction 
Laws,"  Computers  and  Structures,  (to  appear) 

28.  Oden,  J.T.  and  Jacquotte,  0.,  "Stable  Second-Order  Accurate  Finite 

Element  Scheme  for  the  Analysis  of  Two-Dimensional  Incompressible  Viscous 
Flows,"  International  Conference  on  Finite  Element  Methods  in  Fluids,  Tokyo, 
Japan,  (to  appear)  ~~ 

29.  Oden,  J.T.  and  Pires,  E.,  "On  the  Analysis  of  a  Class  of  Contact 
Problems  with  Non-local  Friction,"  Proceedings  of  SECTAM  XI,  (to  appear) 

30.  Oden,  J.T.  and  Pires,  E.,  "On  the  Signorini  Problem  with  Non¬ 
local  Friction,"  Sixtieth  Anniversary  Volume  in  Honor  of  O.C.  Zienkiewicz, 
John  Wiley  and  Sons,  Ltd.,  London,  (to  appear) 

31.  Oden,  J.T.  and  Pires,  E.,  "Nonlocal  and  Nonlinear  Friciton  Laws 
and  Variational  Principles  for  Contact  Problems  in  Elasticity,"  Journal  of 
Applied  Mechanics,  (to  appear) 

32.  Oden,  J.T.  and  Pires,  E.,  "Error  Estimated  for  the  Approximations 
of  a  Class  of  Variational  Inequalities  Arising  in  Unilateral  Problems  with 
Friction,"  International  Journal  of  Numerical  Functional  Analysis  and  Op¬ 
timization,  (to  appear) 

33.  Oden,  J.T.  and  Campos,  L. ,  "Nonquasi-Convex  Problems  in  Nonlinear 
Elastostatlcs,"  (to  appear). 

3.2  Research  Reports 

1.  Oden,  J.T.  and  Bernadou,  M. ,  "An  Existence  Theorem  for  a  Class  of 
Nonlinear  Shallow  Shell  Problems,"  TICOM  Report  80-4,  Austin,  1980. 

2.  Oden,  J.T.  and  Bernadou,  M. ,  "An  Existence  Theorem  for  a  Class  of 
Nonlinear  Shallow  Shell  Problems,"  INRIA  Report  No.  17,  Le  Chesnay,  France, 
1980. 
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3.  Oden,  J.T.,  Kikuchi,  N.,  and  Song,  Y.J.,  "RIP  Methods  for  Problems 
in  Elasticity,"  TICOM  Report  80-7,  Austin,  1980. 

4.  Oden,  J.T.,  Kikuchi,  N.,  and  Song,  Y.J.,  "Reduced  Integration  and 
Exterior  Penalty  Methods  for  Finite  Element  Approximations  of  Contact  Problems 
in  Incompressible  Elasticity,"  TICOM  Report  80-2,  Austin,  1980. 

5.  Oden,  J.T.,  "RIP  Methods  for  Stokesian  Flows,"  TICOM  Report  80-11, 
Austin,  1980. 

6.  Oden,  J.T.  and  Kikuchi,  N.,  "Finite  Element  Methods  for  Constrained 
Problems  in  Elasticity,"  TICOM  Report  81-10,  Austin,  1981. 

7.  Oden,  J.T.,  Kikuchi,  N.,  and  Song,  Y.J.,  "Penalty-Finite  Element 
Methods  for  the  Analysis  of  Stokesian  Flows,"  TICOM  Report  81-11,  Austin,  1981. 

8.  Oden,  J.T.  and  Pires,  E.,  "Contact  Problems  in  Elastostatics  with 
Non-local  Friction  Laws,"  TICOM  Report  81-12,  Austin,  1981. 

9.  Oden,  J.T.  and  Demkowlcz,  L.,  "On  Some  Existence  and  Uniqueness 
Results  in  Contact  Problems  with  Nonlocal  Friction,"  TICOM  Report  81-13, 

Austin,  1981. 

10.  Oden,  J.T,  and  Pires,  E.B.,  "Nonlocal  and  Nonlinear  Friction  Laws 
and  Variational  Principles  for  Contact  Problems  in  Elasticity,"  TICOM  Report 
82-3,  Austin,  1982. 

11.  Oden,  J.T.  and  Campos,  L.T. ,  "Nonquasi-Convex  Problems  in  Nonlinear 
Elastostatics,"  TICOM  Report  82-4,  Austin,  1982. 

3.3  Theses  and  Dissertations 


The  following  Ph.D.  dissertations  and  Masters  of  Science  these  were 
completed  during  the  project. 

1.  Campos,  Luis,  "A  Numerical  Analysis  of  a  Class  of  Contact  Problems 
with  Friction  in  Elastostatics,"  M.S.  Thesis,  The  University  of  Texas  at 
Austin,  January  1981. 

2.  Kim,  S.  J.,  "Interior  Penalty  Approach  to  Contact  Problems," 

M.S.  Thesis,  The  University  of  Texas  at  Austin,  January  1981. 

3.  Song,  Y.  J. ,  "Reduced  Integration  and  Exterior  Penalty  Methods 
for  Finite  Element  Approximations  of  Contact  Problems  in  Incompressible 
Linear  Elasticity,"  Ph.D.  dissertation.  The  University  of  Texas  at  Austin, 
August  1980. 

4.  Aly,  A.,  "A  Finite  Element  Analysis  for  Problems  of  Large  Strain 
and  Large  Displacement,"  Ph.D.  dissertation.  The  University  of  Texas  at 
Austin,  May  1981, 
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5.  Kubrusly,  R.  S.,  "Variational  Methods  for  Nonlinear  Eigenvalue 
Problems  in  the  Post-Buckling  Analysis  of  Unilaterally  Constrained  Elastic 
Structures,"  Ph.D.  dissertation.  The  Unviersity  of  Texas  at  Austin,  August  1981. 

6.  O'Leary,  J.  T.,  "An  Error  Analysis  for  Singular  Finite  Elements," 

Ph.D.  dissertation.  The  University  of  Texas  at  Austin,  August  1981. 

3.4  Oral  Presentations 

Over  50  oral  presentations  were  given  by  members  of  the  research 
team  during  the  contract  period  on  research  results.  A  partial  list  of 
these  is  given  as  follows: 

1.  Oden,  J.T. , "Penalty  Methods  and  Selective  Reduced  Integration  for 
Stokesian  Flows,"  Third  International  Conference  on  Finite  Elements  in  Flow 
Problems,  Banff  Centre,  Banff,  Alberta,  Canada,  June  10-13,  1980. 

2.  Oden,  J.T.,  "Penalty  Methods  and  Reduced  Integration  in  Elasti¬ 
city  Problems,"  Symposium  on  Finite  Element  Methods  for  Nonlinear  and  Singu¬ 
lar  Problems,  University  of  Durham,  Durham,  England,  June  26- July  6,  1980. 

3.  Oden,  J.T.,  "Analysis  of  Incompressible  Elastic  Bodies  by  Exten¬ 
sive  Penalty  Methods,"  U.S. /Europe  Workshop  on  Nonlinear  Finite  Element 
Analysis  in  Structural  Mechanics,  Ruhr  University,  Bochum,  West  Germany, 

July  28-30,  1980. 

4.  Oden,  J.T.,  "Penalty  Methods  for  Constrained  Problems  in  Non¬ 
linear  Elasticity,"  IUTAM  Symposium  on  Finite  Elasticity,  LeHigh  University, 
Bethlehem,  PA,  August  11-15,  1980. 

5.  Kikuchi,  N.,  "Penalty  Finite  Element  Approximations  of  a  Class  of 
Contact  Problems  in  Linear  Elasticity,"  Presentation  made  to  Exxon  Products 
Research  Laboratories,  Houston,  Texas,  August  14,  1980. 

6.  Oden,  J.T.,  "Penalty  Methods  for  Stokesian  Flows,"  Engineering 
Mechanics  Seminar,  The  University  of  Texas  at  Austin,  September  23  and  25, 

1980. 


7.  Oden,  J.T.,  "Constrained  Problems  in  Nonlinear  Elasticity,"  Annual 
meeting  of  the  Society  of  Engineering  Science,  Atlanta,  Georgia,  December,  1980. 

8.  Oden,  J.T.,  "Contact  Problems  with  Friction  in  Elastostatics," 

Annual  meeting  of  the  Society  of  Engineering  Science,  Atlanta,  Georgia, 

December,  1980. 

9.  Oden,  J.T.,  "Finite  Element  Methods  for  Contact  Problems  with 
Friction,"  Department  of  Mathematics,  University  of  Maryland,  February  25,  1981. 

10.  Oden,  J.T.,  "Penalty  Methods  for  Stokesian  Flows,"  Mathematics 
Seminar,  University  of  Texas  at  Austin,  March  23,  1981. 
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11.  Oden,  J.T.,  "Mixed  Finite  Element  Approximations  via  Exterior 
and  Interior  Penalty  Formulations  of  Contact  Problems  in  Elasticity,"  Sym¬ 
posium  on  Hybrid  and  Mixed  Finite  Element  Methods,  Georgia  Institute  of 
Technology,  Atlanta,  Georgia,  April  10,  1981. 

12.  Oden,  J.T.,  "Finite  Element  Methods  for  Contact  Problems  with  Friction," 
MAFELAP  '81,  Fourth  Conference  on  the  Mathematics  of  Finite  Elements  with 
Applications,  Brunei  University,  Uxbridge,  England,  April,  1981, 

13.  Oden,  J.T.,  "Non-Local  Friction  Laws  in  Elastostatics,"  Department 
of  Civil  Engineering,  The  University  of  Swansea,  Wales,  United  Kingdom,  April 
1981. 


14.  Oden,  J.T.,  "An  Introduction  to  the  Exterior  Penalty  Methods  and 
Reduced  Integration  and  Their  Application  to  Constrained  Problems  in  Elas¬ 
ticity  and  Fluid  Mechanics,"  Mathematics  Seminar,  Brunei  University,  Uxbridge, 
England,  May  13,  1981. 

15.  Oden,  J.T.,  "Penalty-Finite  Element  Methods  for  Constrained  Problems 
in  Elasticity,"  Invitational  Symposium  on  Finite  Element  Methods,  Polytechnical 
University,  Hefei,  People’s  Republic  of  China,  May  19-23,  1981. 

16.  Oden,  J.T.,  "Some  New  Results  on  Finite  Element  Methods  for  Contact 
Problems  with  Friction,"  Symposium  on  New  Concepts  in  Finite  Element  Methods, 
Joint  ASME/ASCE  Mechanics  Conference,  The  University  of  Colorado,  Boulder, 

June  23,  1981. 

17.  Oden,  J.T.,  "Variational  Inequalities  in  Finite  Element  Analysis," 
Symposium  on  New  Concepts  in  Finite  Element  Methods,  Joint  ASME/ASCE  Mechanics 
Conference,  The  University  of  Colorado,  Boulder,  June  23,  1981. 

18.  Oden,  J.T.,  "Numerical  Analysis  of  a  Class  of  Contact  Problems  with 
Friction  in  Elastostatics,"  FENOMECH  1981,  Stuttgart,  West  Germany,  August 
26-28,  1981. 

19.  Oden,  J.T.,  "Nonconyex  Problems  and  Phase  Transition  in  Nonlinear 
Materials,"  SECTAM,  Huntsville,  Alabama,  April,  1982. 

20.  Oden,  J.T.,  "Mathematical  Theory  of  Plasticity,"  EM/TICOM  Seminar, 

The  University  of  Texas  at  Austin,  April  1982. 

21.  Oden,  J.T.,  "A  Stable  Second-Order  Accurate  Finite  Element  Scheme 
for  the  Analysis  of  Two-Dimensional  Incompressible  Viscous  Flows,"  4th 
International  Symposium  of  Finite  Element  Methods  in  Flow  Problems,  Tokyo, 
Japan,  July  27,  1982. 


PROJECT  PERSONNEL 


The  following  individuals  worked  on  the  project  during  the  contract 
period: 

1)  Principal  Research  Staff: 

Dr.  J.  Tinsley  Oden,  Principal  Investigator 
Dr.  E.  B.  Becker,  Co-Principal  Inviestigator 
Dr.  N.  Kikuchi,  Senior  Research  Associate 
Dr.  L.  Demkowicz ,  Senior  Research  Associate 

2)  Senior  Research  Assistants: 

J.  O’Leary  L.  Campos 

Y.J.  Song  E.  Pires 

S.J.  Kim  R.  Kubrusly 

Messrs.  T.  Pan,  Jerry  Fine,  R.  Chambers  and  P.  Halamek  devoted  a  very 
small  percentage  of  their  time  to  certain  phases  of  the  project. 

3)  Project  Secretary: 

The  project  secretary  has  been  Mrs.  Ruth  Dye  during  the  past  year. 
During  the  first  year  of  the  project,  Ms.  Nancy  Webster  served  as  project 
secretary.  Some  additional  part-time  secretarial  assistance  was  required 
occassionally  during  the  project. 


4.  THEORY  AND  APPROXIMATION  OF 


CONTACT  PROBLEMS  IN  ELASTOSTATICS 


1.  Introduction 

Perhaps  the  most  primitive  and  intrinsic  feature  of  the  mechanics  of  solids 
is  the  contact  of  one  body  with  another.  Contact,  in  fact,  is  precisely  the 
physical  event  through  which  loads  are  delivered  to  a  structure  and  by  which  a 
structure  transmits  forces  to  its  supports.  Nevertheless,  this  fundamentally 
important  aspect  of  structural  behavior  has,  until  recently,  rarely  been  taken 
into  account  in  practical  structural  analysis  and  design.  Ihe  underlying 
difficulty  is  that  contact  problems  in  solid  mechanics  are  inherently  nonlinear: 
the  area  of  contact  is  not  known  prior  to  the  application  of  loads  and  complex 
physical  phenomena  are  experienced  on  the  contact  surfaces  which  often  require 
special  mechanical  and  mathematical  considerations. 

In  recent  years,  however,  significant  advances  have  been  made  in  the  study 
of  certain  restricted  classes  of  contact  problems  by  finite  element  methods. 
These  classes  of  problems  include  those  adequately  modelled  by  the  so-called 
Signorini  problem  in  elastostatics :  the  behavior  of  a  linearly  elastic  body  in 
unilateral  contact  with  a  rigid  frictionless  foundation.  This  particular  class 
of  problems  can  be  studied  within  the  framework  of  the  theory  of  variational 
inequalities  and;  consequently,  a  great  deal  can  be  established  on  the 
qualitative  behavior  of  the  solutions  and  their  finite  element  approximations. 
In  fact,  there  are  few  nonlinear  problems  in  structural  mechanics  for  which  a 
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i  present  chapter  is  devoted  to  a  general  exposition  on  finite  element 
for  contact  problems  in  elastostatics  that  can  be  characterized  by 
>nal  inequalities.  In  keeping  with  the  general  objectives  of  this 
,  this  class  of  problems  is  chosen  for  consideration  because  of  it6 
>1  importance,  the  richness  of  the  mathematical  foundations  on  which  it 
>ed,  the  level  of  results  one  can  establish  on  the  behavior 
>nvergence  and  numerical  stability)  of  corresponding  finite  element 
nations,  and  because  this  class  of  numerical  problems  provides  the 
lity  to  examine  many  new  algorithms  for  treating  systems  of  linear  and 
ir  inequalities  which  are  drawn  from  linear  and  nonlinear  programming 
ind  optimization. 

i  account  of  the  subject  which  follows  deals  primarily  with  finite 
approximations  of  the  Signorini  problem  for  linearly  elastic  bodies 
lg,  for  the  sake  of  completeness,  formulations  in  which  friction  Is  taken 
:ount.  A  lengthy  treatise  on  this  and  related  subjects  has  been  recently 
*d  by  Kikuchi  and  Oden  [1983],  and  this  work  can  be  consulted  for  a 
list  of  references  on  the  subject  and  for  more  details.  Following  this 
.tory  section,  the  system  of  differential  equations  and  boundary 
ins  for  Signorini's  problem  with  friction  is  derived.  In  Section  4.3, 
variational  formulations  of  contact  problems  are  derived,  all  of  which 
variational  inequalities  given  in  terms  of  an  unknown  equilibrium 
ment  field.  Existence,  uniqueness,  and  approximation  results  for  the 
which  tangential  contact  pressures  are  prescribed  is  taken  up  in  Section 
ese  problems  include  as  a  special  case  the  classical  Signorini  problem 
friction.  Section  4.5  is  devoted  to  contact  problems  in  which  the 
ontact  pressure  is  prescribed.  These  classes  of  problems  are  shown  to 
ndamental  Importance  in  the  analysis  of  problems  with  non-local  friction 


in  Chapter  5.  The  results  of  several  numerical  experiments  are  also  discussed  in 
Sections  4.4  and  4.5;  concluding  comments  are  collected  in  a  final  section. 

2.  Signorini's  Problem  with  Fr iction 

Consider  a  deformable  body  supported  unilaterally  by  a  rigid  body  as  shown 
in  Figure  4.1.,  and  suppose  that  the  contact  surface  of  the  two  bodies  is 
unknown  a  priori.  The  problem  of  finding  the  deformation  and  contact  force  for 
equilibrium  configurations  of  the  body  under  certain  boundary  and  loading 
conditions  is  called  Signorini's  problem.  We  note  that  in  most  works  in 
theoretical  mechanics  Signorini's  problem  is  discussed  in  regard  to  the  special 
case  in  which  no  friction  exists  on  the  contact  surface.  We  first  formulate 
Signorini's  problem  with  friction  for  the  case  of  small  deformations  in  which 
the  displacement  field  is  small  enough  so  that  higher  order  terms  of  the 
displacement  gradient  can  be  neglected  in  the  equilibrium  and 
strain-displacement  equations.  The  stress  tensor  a  =  (8)  ij*  then  the 

Cauchy  stress  tensor  and  no  distinction  between  particles  or  points  is  necessary 
in  defining  its  domain. 

4.2.1.  Equilibrium  Equations 

Let  u  =  u^i^  be  the  displacement  field  that  provides  equilibrium  of  the 
deformable  body  under  given  loading  and  boundary  conditions,  and  let 
e  =  ^-ij^i  ®  denote  the  strain  tensor  for  small  deformation  defined,  for  an 
arbitrary  displacement  v,  by 

*The  summation  convention  is  applied  throughout  this  chapter,  unless 
specifically  noted  otherwise. 


Figure  4.1.  Geometry  of  an  elastic  body  near  a  rigid  foundation 


where 


eij(i>  “  <vi,j  +  vj,i>/2 


v^  j  denotes  the  partial  derivative  of  v^  in  the 
dvj/dx..  Suppose  that  the  stress-strain  relation 


coordinat 


vi.i 
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ft 


ft 


ft 


ft 


°lj  =  EijkiekJl 


(4.2.2) 


is  given,  wherein  the  elasticities  E^^  of  the  material  satisfy  the  conditions 


EijkJl  =  EkJlij  =  Ejik* 


(4.2.3) 


nc  ^  Ei jkJl^i^kJl^i j  ^  rac*i j^i j  *  2  ^ij  —  ^ji  (4.2.4) 


Max  l^^jlc^(x)|  <  M  for  all  i,j,k,  and  A. 
x  Q 


(4.2.5) 


Here  Q  is  an  open-bounded  domain  in  RJ  representing  the  interior  of  the 
deformable  body. 

The  condition  =  Ekjj.ij  guarantees  the  existence  of  the  strain  energy 

function  Eq  such  that 


°ij  “  aVaeij> 


E0  =  J  Eijk^ekJleij* 


(4.2.6) 


(4.2.7) 


Condition  (4.2.4)  is  sufficient  to  imply  strict  convexity  otf  the  strain 
energy  function  Eq,  and  guarantees  the  coercivity  of  the  total  potential  energy 
of  the  body,  as  shown  later.  Condition  (4.2.5)  merely  establishes  boundedness 
of  the  elasticities  these  material  parameters  may  be  discontinuous  in  the 


domain  Q 
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If  o(u)  denotes  the  stress  field  at  the  equilibrium  configuration  given  by 
the  displacement  u,  the  principles  of  balance  of  linear  momentum  and  moment  of 
linear  momentum  imply  the  equilibrium  equations  and  symmetry  of  stress  tensor: 


dzu< 


dt 


r  -  “'shj.j  ■  h 


in  Q 


°ij  =  °ji  J 


(4.2.8) 


Here  p  is  the  density  of  the  body,  f  =  f^i^  is  the  applied  body  force  per  unit 
volume. 


4.2.2.  Boundary  Conditions 

Let  the  boundary  T  of  the  deformable  body  be  divided  into  three  disjoint 
parts  and  Tc  for  each  index  1  <  i  <  3.  The  i  component  of  the 

displacement  is  prescribed  on  the  boundary  rD^,  and  the  i  component  of  the 
traction  is  given  on  1^.  The  part  is  identified  with  the  actual  contact 
surface  and  must  be  large  enough  to  include  the  true  contact  surface  that,  at 
this  point,  is  unknown.  For  simplicity,  we  assume  that  there  are  no  other 
forces  except  those  due  to  contact  on  Tq. 

If  t  *  iiii  denotes  the  applied  traction,  the  two  standard  boundary 
conditions  on  and  Tp^  are 


Uf  -  u*  on  rM 


Fi 


,  1  <  i  <  3, 


(4.2.9) 


o(u)ijnJ  =  tt 


on  r 
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where  n  ■  n^ij  is  the  unit  vector  outward  normal  to  the  boundary  T. 

Suppose  that  the  boundary  is  separated  into  two  parts,  and  as 

shown  in  Figure  4.1,  On  rj.,  the  body  Q  is  in  contact  only  with  the  rigid 
foundation,  whereas  the  interaction  of  two  opposite  surfaces  must  be  considered 
on  Furthermore,  let  a  crack-like  slit  exist  inside  the  domain  Q,  and  let 

its  surface  be  represented  by  T®.  Let  us  construct  a  coordinate  system  0a, 
1  <  a  <  3,  on  each  part  of  the  boundary  related  to  contact,  such  that  two 
opposite  surfaces  can  be  distinguished  from  each  other.  Examples  of  such 

coordinate  systems  are  given  in  Figure  4.1.  Let  ua  be  the  a  component  of  the 
displacement  u  in  the  9a  coordinate  system,  i.e.,  It  uQ  *  u»ia.  Then  the 

kinematical  contact  conditions  due  to  contact  are  given  by 

kapup  “  8a  <  0  on  rC  (4.2.10) 

kipup  +  ka(3up  -  8a  <  0  on  rC  and  rC  (4.2.11) 

where  kfflp,  k*p,  1  <  i  <  2,  and  ga  depend  on  0a,  1  <  a  <  3,  and  the  superscripts 
1  and  2  on  kap  and  Up  indicate  the  two  opposite  surfaces  on  which  these 
quantities  are  defined.  The  functions  kap  and  ga  must  be  obtained  from  the 
geometry  of  the  body  and  the  rigid  foundation. 

For  example,  let  us  consider  the  part  Tq,  and  obtain  the  corresponding 

expressions  for  kap  and  ga.  To  do  this,  suppose  that  the  surfaces  of  the  body 

and  the  rigid  foundation  are  given  by 

0p  =*  $(0jj02)»  =  4*( ® i » 9 2^ »  (4.2.12) 

respectively.  Since  the  body  cannot  penetrate  the  rigid  surface, 
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®3  +  u3  >  4»(0j  +  Up@2  +  U2), 


(4.2.13) 


where  u  *  “aia*  Linearization  of  the  function  <J>  in  and  ^  yields 


03  +  u3  >  4.(©1,02)  +  Va<|>(Oi,02Wuiii  +  u2i2^ 


where  =  ijd/dOj  +  123/S02*  Dividing  by  /I  +  l^a4»H  and  noting  that 


^^■7rr^[^u+^-h) 


(4.2.14) 


is  the  unit  vector  inward  to  the  rigid  surface  ©3  =  <Ji(0  2,02),  we  have  the 


linearized  contact  condition 


NpUp  -  g  <  0, 


(4.2.15) 


where 


g  -  - 1-  -  -  (♦  -  *)N3. 

✓  1  +  |7_<HZ 


(4.2.16) 


Thus,  if 


kop  “  6«pNp*  8l  “  8*  82  "  83 


(4.2.17) 


in  the  general  form  (4.2.10),  the  condition  (4.2.10)  represents  the  linearized 
contact  condition  (4.2.15)  on  Tq.  Since  N  is  the  same  as  the  unit  outward 
normal  to  the  surface  of  the  deformable  body  on  the  contact  surface  after  the 
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That  is,  the  deformation  of  the  point  JP  is  always  on  or  above  the  tangent  plane 
A  at  the  point  Q  of  the  rigid  surface  &2  *  • 

In  almost  all  literature  on  contact  problems  dealing  with  theoretical 
aspects  of  these  problems  the  linearized  form  (4.2.15)  of  (4.2.13)  is  not  used. 
The  condition  more  commonly  applied  is 

n0u0  “  8  <  0,  (4.2.18) 

1  * 

where  n  =  npip  is  the  unit  vector  outward  normal  to  the  boundary  Ig»  and  g  is 

A 

the  projection  of  (4i  -  $)  in  the  n^-direction;  g  =  (4#  —  ^n^.  The  linearized 
condition  (4.2.18)  is  represented  by  the  plane  B  in  Figure  4.2.  In  this  case, 
the  deformation  of  the  point  £  is  always  on  or  above  the  plane  B.  Of  course, 
since  we  are  dealing  with  infinitesimal  deformations,  the  contact  surface  must 
be  very  close  to  the  foundation,  and  the  differences  between  n  and  N  are, 

A*  » 

accordingly,  small.  We  shall  apply  (4.2.15)  in  this  chapter. 

Following  a  similar  approach,  a  linearized  contact  condition  on  T®  and 
can  be  obtained.  Let  the  particles,  the  position  vectors  of  which  at  the 
initial  stage  are  0^  ■  0*ia  an<*  6^  “  ®a*a  w^ere  ®3  *  an<*  ®3  =  ’K®!*®^ 

be  in  contact  after  deformation.  That  is,  let 

+  ui  "  °a  +  ua *  a  =  l>2  (4.2.19) 

ii  9  ? 

where  u1  =  ua*a  an^  u  =  uo^a  are  t*ie  displacement  vectors  of  the  particles 
identified  by  9*  and  0^.  Here,  the  two  opposite  surfaces  are  assumed  to  be 
represented  by  the  equations 

9^  3  <ind  9^  3  (4.2.20) 

for  a  proper  coordinate  system 
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The  "no  penetration"  condition  is  then  given  by 

®3  +  u3  >  e3  +  u3  (4.2.21) 

Using  (4.2.19)  we  now  express  the  condition  (4.2.21)  In  terms  of  0*.  Neglecting 
higher  order  terms,  we  have 

♦(e},e|)  +  u$(e},e|.4(e{,ej))  >  *<e{,e£) 

+  ^L(e},eJ)(  u}(011,ej,«(e{,e^))  -  ) 

+  ^(0l1,02)(  ^(011,0j,4,(011,0j))  -  u{(011,0j,<|;(011,021))  ) 

+  u§(011,0j,<p(011,0])).  (4.2.22) 

For  simplicity,  we  shall  express  (4.2.22)  as 

♦  +  u3  >  *  +  ^  <u!  -  ul>  +  ^  <u2  -  u2>  +  u3 

Applying  the  relation  (4.2.14)  of  the  unit  vector  inward  normal  to  the  surface 
03  *  <K0i*02)»  we  can  represent  (4.2.22)  as 

VUB  "  UB^  "  «i  <  0, 


(4.2.23) 
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gj  ■  (4»  “  ♦)n3  • 


(4.2.24) 


Thus,  if 

kip  *  -  klp  *  V  6l  -  8 

ul  «  V  2 

op  op  =  6ap  for  a  -  2,3,  g2  -  g3  -  + 


(4.2.25) 


then  the  linearized  condition  (4.2.23)  is  given  by  the  generalized  form 
(4.2.11). 

As  far  as  the  linearized  forms  are  concerned,  the  contact  conditions 

(4.2.15)  and  (4.2.23)  have  a  similar  form.  More  precisely,  if  we  replace  the 

displacement  u  in  (4.2.15)  by  the  relative  displacement  uR  =*  u*  -  u^,  the 

fl  7 

condition  (4.2.15)  can  represent  (4.2.23)  for  the  boundaries  and  Iq. 

Therefore,  for  the  mathematical  development  of  contact  problems,  we  need  only 
consider  the  conditions  on  boundary  r£.  We  thus  describe  the  general  kinematic 
contact  condition  as 


kapup  ~  «a  <  0  on  rC» 


(4.2.26) 


or  more  precisely 


NpUp  -  g  <  0  on  Tq. 


(4.2.27) 


These  represent  only  kinematical  requirements.  In  addition,  the  stresses  must 
satisfy 
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Ojj  *  o»N  <  0  if  there  is  contact 

(A. 2. 28) 

Ojj  *  0  if  there  is  no  contact 

where  o  »  (°ijnj);Li  Is  the  traction  on  the  boundary  Tq.  Combining  (A. 2. 27)  and 
(A. 2. 28)  yields 

Ujj  -  g  <  0,  Ojj  <  0,  °n^uN  —  =  0»  (A. 2. 29) 


where  Ujj  *  N^u^  =  u»N. 

If  the  contact  surface  is  well  lubricated,  friction  might  be  neglected.  In 
this  case,  on  the  boundary  T^.,  the  additional  boundary  condition, 

aT  =■  a  -  aN?J  =  0  (A. 2. 30) 


must  be  applied,  where  a  is  again  the  traction.  However,  in  many  cases  In 
practice,  friction  has  an  important  role  to  stress  analysis.  Here  we  shall 
simply  assume  that  Coulomb's  friction  law  holds  pointwise  on  the  contact 
surface.  More  general  friction  laws  have  been  considered  recently  by  Duvaut 
[1981],  and  Oden  and  Pires  [1981]  and  we  explore  one  such  class  of  friction 
problems  In  the  next  chapter.  The  tangential  velocity  is  assumed  to  be  governed 
by  the  following  form  of  Coulomb's  law: 


I  2t I  <  "  then  uT  *  0 


|oTl  *  -  pOjj,  then  for  some  \  >  0,  u-j.  *  -  Xa^ 


(A. 2. 31) 
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These  conditions  hold  on  the  surface  in  which  cr^  <  0.  Here  u  =  du/dt  is  the 
velocity  and  u^.  is  its  tangential  component  defined  by 

•  •  •  •  • 

uT  =  u  -  u^N,  Ujq  =  u  *  N. 


4.2.3.  Boundary-Value  Problems 

In  summary,  we  have  the  following  initial-boundary  value  problem 


pu  -  div  c(u)  *  f  in  Q 


o(u)  =  t  on  I\ 


u  =  u  on  rT 


1 2x1  <  “  then  “t  =  ® 


1 2x1  =  “  then  X  >  0  s.t.  uT  = 


with  the  initial  conditions 


>  (4.2.32) 


uN  -  g  <  0,  oN  <  0,  ctn(un  -  g)  -  0  ^ 


r  °n  rC 


-  XOrj,  / 


J 


u  ■  Uq  and  u  ■  Vq  at  t  *  0, 


(4.2.33) 


for  the  Signorini  problem.  Here  a  *  0£ji^0ij,  and  div  is  the  divergence 


operator 
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s 

1 


( 


If  the  motion  is  slow,  then  the  aceelerat ion  is  small  and  the  process  of 
the  deformation  of  the  body  represented  by  (4.2.32)  and  (4.2.33)  can  be 
approximated  by  the  following  incremental  form  using  an  artificial  monotonically 
increasing  time  parameter  t.  Let  uC  be  the  deformation  at  time  t,  and  let  Au  be 
a  small  increment  of  u  produced  during  the  time  interval  At,  i.e., 
Au  *  ut+^t  -  uc.  Suppose  that 


-  div  o^1-)  =  f11  in  Q 


o(ut)  «*  tc  on  r; 


(4.2.34) 


ut  =  uC  on  Tq, 


where  f1",  t1-,  and  uc  are  the  values  of  f,  t,  and  u  at  time  t,  respectively, 


Introducing  similar  quantities  Af,  At,  and  Au  to  Au,  we  arrive  at  the  boundary 
value  problem 


"S 


-  div  o(Au)  =  Af  in  Q 


a(  Au)  =  At  on  fp 


Au  =  Au  on  T, 


>(4.2.33) 


AuN  -  Ag  <  0,  AaN  +  oN  <  0,  (AoN  +  cN)(AuN  -  Ag)  =  0 


|AOj.  +  d~\  <  -  p(AoN  +  5«) ,  then  AuT  =  0 
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|Aot  +  oTl  =  -  p(AoN  +  aN),  then  X  >  0  s.t. 

Auij>  ~  (J^)  i 

the  time  interval  [t,t+At].  Here 
Ag  =  g  -  u^,  Aon  =  aN(Au),  oN  =  Oj^u1) 


ot(Au) , 


i 


°n  rc 


J 


(4.2.36) 


eed  not  be  same  as  X  in  (4.2.32). 

refer  to  problem  (4.2.32),  (4.2.33)  as  the  dynamical  form  of  Signorini's 
with  Coulomb  friction.  The  system  (4.2.35),  thus,  defines  an 
ntal  form  of  this  problem  valid  for  slow  deformation  or  quasi- 
^nditions. 


3.  Variational  formulations  of  Signorini's  Problem 
shown  in  the  previous  section,  Signorini's  problem  with  friction  for 

i  in  elastostatics  assumes  the  form  of  the  boundary  value  problem 
.  We  shall  consider  a  variational  formulation  of  the  problem  that 
i  the  case  (4.2.35)  as  a  special  case,  and  shall  develop  a  mathematical 
»f  such  problems  subject  to  a  few  additional  hypotheses. 


4.3.1.  Abstract  Signorinl's  Problem  with  Friction 


Let  us  consider  the  boundary  value  problem  of  finding  a  displacement  field 
u  such  that 


-  div  o(u)  =  f  in  Q 


o(u)  =  t  on  rF 


u  =  u  on  r  r 


>(4.3.1) 


UN  “  8  ^  "b  Ojj  <  0,  ( cr^j  +  —  g)  -  0 


|oT  +  oTl  <  -  p(oN  +  <*N),  then  uT  =  0 


|aT  +  aTl  =  -  p(oN  +  aN),  then  X  >  0,  s.t.  uT  =  -  X(aT  +  aT) 


J 


for  given  functions  f,  t,  u,  cfy,  a^,  and  g. 

Multiplying  the  first  member  of  (4.3.1)  by  an  arbitrary  smooth  test 
function  v  such  that  v  =  0  on  Tp,  and  integrating  by  parts,  we  have 


Jq  a(y>ijvi,jdx  "  Jra(ii>ijnjvids  *  / Q  £ividx* 


Using  the  relations 
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a.b  =  aNbN  4  aT*bT,  aN  =  a*N,  aT  =  a  -  a^N,  etc. 

for  the  unit  vector  N,  and  applying  the  boundary  conditions  (4.3. 1)2*  and 
(4.3. 1)3,  we  obtain 

Jq  o(u)ijvi>jdx  =  Jq  fiVldx  4  Jr  tiVlds 


4  /p  (o(u)NvN  4  a(u)j*v^)ds,  V  v  s.t.  v  =  0  on  T 

c 


D* 


(4.3.2) 


where  vN  =  v«N  =  vaNa,  va  =  vT  =  v  -  vT,  a(u)^  =  a(u)  *N,  and  so  on. 

Let  w  be  an  arbitrary  function  that  satisfies 


w  =  u  on  rD,  wN  -  g  <  0  on  Tq. 


(4.3.3) 


It  is  noted  that  (4. 3. 1) ^-(4 . 3. 1)^  are  equivalent  to  the  two  inequalities 
o(^N(wN  "  UN>  =  +  5nHwN  ~  8  “  uN  +  s) 

-  °n(wN  "  UN^  >  ~  5N^wN  _  UN^» 


and 


o(u)T*(wT  -  uT)  >  n(a(u)N  4  on)(|wt|  -  |uTl)  -  §T» (wT  -  uT) 
Thus,  (4.3.2)  can  be  written  as 

JQ  °(H)ij(wi,j  -  ui,j)d*  +  hc  +  5n)(,2ti  “  l-T1 ) ^ds 
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>  Iq  ii(vl  -  uA)dx  +  /pF  t1(wi  -  Ui)ds 
-  / rc  5i(wi  “  ui)ds» 

for  every  w  satisfying  (4.3.3),  where  the  traction  vector  a  is  defined  by 
o  =  onN  +  oT  on  rc. 

Let  us  introduce  the  following  standard  notation: 

B(u,v)  =  /q  cr(u)i jeij^)dx 

f(v)  =  Jq  f^v^dx  +  /p  t^Vpds  +  / p  o^v^ds  (4.3.4) 

j(u;v)  =  / r  {-  p(o(u)N  +  oN)}|vT|ds, 

Then  the  problem  (4.3.1)  yields  a  variational  form 

B(u,w-u)  +  j(u;w)  -  j(u;u)  >  f (w-u) ,  (4.3.5) 

for  every  w  satisfying  (4.3.3),  (since  the  stress  tensor  a  is  symmetric). 

Using  (4.2.3)  and  (4.2.5),  we  have 

B(u,v)  -  Jq  EijkJluk,Avi,jdx 

<  MI ul jlvl j  (4.3.6) 


where 


1  V>  1  *  t Jq  +  viv1)dx} 


(4.3.7) 


If  the  given  data  f,  t,  and  o  are  smooth  enough,  i.e.,  if 

f±  6  L2(Q),  tt  £  L2(rF),  5N  €  L2(rc),  aTi  €  L2(rc), 


(4.3.8) 


then 


f(v)  <  (lflQ  + 


'-"o.rp 


(4.3.9) 


where 


l£»0  =  {J0  f^dx}1'*, 

(4.3.10) 

"i'o.Tp  =  {/rF  titidsJ 

and  lol0>rc  *s  defined  similarly  to  (4.3.10)2*  Thus  the  bilinear  form  B(u,v) 
represents  the  virtual  work  of  the  elastic  body,  and  the  linear  functional  f(v) 
of  the  work  done  by  the  virtual  displacement  v.  Both  are  well-defined  on  the 
Sobolev  space  H*(Q): 

H1^)  «  {v  -  vii1  :  vt  €  L2(Q),  Vljj  6  L2(Q),  1  <  i,j  <  3),  (4.3.11) 


where  v^  j  is  the  generalized  derivative  of  v^.  However,  the  term  j(u;v)  is  not 
well-defined,  since  the  normal  stress  o(u)^  cannot  be  defined  for  a  function 
u  £  H^(Q).  The  contact  stress,  as  It  enters  the  present  formulations,  cannot  be 
defined  in  the  appropriate  space  because  we  cannot  infer  sufficient  smothness  of 
the  displacement  u.  Thus,  in  this  sense,  the  variational  formulation  (4.3.5)  of 
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the  Signor Ini  problem  with  f r iction  may  not  be  meaningful  in  the  sense  of  this 
variational  principle. 

For  a  particular  case  that  the  friction  coefficient  p  is  sufficiently 
small,  there  are  several  approaches  one  can  use  to  show  the  existence  of 
solutions  to  the  variational  formulation  (4.3.5)  within  the  context  of  the 
Sobolev  space  H*(Q).  Details  of  such  treatments  can  be  found  in  e.g.,  Necas, 
Haslinger,  and  Tuisuk  [1980],  Duvaut  [1981],  Oden  and  Demkowicz  [1981],  Oden  and 
Plres  [1981],  and  others.  In  the  present  study,  we  shall  not  discuss  these 
mathematical  technicalities,  rather,  we  shall  address  the  problem  of  obtaining  a 
solution  of  (4.3.5)  using  finite  element  methods.  To  do  this  the  following  two 
partlclar  cases  of  (4.3.5)  will  be  studied  in  detail. 

4.3.2.  Special  Case  I:  Prescribed  Tangential  Stress 

Suppose  that  the  tangential  stress  is  known.  That  is,  the  tangential 
stress 


o*j*  —  t  «jt  on  r*£ 


(4.3.12) 


is  prescribed,  and  the  Coulomb  friction  law  is  abandoned  on  the  contact  surface 
rc.  Then  the  variational  formulation  (4.3.5)  reduces  to 


B(u,w-u)  >  fj(w-u). 


(4.3.13) 


for  every  w  satisfying  (4.3.3),  where 


f,(v)  *  f(v)  +  /p  £T»vTds. 


(4.3.14) 
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tT1  €  L2<rc),  1  <  i  <  3,  (4.3.15) 

the  linear  fora  fj  is  bounded  in  H*(Q),  i.e. , 

<  i'l'o  +  l£lo,rF  +  1 21  o, +  ,£Tlo,r(J,2,i*  (4.3.16) 

In  this  case,  there  are  no  ambiguous  terms  such  as  j(u;v)  in  (4.3.5)  in  the 
variational  formulation  (4.3.13). 

4.3.3.  Special  Case  lit  Prescribed  Normal  Stress 

Suppose  that  the  normal  stress  oN  is  known  at  this  time.  That  is,  the 
normal  stress 

oN  =  tN  €  L2(rc)  (4.3.17) 

is  prescribed,  and  only  the  Coulomb  law  of  friction  has  to  be  considered  on  the 
contact  surface  Tc,  while  the  unilateral  contact  condition  (4.3. 1)^  is 
abandoned.  Then  the  variational  formulation  (4.3.5)  becomes 

B(u,w-u)  +  j(w)  -  j(u)  >  f2(w-u) ,  (4.3.18) 


for  every  w  such  that 

wj  -  0  on  rD^,  1  <  1  <  3, 


(4.3.19) 


where 


J(v)  “  hr  ^(tN  +  5N)Jl£T,d8* 

v 


(4.3.20) 


V* ~'W  — T  ■»  \  » v  »  "  '  ■»  .  V^-  .  , •: — ■'.-'"V"  's-vi  i  m  w  -m--  »  ■■  i—  - - - -  »  -  w  — 
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-W 

and 

n 

•» .  > 

“  f(v)  +  Jr-  tNvNds* 

*°  V 

(4.3.21) 

J 

■ 

.  o 

Since  both  a jj  and  tjj  are  In  L  (T^),  the  convex  functional  j  Is  bounded 

in  H^Q), 

• 

and  Is  continuous.  Indeed, 

-i 

1 

|j(v)  -  j(w)|  <  Jr  |-  n(tN  +  oN)|(|v|T  -  |wT|)ds 

c 

3 

"1 

<  »-n(tN  +  oN)i0>rcJvT  -  wTi0trc 

■» 

<  ll-ti(tN  +  erN)l0jpclv-w|  j* 

(4.3.22) 

$ 

>.*  ■  ’ 

Thus,  the  variational  formulation  (4.3.18)  is  well-defined  in  the 

Sobolev  space 

a 

H^CQ),  since  f2  is  clearly  linear  and  continuous. 

-  : 

*  1  •  ; 

A  solution  u  to  the  general  problem  (4.3.5)  might  be  obtained  by  the 

iteration  process 

H 

(i)  Solve  the  Special  Case  I  by  assuming  the  tangential 

( 

-  - 

-  *  • 

stress  Oj.  As  a  result,  the  normal  stress  aN  is 

computed. 

• 

1 

■  • ' 

(11)  Using  the  computed  normal  stress  o^,  solve  the  Special 

> (4.3.23) 

r_  - 

- 

(4.3.18)  Case  II.  As  a  result  the  tangential  stress 

r 

Of  is  computed. 

F» 

( 

T  ' 

(ill)  Check  the  convergence  of  the  solution.  If  convergence 

is  not  attained,  repeat  (i)  and  (ii). 

> 

We  know  of  no  convergence  proofs  for  this  iteration  process, 

although  the 

• 

i 

• 

( 

t~  - 

► 

►  . 

>  '  . 

1 

.  1 
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process  has  proved  to  converge  to  reasonable  approximate  solutions  to  many 
problems . 


4.  Special  Case  I:  Prescribed  Tangential  Stress 
As  mentioned  earlier,  we  shall  solve  the  Signorini  problem  (4.3.5)  with 
friction  by  using  the  iteration  process  (4.3.23)  that  is  the  sequential 
approximation  of  the  problem  (4.3.5)  by  the  two  special  cases  (4.3.13)  and 
(4.3.18).  In  this  section,  the  first  case  (4.3.13),  which  corresponds  to 
situations  in  which  the  tangential  stress  oT  (possibly  due  to  friction)  is  known 
as  a  function  of  L^(TC),  i.e., 

2t  “  S,T»  St  €  L2<rc),  (4.4.1) 

where  L2(rc)  =  (t  -  tpip  :  tp  6  b2(rc)}.  As  shown  in  (4.3.13),  a  variational 
formulation  to  this  special  case  is  governed  by 

u  €  K  s  B(u,w-u)  >  f i<w-u) ,  V  w  6  K,  (4.4.2) 

where 

K  -  (v  6  V  :  vN  -  g  <  0,  a.e.  on  rc) ,  (4.4.3) 

V  =»  (v  £  H*(Q)  :  V£  -  uif  a.e.  on  r^,  1  <  i  <  3) ,  (4.4.4) 

and  H*(Q)  is  the  Sobolev  space  defined  by  (4.3.11).  Since  the  bilinear  form 
B(»,»)  and  the  linear  form  f(»)  are  well-defined  on  H*(Q),  we  can  define  the 
variational  formulation  (4.3.13)  by  using  the  Sobolev  space  H*(Q)  as  (4.4.2). 
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The  subset  K  of  H^(Q)  is  the  so-called  admissible  set  of  all  possible  functions 
satisfying  the  given  boundary  condition  on  and  the  kinematic  contact 

condition  on  and  producing  finite  energy.  The  set  V  consists  of  all 

functions  that  satisfy  only  the  boundary  condition  on  r^. 

4.4.1.  Existence  of  a  Unique  Solution 

Because  of  the  assumption  on  the  elasticity  constant  (4.2.4),  Korn's 
inequality  yields  the  property 

B(v,v)  >  mlvlj,  v  £  V  (4.4.5) 

where  m  is  a  constant  independent  of  v  V.  Details  of  Korn's  inequality  can  be 
found  in  e.g.,  Hlavacek  and  Necas  [1970].  We  also  recall  that  (see  (4.3.6)  and 
(4.3.16)) 

B(v,w)  <  MllvB^lwllp  v,  w  £  H*(Q)  (4.4.6) 


and 


If  l(v)  I  <  »f  1*  v  €  H1(Q), 


(4.4.7) 


where  Ifjlj  -  lfl0  +  Itl^  +  l2«0,rc  +  "St'o.IV 


These  properties  are  sufficient  to  guarantee  the  existence  of  a  unique  solution 
u  (  K  to  the  problem  (4.4.2)  since  the  set  K  is  closed  and  convex  in  H*(Q)  under 
the  assumption  of  the  Lipschitz  domain  Q. 

Theorem  4.1.  Let  the  domain  Q  is  Lipschitzian,  and  let  (4.2. 3)— (4.2. 5) 
hold.  Suppose  that  (4.4.1)  holds.  Then  there  exists  a  unique  solution  u  £  K  to 
the  variational  inequality  (4.4.2): 
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u  £  K  :  B(u,w-u)  >  fj(w-u),  w  €  K 

Proof .  See  e.g.  Lions  and  Stampacchia  {1967].  0 

4.4.2.  Penalty  Resolution  of  the  Contact  Condition 

In  the  formulation  of  (4.4.2),  there  is  an  inequality  constraint  Vjj  -  g  <  0 
a.e.  on  Tq,  that  leads  to  the  inequality  form  (4.4.2)  of  the  variational 
formulation.  We  shall  resolve  the  condition  vN  -  g  <  0  by  using  exterior 
penalty  methods.  Introducing  penalty  parameter  e  >  0  such  that  e  ♦  0,  the 
variational  inequality  (4.4.2)  is  approximated  by 

ue  €  V  :  B(u£,v)  +1  (P(u£n)«vn)  =  X  €  VQ,  (4.4.8) 

where 

P(<P)  =  (4»  -  g)+,  <t>+(x)  =  sup{<)»(x)  ,0} ,  a.e.  on  Tc,  (4.4.9) 

(•,•)  is  the  inner  product  of  L^(T^)  such  that 

( <t>  ,4»)  3  Jr  4>4><ir ,  4>  €  L2(rc)  (4.4.10) 

and  Vq  is  a  subspace  of  H*(Q)  defined  by 

V0  -  {v  €  HX(Q)  :  vt  -  0  on  1  <  i  <  3} .  (4.4.11) 

We  expect  that  u£  converges  to  u  as  e  0. 
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The  physical  idea  behind  the  penalty  approximation  (4.4.8)  is  that  the 
rigid  support  can  be  approximated  by  the  Winkler  foundation  consisting  of 
continuously  distributed  springs. 

Theorem  4.2.  Under  the  same  conditions  of  Theorem  4.1,  the  sequence  u£  of 
the  solutions  to  (4.4.8)  converges  to  the  solution  of  the  variational  inequality 
(4.4.2),  as  e  +  0. 

Proof .  We  first  note  monotonicity  of  the  operator  0.  Indeed, 

(P(a)  -  0(b) ,  a-b)  =  (0(a)  -  0(b),  (a-g)  -  (b-g)) 

-  Jr  {(a-g)+  -  (b-g)+H(a-g)  -  (b-g)}ds 

c 

>  Jr  {(a-g)+  -  (b-g)+H(a-g)+  -  (b-g)+}ds, 


i  .e. 


9 


(0(a)  -  0(b),  a-b)  >  l(a-g)+  -  (b-g)+l£  r  >  0  (4.4.12) 

for  every  a,  b  6  L^(TC).  Here  it  has  been  applied  that  <|>  =  $+  -  4»~, 

<>~  ■  sup{-$,0},  4>+$-  =  0,  and  >  0,  for  every  cp  C  R» 

Since  0(w^)  =  0  for  any  w  £  K,  (4.4.8)  implies 


ie>  +  7  (P<ueN> *wN_ueN)  “  fl^"lie> 


and 
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"7  (P<WN>  “  *wN_ueN)  =  (4.4.13) 

Applying  the  monotonicity  of  the  operator  0  yields 
B<“e»£~!P  *  fl(5“He)*  w  6  K. 

Applying  <4.4.5)— (4.4.7) ,  we  have 

raHUglj  <  M»ug»  Llwl  t  +  lug.lj) 

Using  Young's  inequality 

<J>4,  <  6<}>2  +-^  <J>2»  6  >  0,  (4.4.14) 

it  can  be  concluded  that 

Bu^l ^  <  C(m,M, I f i* 1 » Iwl p  <  +*. 

In  other  words,  the  sequence  {u£}  of  the  solutions  of  (4.4.8)  is  uniformly 
bounded  in  H*(Q)  in  e.  Thus,  sequential  compactness  of  the  Sobolev  space  H*(Q) 
implies  the  existence  of  at  least  one  subsequence  of  {u£},  still  denoted  by 
{ue},  that  converges  weakly  to  a  limit,  say  u  6  H^Q). 

Now  we  recall  that  convexity  and  Gateaux  differentiability  of  the 
fifnctional  v  -*•  B(v,v)  in  H^Q)  yields  weak  lower  semicontinuity  of  the 
functional.  Thus,  to  such  a  weak  weak  convergent  subsequence  {u£},  we  have 

B(u,w-u)  >  lim  sup  B(ue,w-ue)  >  fj(w-u), 


B(u,w-u)  >  fj(w-u),  w  6  K. 


(A. A. 15) 


On  the  other  hand,  relation  (A. A. 12)  and  equation  (A. A. 13)  imply 


<ueN-g>+,o,rr  <  ~  f<rHe>>- 


(A. A. 16) 


for  any  w  6  K.  Continuity  of  the  norm  and  uniform  boundedness  of  lu£lj  yield 


0 (uN~8)+lo,r  =  °*  i,e*‘  un  “  8  <  0  a*e*  on  r( 


Since  u£  V  and  since  V  is  closed  in  H1(Q),  the  limit  u  of  {u£}  belongs  to  K. 
Thus,  the  limit  of  a  subsequence  {u£}  is  a  solution  of  (A. A. 2). 

Because  of  the  uniqueness  of  the  solution  to  (A. A. 2),  every  convergent 
sequence  must  have  a  common  limit.  Therefore  the  sequence  {u£}  itself  as  well 
as  convergent  subsequences  has  to  converge  to  the  solution  u  6  K  of  (A. A. 2).  Q 
Let  the  approximation  of  the  contact  pressure  (normal  stress)  o£N  be 
defined  by 


aeN  "  -  7  P<ueN>  in 


(A. A. 17) 


It  can  be  shown  that  (A. A. 17)  is  equivalent  to  the  variational  inequality  form 


°eN  €  M0  :  (xN  "  °eN*  ea£N  +  ueN  "  8>  >  0,  tN  €  Mo,  (A. A. 18) 
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M0  =  {tN  6  L2(rc)  :  1N  <  0  a.e.  on  Tc} .  (4.4.19) 

w  (4.4.18),  it  suffices  to  note  that  for  x,y,  and  f  (  R,  the  solution 
o  the  problem 

x  <  0  :  (ex  +  f)(y  -  x)  >  0,  y  <  0 

n  by  x  =  —  f  . 

us  the  exterior  penalty  formulation  (4.4.8)  of  the  variational  inequality 
is  equivalent  to  the  following  perturbed  Lagrangian  multiplier 
tion: 

<!ie»0eN>  ^  Y  x  Mf)  : 


B(u£»v)  ~  ^°eN,vN^  =  v  £  V0  (4.4.20) 

(xN  "  °eN*  eoeN  +  ueN  "  8)  >  °»  XN  ^  **0 

here  are  no  guarantees  that  contact  pressures  belong  to  L2(TC),  the 
ance  argument  of  (4.4.20)  in  H^(Q)  *  L2(T^)  does  not  have  any  sense  here. 
Lf  a  flat  rigid  punch  is  indented  into  an  elastic  foundation,  the  contact 
a  Ojy  is  generally  not  in  L  (T^).  Thus,  we  need  to  consider  convergence 
3eN)  in  larger  spaces. 

the  fractional  Sobolev  space  HB(D,  0  <  S  <  1,  is  defined  by  the 
»  of  L2(D  equipped  with  the  norm 
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6$;s,r  =  1 11  ^ a o , r  +  h  h  n+2sm~  drdr  ^ 

i^-^i 


(4.4.21) 


where  QC  RN  is  Lipschizian,  the  extension  of  a  function  v  £  H*(Q)  to  the  boundary 
1 /2 

belongs  to  H  (T)  by  the  trace  theorem,  see  Necas  [1967,  Chapter  2].  suppose 
that  the  surface  of  the  rigid  foundation  is  polygonal.  That  is,  assume  that 


1/2 

N  €  H  -6(rc),  6  >  0 


(4.4.22) 


for  a  small  positive  number  6.  Then  the  normal  displacement  vN  of  v  6  H^(Q)  is 


in  the  space 


1  /2_* 

Q  =  H  (rc),  6  >  0. 


(4.4.23) 


Let  Q'  be  the  topological  dual  of  Q,  and  let  [•,*]  be  the  duality  pairing  on 
Q'  x  Q.  The  partial  ordering  in  Q'  is  defined  by 


x  £  Q'  ,  x  <  0  if  and  only  if  xn  €  Mq  s.t.  xn  -*■  x  in  Q'  (4.4.24) 


Then  the  convex  subset 


M  =  {x  £  Q'  :  x  <  0} 


(4.4.25) 


is  closed  in  Q' .  Note  that  since  Q  is  densely  imbedded  in  L^(Tq),  the  duality 
pairing  [*,*]  has  the  property  that 


[ P ,4> ]  =  (t,4>)  =  /r  x$ds,  x  £  L2(TC),  6  Q. 

c 
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thus  the  L^(r^)-inner  produce  (•,•)  in  (4.4.20)  can  be  replaced  by  (•,•]  as  long 
as  g  6  Q. 

Theorem  4.3.  In  addition  to  the  conditions  in  Theorem  4.1,  suppose  that 
(4.4.22)  holds.  Then  the  sequence  t(ue*°eN^  €  V  x  Mq  of  the  solutions  to 
(4.4.20)  converges  to  the  solution  (u,a^)  (  V  x  M  of  the  problem 


B(u,v)  -  [oN,vNJ  =  fL(v),  v  £  V0 


hN  "  °N»  uN  -  g]  >  0,  xN  €  M 


> 

/ 


(4.4.26) 


as  e  +  0. 

Proof.  It  suffices  to  show  that  fl q'  *s  uniformly  bounded  in  Q'  and  in 

1/2 

e.  To  do  this,  we  note  that  since  the  trace  map  y  :  H^Q)  +  H  (T)  is  linear 

1/2 

continuous  and  surjective,  and  since  the  map  £  *  ~  <J>*N  from  H  (Tq)  iuto 

H  /  _6(rc)  is  linear  continuous  and  surjective,  it  can  be  found  an  element 
v  £  Vq  for  every  <J>  £  H1 '  -^(rc)  such  that 

v^  =  4>  and  Bvllj  <  C>  <J>>  q,  (4.4.27) 

for  some  constant  C  >  0  independent  of  the  choice 
Then 


,aeN,Q' 


sup 

D  €  Q 


t°eN*^ 


I*!, 


<  C  sup 

-Q 


to 


eN»vNJ 

Ivl  j 


where  Vq  is  the  set  of  v  satisfying  (4.4.27).  Applying  (4.4.20)p  we  have 
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*oeN#Q'  <  C{M«ue«  j  +  llfjl*}. 

Since  the  uniform  boundedness  of  II u£ I  ^  is  obtained  similarly  to  Theorem  4.2,  the 
uniform  boundedness  of  is  assured.  Q 

4.4.3.  Finite  Element  Approximations 

There  are  at  least  three  kinds  of  finite  element  formulations  of  the 
variational  problem  (4.4.2).  If  aproximation  is  performed  through  the 
variational  inequality  (4.4.2)  itself,  a  finite  element  formulation  of  the 
so-called  displacement  type  can  be  obtained.  This  method,  however,  does  not 
involve  the  contact  pressure  and  requires  an  auxiliary  scheme  to  compute  it 
with  sufficient  accuracy.  A  second  approach  is  to  obtain  a  mixed  finite  element 
approximation  through  the  formulation  (4.4.26)  using  the  Lagrangian  multiplier 
method.  In  this  case,  the  contact  pressure  itself  becomes  an  additional 
unknown  to  be  calculated  along  with  the  nodal  displacements.  Naturally, 
information  on  the  deformation  and  the  contact  pressure  can  be  easily  obtained 
by  this  method.  This  approach,  of  course,  involves  the  approximation  of  two 
independent  fields,  (u  and  oN)  and  may  be  computationally  inefficient.  A  third 
approach  is  to  use  penalty  methods  to  handle  the  contact  constraint.  Then  again 
only  one  unknown  field  appears  in  the  variational  problem;  the  penalized 
displacement  field  u£.  An  approximate  contact  pressure  can  be  obtained  later 
using  (4.4.17)  after  solving  the  displacement.  Thus,  in  this  case,  the  contact 
pressure  is  computed  by  a  simple  post-processing  procedure.  In  the  present 
study,  the  penalty  method  will  be  applied  to  solve  the  problem  (4.4.2). 
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The  following  finite  element  analysis  will  be  performed  for  the  case  that 
the  domain  Q  is  a  bounded  polygon  in  R  .  That  is,  all  deformations  of  the  three 
dimensional  body  are  projected  in  the  plane,  and  the  domain  Q  can  be  covered  by 
finite  elements  that  have  straight  edges.  Extension  to  three-dimensional  and 
curved  finite  elements  is  straightforward,  although  technically  very  complex. 

Let  the  domain  Q  be  covered  by  E-finite  elements,  i.e., 

E 

5  -  U  0e' 

e=l 

As  usual,  we  suppose  that  each  element  Qg  contains  Ne  nodal  points  on  which 

N 

shape  functions  are  assigned  so  that 

Ne 

4»®(xj)  =  6^,  l  4<i(x)  =  1  x  €  Qe 
i=l 

Ne 

where  is  the  set  of  nodal  points  in  Q£.  The  functions  u£  and  v  in 

(4.4.8)  are  approximated  by  piecewise  polynomials  u£  and  vh  defined  by 

uei 1  Qe  =  vj'Qe  =  vj<^~)  (4.4.28) 

in  each  element  Qe,  respectively.  Here  u^Iq  means  the  restriction  of  the  i 
component  of  u^  into  the  subdomain  Qg.  Approximation  (4.4.28)  leads  to  the 
following  finite  element  approximation  of  V  and  Vq: 

Vh  -  {vh  :  vh  =  v}^,  vj  €  C°(Q),  vJ|Q  -  vj(^(x), 

e 

vi  *  ^(xj^  on  rDi»  1  <  i  <  2},  (4.4.29) 

L  L  L 

and  Vq  is  a  subspace  of  Vn  containin  functions  such  that  ^(x^)  -  0. 


D 

It  is  worthwhile  to  note  that  the  term  3(U£n)  may  not  be  a  piecewise 
polynomial  even  though  is,  because  of  the  special  operation  41  +  4>+.  Then  the 
penalty  term  cannot  be  integrated  as  the  other  terms  when  the  stiffness  matrix 
and  the  generalized  load  vector  are  computed. 

To  resolve  this  difficulty,  let  us  introduce  a  quadrature  rule  to  evaluate 
the  penalty  term  (recall  Chapter  3  of  Volume  II): 


(U<^n),v>5)  i  I  (fS<UgN)vg) 


h  \„h' 


E' 


e=l 


> 


Xe ( ^ ^ucN^vN )  “  l  wjP(ueN^iJ^vN^J^' 
J=1 


J 


(4.4.30) 


Here  {wj}  and  {xj}  are  the  sets  of  weight  and  integration  points,  G  the  number 
of  integration  points,  and  E'  is  the  number  of  finite  element  edges  covering  the 
contact  boundary 

The  penalty  formulation  (4.4.8)  is  thus  approximated  by  the  following 
reduced-integration  penalty  finite  element  method: 


uhe  6  _Vh  :  B(u£,vh)  +1  I(P(u£n)v5J)  =  f(vh), 

vh  i  (4.4.31) 

Because  of  the  term  problem  (4.4.31)  is  nonlinear,  and  the 

stiffness  matrix  due  to  the  penalty  term  is  not  readily  computed.  The  actual 


computation  depends  upon  the  method  we  use  to  deal  with  the  nonlinear  term 
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To  illustrate  one  alternative,  suppose  that  1ujl  is  the  approximation  of  the 


i  step,  and  let 


PiC1^*)  = 


lueN  "  8  if  1_lueN  "  8  >  0 


lf  i_lueN  *  8  <  0. 


(4.4.32) 


wherein  i  is  not  summed.  Then  the  nonlinear  equation  (4.4.31)  can  be  solved  by 
the  successive  iteration  scheme 


iuhE  6  Vh  :  B(iu^ ,yh)  +  1  l(pi(i»eN)v}))  =  f(vh),  vh  6  vg  (4.4.33) 


h  r  ,,h 


for  i=l,2,3,...,  and  the  initial  assumption  ug.  Since  (4. '.33)  is  linear  for 
each  i,  the  stiffness  matrix  and  load  vector  can  be  obtained  by  stancard  finite 
element  procedures  for  linear  problems  discussed  in  great  detail  in  Vols.  I  and 
III.  Indeed,  for  each  finite  element  we  have 


B|a  ‘  viKi«ipuj 


Ua  -  c8i«) 


(4.4.34) 


flo  (vh)  =  v?ff. 


where 


^ia  jp  “  -^Q  E i k  j  ifta , k^p  , 1 


1 


54 


The  assembly  and  solution  of  these  linear  equations  is  now  done  following  any  of 
the  standard  procedures  described  in  Vol.  III. 

Since  the  coercivity  of  the  bilinear  form  B(»,»)  and  continuity  of  B(»,») 
and  the  linear  form  f(«)  are  independent  of  the  choice  of  the  penalty  parameter 
e  >  0  and  the  ste^  i  of  the  iteration  (4.4.33),  arguments  similar  to  those  in 
Theorem  4.2  yield  the  following  result. 
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Theorem  4.4.  Under  the  same  conditions  in  Theorem  4.1.,  if  the  weights 
{wj}  of  the  quadrature  rule  (4.4.10)  are  all  positive,  then  the  sequence  {1u^} 
of  the  solutions  to  (4.4.33)  converges  to  the  solution  u^  of  (4.4.31)  as  i  + 
<®.  In  addition,  the  sequence  {u^}}  of  limits  obtained  for  each  e  >  0  converges 
as  e  tends  to  zero  to  the  solution  uh  £  K*1  of  the  variational  inequality 

uh  €  Kh  :  B(uh,vh-uh)  >  f(vh-uh),  vh  £  Kh  (4.4.38) 


where 


Kh  =  {vh  £  vh  :  (vj  -  g)(xj)  <  0  on  Tc} .  (4.4.39) 

□ 

It  is  noted  that  the  iteration  process  finding  uj?  for  each  e  >  0 

converges  rather  quickly  for  example  problems  solved  later,  although  the  results 
in  Theorem  4.4  imply  only  the  convergence  of  the  sequence  {^u^}  as  i  >  ®. 

4.4.4  Convergence  of  the  Penalty/Finite  Element  Method 

We  shall  define  an  approximation  of  the  contact  pressure  o^N  as  in  (4.4.17): 

f  Qh  and  °£N<5j)  '  -  £  6  (u£N>  'Jj'- 
J  =  1 . G,  on  rj,  e  =  1,  2 . E',  (4.4.40) 

v> 


where 


^h 


TgMg,  1  <  e  <  E'  ) 


(4.4.41) 


and  {M„}  is  the  set  of  shape  functions 
P 

the  quadrature  rule.  Each  function 

h 

boundary  T  .  Continuity  of  T  depends 
Applying  (4.4.40)  to  the  equation 


associated  to  the  integration  points  of 
in  might  not  be  continuous  along  the 
upon  the  choice  of  the  quadrature  rule. 
(4.4.31)  yields 
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j  h  h  f  .,h  .  .h 

<V  V  *  l  x  M  ; 


_  »  h  h*  t  /  ^  bv 

B<V  V  )  -  I(<,eN  v„> 


f  <vh)  ,  V  vh?  vh 


(4.4.42) 


I  (  (t  h  .  ,  h  ,  h  O'  ^  n  u  h  £  h 

^  N  -  °eN  (e0£N  +  ueN  "  *>)  >  °*  V  TN  C  M  ’ 


where  M°  is  a  subset  of  Q*  defined  by 

ft 


Mh  =  {xh  €  Q;  :  Th  (Xj)  <0,  1  <  J  <  G  on  r®  ;  1  <  e  <  E*  }- 


(4.4.43) 


Lemma  4.1  Let  (u,  o„)  V  x  M  and  (u*\  a*\,)  f  V*1  x  M*1 

-  ^  N  ^  £  N  ^ 

be  the  solution  of  (4.4.26)  and  (4.4.42),  respectively.  Let  E^.  be  the 
quadrature  error  of  the  numerical  integration  introduced  in  (4.4.30): 


Ej  (f,g)  =  (f,g)  -  I(fg) 


(4.4.44) 


for  properly  defined  functions  f  and  g  so  that  the  quadrature  rule 

h  h  b  h 

(4.4.30)  makes  sense.  Then  for  every  v  €  K  and  Tjj  €  M 


_  /  h  hv  n  /  b  hv  ,  r  h 

B(u  -  u£,  u  -  u£)  <_  B(u  -  u£,  u  -  v  )  +  [oN>  vN 


+  ^°N  *  tN’  UN  “  UeN^  +  [tN  "  °N’  UN  * 


-  h  (tn  <ucN  -  8))  +  1  (  (tn  ■  tcB> 


h  .  h\ 

_v.)  t  e, 


(4.4.45) 
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1 


under  the  assumption  that  the  "gap"  function  g  is  well  defined  so  that  the 
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Then  the  estimate  (A. A. 45)  follows  from 


h  hN  .  h  - 

3en  -  V  <U£N  -  s> 


(<0^N  - 


/  h  v 

V  ‘-“sV 


£  -  I  ((»£  -  tn>  * 


S)E  •  D 


It  is  noted  that  the  estimate  (4.4,45)  is  obtained  within  the 
restricted  sets  K*1  and  M*1.  The  first  four  terms  of  the  right  hand  side 
of  (4.4.45)  are  related  to  the  interpolation  error  of  function  by  finite 
element  methods,  the  fifth  term  comes  from  the  ir tegration  error  for 
the  penalty  term,  and  the  last  one  is  from  the  method  of  penalty.  It 
is  expected  that  if  £  -+■  0,  the  last  term  goes  to  zero.  We  shall  study 
this  more  precisely. 

Toward  this  end,  let  us  introduce  an  approximation  of  the 

normal  trace  operator  v  -*■  v  on  H^(Q): 

n 


'V*h  (*h)1  ■  (v  Ah  (’h> )  ■ 


N’  h 


T  ,  h  h. 

1  <T  V 


(4.4.46) 


c  h  s  h 

for  every  v  f-  V  . 


Lemma  4.2  Let  (u,  0  )  £  V  x  M  and  (u^,  )  £  Vh  x  Mh  be  the  solutions 

-  ~  nv'~  ~  ~ 

of  (4.4.26)  and  (4.4.42)  respectively.  Suppose  that  there  exist  a  positive 
number  independent  of  tjj  and  an  element  v*1  £  such  that 
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\<?h>-^  “k  H-h  »i  i  Iloilo.  rf 


(4.4.47) 


for  a  given  €  Mh  R  cy,  where  (A^)  is  the  range  of  the  operator 
A^.  Then,  for  every  M*\ 


(<$  -  »CN>2)  I  <  "II  «  -  &H  ill  - 


°cn'I  o,r  /a  h 

V 


+  II TN  ~  aN II  0,  r  II  TN  °eN  II  0,  r_  +  iEI  l  tn  (t  °eN)) 


(4.4.48) 


Proof .  From  (4.4.26)  and  (4.4.42), 


h  h  „  h  N 


1  \(tN  “  aeN)vN 


1  ((aN  0gN)Vn)+  1  ((XN  "  °n)vn) 


+  [V  VN]  '  1 (tN  +  I(aN  VN} 


"  loH.  vj]  +  i((tJ  -  aN)  v^ 


n f  Vi  Vi  *  t  i  .  _  /  Vi  Viv 

B(u  -  u£f  v  )  -  [0N,  vN]  +  I  (tn  vn) 


_  f  h  h»  r  ti  h «  p  /  h  tiv 

B(u  -  u£,  V  )  +  [tn  -  o„,  VN]  -  Ej  (T„.  v„> 


Applying  (4.4.47)  and  taking  absolute  value  yield  the  result.  0 
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Let  us  now  apply  the  interpolation  results  of  conforming  finite 
element  approximations  on  the  restricted  sets  K^1  and  M*1  following  Falk 


[1974] 


there  is  a 


h  r 

?  i 


K  such  that 


vh  -  v||  r  <  Cx  hpl||  v  Ils 


for  every  v  ..  K,  (Q) ,  r  _<  s,  and  s  _>  2 ,  where 


=  Min  {k  +  1  -  r,  s  -  r}. 


and  k  is  the  order  of  the  complete  piecewise  polynomials 
contained  in  finite  element  approximations  of  v. 


(4.4.49) 


there  is  a  tJJ  M^1  such  that 

II  tn  '  tn II  P,r  -  c2h 

V 


N11  q,rr 


for  every  xN  ~  MOH^Cr^.),  p  _<  q,  and  q  >  ^  +  5, 

6  >  0,  where 

P2  =  Min  (t  +  -  -  p,  q  -  pi 

and  t  is  the  order  of  complete  picewise  polynomials  contained 
finite  element  approximations  of  T^. 


> 


(4.4.50) 


We  shall  suppose  that  the  parameters  r  and  p  in  (4.4.49)  and  (4.4.50), 
respectively,  can  be  negative  numbers  (see  Babuska  and  Aziz  [1972,  p.95j) 
and  that  M*1  need  not  be  conforming. 


Theorem  4.5.  Suppose  that  the  quadrature  rule  I  satisfies  for  some 


ive  constants  C^,  i  -  3,  ....  6,  and  parameters  and  A  the  in- 
ities 


KvJ, 


)  >  c3| 


I2 

'o,r 


T  /  h  Wv 

T^TN  vN)  -  c4 


h  ,| 

tmI! 


o,rr 


n  o,r 


(4.4. 31) 


Ei(  V  V  I  -  Sh  II  tn  II  q , r  I!  v5 II  s-l/2,rr  ’ 


In  _h  ii  m  h i 

C 


EI  (V  tN)  I  -  c6  h  II  TN II  q,  r  II  fN  1 1 0  ,  T 

c  c 


EI  (V  8>  I  -  C7h  II  TN  I!  q,r  II  gll  s-l/2,r 


(4.4.52) 


q  and  s  are  the  parameters  in  (4.4.49)  and  (4.4.50). 
i  under  the  regularity  assumption 


<«)»  s  >  2,  0N  £  HS_3/2(r  ) ,  and  g  £  HS_1/2  (T  ) 


(4.4. 53) 
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the  following  error  estimates  are  obtained: 


!  '  ^  1  -  C  (MI  3*»  °N«  s-3/2, r  •ll*Ha.l/2,  >  <h  3  +  E/Ofc  ■*  ^  S 


°N  °£H  I'  0,r  -  c(!l  If  s-3/2 ,  r  ^  (h  4  +  II  u  -  u^ll  1/c l  ) 

c  *  C  "  1 


(4.4.54) 


where 


y3  =  min  {k,  s-1,  (k  +  s-l)/2,  t,  (t  +  s) /2 ,  X^l ,  X^/2) 


y^  =  min  {t  +  1/2,  s-3/2,  X^ }  . 


(4.4.55) 


Proof.  Applying  (4.4.51)  and  (4.4.52)  into  (4.4.48)  yields 


S»  t 


"w11  o,rci  M  l!“  -  "ell  i/ah 


h  I 

tn  -  V 


+  c6  h  2 11  ’S"  ,.rf 


Applying  the  interpolation  theorem  (4.4.50),  we  have 


h  h  i- 

TN  "  aenllo,r_  -  L 


(  -  “eHl/ah  +  U  +  l/2,q,X2}  ||  TJ(|  r>  (4.4. 56) 
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On  the  other  hand,  from  the  estimate  (4.4.45),  we  obtain 


m 


I!  u  -  u*||  l  <  M||  u  -  »J||  Jl  «  -  vh||  t  +  ||  aN||  ^^11  uN  -  vj||  ( Hq(rc))  1 

+  !!°N  ‘  TNH  ^Hl/2(r^}^  J1  UN  "  UeN!i  1/2,  Tc  + 


°N  “  TNI!  (Hs-l/2(r^))  ,'1  UN 


8'L-l/2  +  C5h  l!'TN^q,r  ^UcN^s-l/2,  P 


+  C7h  H  TN  ^  q,T  H  8 II  s-l/2,I\ 


+  c4  ell  XN  -  °eNll  0,r  HtnII  o,r 

L*  C 


Substitution  of  (4.4.56)  and  applying  Young's  inequality  imply 


U  -uj||2i  <  c  ||u  -  vh||2i||  + 


ctnII  q,rc II  un 


-  "ill/  X 


+  "°H  -  (s1/2(rc>) .  +  "  "■  -  8  «-x/».rcP  ^s.1/2(  J  f 


+  Sh  'll  th'I  q,rcll  u"»ll  s-i/2, rc  +  c7h  'll  T»H  q,rcll  *  H5-i/2,r 


3  it  h, 


-  1 
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2 {Min  t  +  1/w,  g,A  }  ,  „ 

i-  I!  J*  I!  *■ 


*  ca  (E/oth»2 II  tnII  o,r  +  2  WU 


*  T 


n  o ,  r 


putting  q  =  s  -  3/2,  we  arrive  at  the  bound 


II  ,  5  c(ll  Ml  s.n  oN||  s_3/2iI.c,  IU  II  ,  .  1/2, gfh"1"  f-5'1} 

in  jk+s-l,2s-2^/2  +  hmin  |t,s-lj  +  hmin  |^+s , 2s-2j /2 
\/2  V2  ,  .  min  (t+1/2 ,  s-3/2, 

+  h  +  h  +  (e/ah)  +  eh  J  ( 


l“"Helll-C0^^s,^ON  IU-3/2,r  ’  5-1/2,  r  )(h 

v  c  c# 


+  e/c^  +  eh 


where 


=  min 


|k,s-l, 


(k+s-1) / 2 ,  t,  <t+s)/2 


,  A1/2,  X3/2j 


=  min  ^t+1/2,  .  3/2,  A2^  • 
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Now  an  estimate  for  the  approximate  contact  pressure  is  derived 
immediately  from  (4.4.56): 


ii  II  u  i« 

II aN  ~ °eN  II o, r  -  HaN  “  tn  Ho,r  +  Htn  “  aeN  ||o,r 

c  c  11  c 

i  C  (llaN  lls-3/2, r  )(h  4  +  11“  "  ~e  lll/ah)  •  D 


Tt  is  clear  from  these  results  that  three  kinds  of  estimates  of  the 
error  introduced  by  the  numerical  integration  I  are  needed  as  well  as  an 
estimate  of  the  parameter  of  (4.4.47)  for  each  given  choice  of  integra¬ 
tion  rule  and  finite  elements  basis.  We  shall  discuss  these  in  the 
subsequent  subsections. 


4-Node  Isoparametric  Element  and  Trapezoid  Rule 
One  of  the  simplest  finite  elements  to  the  contact  problem  is  the  4-node 
isoparametric  quadrilateral  element.  This  element  yields  piecewise  linear 
approximation  of  the  displacement  on  the  boundary.  Moreover,  the  edges 
are  just  straight  lines.  For  the  integration  rule  "I"  ,  let  us  apply 
the  trapezoid  formula  on  each  boundary  element  that  correpsonds  to  an  edge 
of  4-node  element.  The  "normal"  component  v^  of  v*1  on  the  boundary 
element  is  constructed  by  taking  the  same  interpolation  to  v^  by  using 
the  value  v*1  •  N  at  each  nodal  points  on  the  boundary.  That  is,  vj!  is 
a  piecewise  linear  polynominal  on  which  is  the  same  degree  as  v\ 

On  the  other  hand,  the  trapezoid  rule  on  each  boundary  element  yields  the 
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same  approximation  of  the  contact  pressure  T  of  the  "normal"  displacement 

.  That  is,  both  functions  vj^  and  are  linear  polynomials  on  each 


N 


boundary  element  T 


Lemma  4.3.  For  the  functions  v!l  and  defined  in  above,  we  have 
-  N  N 


‘(VN  )±llVNllO 

*c 

[(TtS  vj)  ^  Iloilo, ,l|vhNll0>, 


(4.4.57) 


E' 


Proof.  Noting  that  ||vjj  ||*  =  I|vJ  \\2Q  , 

’  r  J  o 


e=l 


Te  =  3fljlrc  ,  we  shall  derive  the  inequalities  of  (4.4.57)  on  each 
boundary  element  1^  .  Since  vjj  and  are  linear  polynomials  on  T  , 

we  can  express  these  as 


h  h 

v__  =  a  +  bs  and  T„  =  c  +  ds  , 
N  N 


(4.4.58) 


where  s  is  the  local  coordinate  along  the  boundary  element,  and  {a,b,c,d} 
is  the  set  of  numbers  defined  by  the  values  of  vjl  and  ill  at  nodes  and 
the  length  of  the  element  fe  .  Then  we  have 


ibSC 


(a  +  bs)  ds 


-h(, 


2 

h(  a2  +  abh  +  h^ 


) 


L 
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‘(VN2)  '  h(‘ 


a2  +  nbh  +  y  h2 


)t  HvnII  o,r 


(tnvS)  i 


h2 

I(TNh  ) 


h2 

KvJ  ) 


-  3  Htn  ^o, r  ^vn  ^o.r 
c  c 


We  now  show  error  of  the  quadrature  rule  I  for  the  choice  of  the  4-node 

element  and  the  trapezoid  rule  to  the  penalty  terra* 

t  /.  /.  T? _ i _ r _ • _  h  h  h  .  .  _ _ 


Lemma  4.4. 

For 

the  functions  x 

estimates: 

1  Ex  (tn 

•  vn)| 

iCih2  ||tJ||1/2 

|Ei(TN. 

tn)| 

<  C2h  HtnIIi/2,: 

IEi(TN 

.  e)  l 

-  C3h  ^  TN  ^1/2, 

rc  TN'ii/2,rc 


(4.4.59) 


where  C^,  ^2*  and  are  positive  constant  independent  of  the  mesh  size  h. 


Proof .  Noting  that T ^  and  are  piecewise  linear  on  the  boundary 


element. 


lEi(vvS)l  1  *  (tS)(vn)'  1,31  i  £  h2  ll(TN>ll<H1/2<r  , 


n/  1*1/2, r.  -  cih  IItn  II i/2,r  l|vN  II 3/2,r  . 


Simi lar ly. 


E’ 


h  "i,r 


Applying  the  inverse  inequality 


1,rc  ”  Ch 


■1/2 


n  ni/2,r 


(4. A. 60) 


we  obtain 


< 


i/2, r 

c 


i/2, r 

c 


The  last  inequality  in  (4.4.59)  is  obtaint  oy  the  following  triangular 
inequality 


iEl(TN  >g)l  1  lEl(TN  ’  gh)l  +  iEl(TN  * 

by  using  the  piecewise  linear  interpolation  of  g.  The  first  part  in  the 
right  hand  side  is  estimated  by 

lEl(TN  ’  g  )l  1  C3h  HTN^l/2,r  II  3/2, T  (4.4.61) 

similarly  to  the  first  inequality  of  (4.4.59).  The  second  part  is  estimated 
by  using  the  result  of  interpolation: 

lEl(TN  ’  8-8  )  I  1  c  1 1 TN  ll  1/2 ,  r  llg“8  II  (H1/2(r_))' 
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£Sh2HT«H 


n  "1/2, r 


Is! 


3/2, r 


(4.4.62) 


Thus  combining  (4.4.61)  and  (4.4.62)  yields  the  inequality  that  we  need.  Q 


The  last  preliminary  to  the  application  of  the  general  result  (4.4.54) 
to  specific  the  error  estimates  for  the  penalty  finite  element  approximation 
is  to  find  the  stability  constant  ot^  of  the  approximate  contact  pressure 
stated  by  (4.4.47). 


Lemma  4.5.  For  the  choice  of  4-node  isoparametric  elements  and  the 
trapizoid  rule,  we  have  =  o^h^2  and  vh€\/hfor  a  given  tJJ  such  that 


(4.4.63) 


where  a  is  a  positive  constant  independent  of  mesh  size  and  T 


N 


Proof.  For  the  case  xN  5  0,  the  inequality  (4.4.63)  is  obvious  for  the 
choice  v^  =  0.  Thus  we  need  to  consider  the  case  xl"!  $  0,  i.e.,  xjl €  ker(oN)  » 
where  ker(oN)  is  the  kernel  of  the  normal  trace  operator  such  that 
aNv  =  v  .  N  on  the  boundary  for  a  function  v  £  C(ft).  Applying  the  Banach 
inverse  theorem  to  the  continuous  linear  operator  o^,  we  may  condlude  that 
there  exists  a  positive  constant  C>0  such  that 


Ml  ! 


1  c  IIvn  II 1/2 , rr 


(4.4.64) 


70 


for  v*1  £  ker(o,,)Ov^  ,  where  C  is  independent  of  h. 
NT  N 


Now,  taking  the  special  case  that  at  each  integration  point. 


we  have 


■(*5)  ■ 


h2  i1/2f  h2^2 

j'*S>  r<vS> 


Applying  the  results  in  Lemma  4.3  implies 


(v»)-  11  tn  Uo.r  l|v»llo,i 


The  following  inverse  inequality  holds  under  the  assumptions  of  the 


chl/2  llvN  lli/2,r  -  llvN  ll0,r  » 
C  c 


(4.4.65) 


Thus,  we  have 


‘(vlO^^Uo.r  Hhn.> 


Finally,  applying  (4.4.64)  yields 

iW)l<*1,2Kh,rXh  ' 

V 

The  inequalities  in  (4.4.63)  now  follow  from  this  result.  Q 


Using  the  all  results  in  Lemma  4.3,  4.4,  and  4.5  and  Theorem  4.5,  we 
can  obtain  an  error  estimate  for  the  penalty/finite  element  approximation 
(4.4.31)  for  4-node  quadrilateral  isoparametric  elements  and  the  trapezoid  rule. 


Theorem  4.6.  Let  the  hypotheses  of  Lemmas  4. 4-4. 5  and  Theorem  4.5  hold 


and  suppose  that  the  solution  (u.o^)  to  the  original  problem  is  smooth  enough 
so  that 


u£H2(fi)  and  aN£H1/2(rc)  . 


(4.4.66) 


3/2 

Then,  for  g£H  (f^  )  we  ha 'e 

llu  -  uh  llj  <  c/h1  +  h_1/2E  +  h1/2S) 

l|os‘0Nllo.rc5c2(hl/2  +  h'l£  +  £) 


> 


(4.4.67) 


where  and  C2  are  positive  constants  independent  of  h.  D 


4,4.5  Numerical  Examples.  The  results  of  numerical  experiments  on  two 
example  proteins  are  given  to  demonstrate  the  performance  of  the  penalty/ 
finite  element  approximation  (4.4.31).  The  first  example  is  a  classical 
Hertzian  contact  problem  in  which  two  identical  circular  cylinders  are  pushed 
into  each  other  by  applied  line  forces  P  on  the  exterior  surfaces.  Symmetry 
of  this  problem  yields  a  Signor  ini-type  problem  since  the  deformation  is  the 
same  as  that  of  a  circular  cylinder  at  rest  on  frictionless  flat  rigid  sur¬ 
face  and  is  subjected  to  the  force  P  applied  on  the  top.  If  the  circular 
cylinder  is  long  enough,  the  problem  can  be  considered  as  a  plane  strain 
problem  on  a  circular  domain. 

As  a  further  simplification,  we  shall  solve  the  problem  only  on  the 
quadrant  applying  a  uniformly  distributed  force  instead  of  the  "point"  force 
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P  on  the  top.  Forty-four  elements  are  used  for  the  discretization  of  the 

quadrant  of  the  cross  section  of  the  cylinder  as  shown  in  Fig.  4.3.  As  an 

example,  let  the  radius  R  =  8  cm.  Young's  modulus  200  kg/ciiT,  Poisson's 

ratio  v  =  0.3,  and  the  applied  force  P  =  156  kg  on  the  top.  The  iteration 

-4 

scheme  (4.4.37)  converges  at  the  3rd  iteration  within  10  tolerance.  The  com¬ 
puted  deformed  configuration  is  shown  in  Fig.  4.4,  and  the  process  of  conver¬ 
gence  of  the  iterative  method  (4.4.37)  is  described  in  Fig.  4.5  using  the  con¬ 
tact  pressure  p^  obtained  by  (4.4.40). 

The  second  example  is  a  rigid  punch  problem  in  which  a  rigid  circular 

cylinder  is  indented  into  an  elastic  foundation,  as  shown  in  Fig.  4.6.  The 

2 

size  of  the  punch  is  R  =  8  cm,  Young's  modulus  E  =  1000  kg/cm  ,  Poisson's 
ratio  v  =  0.3,  and  the  depth  of  indentation  is  0.6  cm.  The  width  of  the 
elastic  foundation  is  8  cm,  and  its  height  is  4  cm.  We  again  assume  plane 
strain.  For  this  problem  the  relative  errors  of  the  total  strain  energy 
defined  by 


are  computed,  where  indicates  fixed  values.  For  E  ,  we  take  h  =0.8  cm 

—  /  —A  A 

and  e  =  10  ,  and  for  E^  we  take  6  =  10  and  h  =  0.8cm.  Figures  4.7  and  4.8 

contain  the  computed  result  values  of  E^_  and  E^.  Since  the  erroi  estimates 

2 

(4.4.67)  indicate  F.  =  0(e)  and  Eh  =  0(h  ),  we  have  a  slight  gap  in  numerical 


and  theoretical  results. 


Linearly  Elastic 


Pressure  Distribution  A  1-st  iteration 
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Figure  4.5  Convergence  Procedure  of 


the  Iteration  Scheme 


Deformed  Conftfur  ,it  to  n 


Figure  4.6  Indentation  Problem  for  Example  2 


Relative  Error  In  the  Strain  Energy 


4.5.  Special  Case  II.  Prescribed  Normal  Stress.  As  shown  in  the 
previous  section,  unilateral  contact  problems  can  be  solved  by  the 
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penalty  method  if  the  tangential  stress  along  the  candidate  contact 
surface  is  known.  In  this  section,  we  shall  study  the  reverse 
situation:  we  shall  assume  that  the  normal  stress  on  and,  therefore, 
the  actual  contact  surface  are  known  a  priori. 

Let 


be  given  on  the  contact  boundary  .  Then,  as  shown  in  Section  4.3, 
a  variational  formulation  of  the  corresponding  equilibrium  problem 
assumes  the  form 


u  £  V  :  a(u,v-u)  +  j(v) 


j(u)  >  f ,  V  v  €  V 


(4.5.2) 


where 


j(v)  -  /p  (  pi t jq )  |vTl  dT  with  tjj  —  tj^  +  Ojj,  (4.5.3) 

and  the  space  V  is  defined  in  (4.4.4).  Because  of  the  assumption 
(4.5.1),  the  functional  j(’)  is  well  defined.  The  fact  that  j(»)  is 
continuous  on  H*(ft)  follows  from  (4.3.22).  The  functional  j(*)  is 
also  convex.  To  see  this,  we  need  only  recall  that  the  RN-euclidean 
norm  is  convex. 


|  (l-e)vT  +  0w.p|  <  (1-0)  |Vp|  +  0|wp| 


(4.5.4) 


Thus,  a  direct  application  of  well-known  results  from  BREZIS 
[1968]  on  variational  inequalities  leads  immediately  to  the  following 
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existence  and  uniqueness  theorem. 

Theorem  A. 7.  Let  the  domain  0  be  Lipschitzian,  and  let  (A. 3. 8), 
(4. A. 5) , (4.4.6) ,  and  (4.5.1)  hold.  Then  there  exists  a  unique  solution 
u  £  V  to  the  variational  inequality  (4.5.2). 

4.5.1  Regularization  of  the  Functional  j(«)  The  particular  class  of 
problems  characterized  by  (4.5.2)  may  seem  to  be  of  limited  practical 
interest  since  it  depicts  a  rather  rare  situation  in  which  normal 
stresses  are  prescribed  on  a  surface  with  unknown  frictional  stresses. 
However,  we  shall  now  describe  a  regularization  of  this  problem  which, 
interestingly  enough,  can  depict  in  a  much  more  realistic  way  certain 
features  on  friction  mechanisms  observed  in  experiments  with  metals. 

We  begin  by  observing  that  the  functional  j(*)  in  (4.5.3)  is  not 
Gateaux  differentiable  at  the  origin.  This  is  not  surprising  because 
the  Euclidean  norm  | • |  in  1RN  is  not  differentiable  at  the  origin, 
nor  is  the  function  x  -*•  |x|  in  CR^  •  However,  the  source  of  the 
non-differentiability  of  j(*)  is  the  fact  that  this  particular  model 
of  friction  depicts  the  separation  of  the  sliding  and  full  adhesion 
portions  of  the  contact  surface  as  a  point  or  line.  Physical 
experiments  on  friction  on  metallic  surfaces  show  that  no  such  line  of 
separation  of  sliding  and  full  adhesion  exists;  rather,  there  is  a 
boundary-layer  between  regions  of  full  "stick"  and  "slip."  This,  in 
turn,  suggests  that  an  approximation  of  j(*)  differentiable  at  the 
origin  can  be  designed  which  leads  to  a  representation  of  such  a 
boundary  layer  but  which,  at  the  same  time,  can  be  made  arbitrarily 
close  to  the  functional  j(»)  in  some  sense. 


Toward  the  construction  of  such  an  approximation,  we  first 
consider  smooth  approximations  of  the  function  x  |x|  .  As  examples, 
consider  the  functions, 


I  x  -  ie 


( 


(x)  =  ^  x2  ;  -^-(x)  =<(  -r  x  if  lxl  <  c 


if  x  >  e 


-x  +  — £ 

V  2 


(4.5.5) 


-1  if  x  <  -e 


( 


<i>p(x)  = 


E+l  e  dx 


(-) 

e 


if  x  >  0 


(4.5.6) 


^  -(.2L)e  if  x  <  0 


and 

^(x)  =  -(x  tan_1^-i*n  (x2+e2));  ^E(x)«i  tan-1(l).  (4.5.7) 

x  £  2  dx  it  t 


It  is  easily  verified  that  all  of  these  approximations  converge  to  the 
absolute  value  function  as  e  ♦  0  .  The  most  popular  approximation  is 


the  first  one  which  exhibits  a  piecewise  linear  first  derivative 
However,  the  other  two  approximations  have  nonlinear  first  derivatives 
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It  is  also  worth  mentioning  at  this  point  that  an  analogy  of  the 
friction  problem  to  the  perfect  plasticity  can  be  realized  by  examining 

1  O 

at  the  graph  of  the  first  derivatives  of  <j>*  ,  *  and  $£  •  The 

first  resembles  the  stress-strain  curve  of  an  elastic-perfectly  plastic 
material  whereas  the  second  resembles  that  of  a  material  with  strain 
hardening. 

Another  approximation  which  is  reminiscent  of  plasticity  is  the 
function 


Here  e  and  £  are  two  positive  parameters  such  that  <f>£(x)  +  |x|  as  e  ♦ 
0  and  £  >  0  .  In  this  case,  we  have  a  sort  of  strain  hardening 
effect  after  the  slipping  (i.e.,  yielding). 
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To  introduce  similar  approximations  for  the  functional  j(*)  in 
friction  problems,  we  first  introduce  a  regularization  of  the  function 
4>(x)  =  | x | ,  where  |x|  is  the  Euclidean  norm  for  IRn  :  |x|  =  /x*x  . 

Then  its  approximation  similar  to  (A. 5. 8)  and  (4.5.9)  and  is 

given  by 


/ 


I  x  | 


.  e 

+  — x*x 
2 - 


£e2) 


if  | x  |  >  e 


=  < 


_I_ 

2e 

\ 


x*x 


if  | x  |  <  e  , 


(4.5.10) 


where 


f  (-j^-p  +  if  | x |  >  e 


d  X 


(x)*y  =  < 


(4.5.11) 


v  k* 


if  |x|  <  e 


Let  us  now  make  aa  estimate  of  the  difference  of  <{>E  and  <f> 


Note  that 
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|  |x*x  -  i(e  +  £e2) |  if  |x|  >  e 

k£(x)  -  4>(x)  |  = 

j-^-x’X  -  |x|j  if  I X |  c  e 

if  I  x  |  >  e 

if  |x|  <  e 


<  ^ixj2  +  c) 


(4.5.12) 


i  .e , 

|<t>e(x)  -  <f>(x)  l  <  ~-(  I x |  2  6  +  e).  (4.5.13) 

Thus,  as  the  parameters  t  and  e  tend  to  zero,  the 

approximation  of  the  nondif f erentiable  function  <KX)  = 

|x|  converges  to  <J>  at  a  rate  0(e  +  t) .  The  differentiable  function 
4>e  is  called  a  regularization  of  the  nondif f erentiable  function  4>  . 
As  shown  above,  there  are  infinitely  many  choices  of  such 

regularizations.  However,  for  definiteness,  we  shall  consider  only  the 
particular  choice  (4.5.10)  in  the  following  discussions.  Other  choices 
of  regularization  and  their  quality  in  computation  are  studied  in  e.g., 
KIKUCHI  [1982]. 


< 
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We  now  return  to  the  work  done  j(* )  by  the  friction  force  and 
its  regularization.  In  view  of  (4.5.13),  we  take  as  a  regularization 
of  this  functional. 


(4.5.14) 


We  easily  verify  that  with  this  approximation 


lje(v)  j(v)|  <  y11  ~'jtNn  0,  T  r  (  l^11 0 ,  A ,  f  re  +  (mes  rc)e), 


(4.5.15) 


where  mes  Tc  is  the  measure  of  the  boundary  rr  ,  and  |.Jn  _  is 

u  u,4,rc  ' 

the  norm  of  L4(TC)  defined  by 


V.rc  ‘  {' 


•/'pc(v*v)2ds^ 


(4.5.16) 


Note  that  we  have  used  here  the  facts  that  v-  £  L4(Tr)  for 

1  c* 

vi  ^  L1  if  n  is  Lip-Chitzian  for  ICICn,  and  that 

2  *  v  *  IT  *  Vt  +  • 

jl*-e?_l'£!E  jjJi*  Let  Q  be  Lipschitzian  and  let  "tN  €  L^(rP). 

IN  O 

Suppose  that  the  conditions  (4.3.8),  (4.4.5),  (4.4.6),  and  (4.5.1)  hold. 
Then  there  a  unique  solution  to  the  regularized  problem 

2?  :  a(“e>  v)  +  <Dje(ue),  v>  =  f2(v)  ,  Vv€vq  (4.5.17) 


for  a  given  pair  (e.6),  where  VQ  is  the  space  defined  in  (4.4.11), 


and 


>1 


(4.5.18) 


<Dje^e)*  l>  =  ‘  lTds  * 


Furthermore,  the  sequence  of  the  solutions  {u£}  to  the  regularized 
problem  (4.5.17)  converges  to  the  solution  u  of  the  friction  problem 
(4.5.2)  as  €  and  t  tend  to  zero.  Indeed,  there  exists  a  positive 
constant  C  ,  independent  of  e  and  £,  such  that 


II  ue  -  ull  l  <  C(/F  +  /?), 


(4.5.19) 


Before  proving  the  above  assertion,  let  us  note  that  the 
regularized  problem  (4.5.17)  is  equivalent  to  the  problem 


ue  €  V  ;  a(ue,  v  -  uE)  +  j£(v)  -  je(ue)  > 
Vv  €  V  , 


(4.5.20) 


since  j  (•)  is  a  convex  continuous  and  differentiable  functional  on 


Proof  of  Theorem  4.8.  Because  of  the  equivalence  of  (4.5.17)  to 
(4.5.20)  the  existence  of  a  unique  solution  u^  follows  from  reasoning 
similar  to  that  used  to  d  rive  (4.5.2).  Thus,  it  suffices  to  merely 
verify  the  inequality  (4.5.19). 

Applying  (4.5.2)  and  (4.5.20),  we  have 


a(u-ue,  u-ue)  <  j£(u)  -  j(u)  -  je(ue)  +  j(u£) 

<  c»-ptNll0  p  «lul§  4  r  +  «uelg  4  r  )£  +  (mes  rc)c) 


Because  of  (4.4.5),  the  solution  * ue  11 0  4  T  unif°rn,ly  bounded  in  £ 

*  *  C 

and  e.  We  can  thus  conclude  that  (4.5.19)  holds  for  sufficiently  small 
£  and  e  .  0 


Let  us  now  investigate  in  more  detail  the  physical  meaning  of  the 
regularization  method.  Toward  this  end,  the  variational  form  (4.5.17) 
is  considered.  If  the  generalized  Green  theorem  introduced  by  Aubin  ll979 
,  Chapter  13]  is  applied,  and  if  the  normal  and  tangent  stresses  oeN 
and  aej  resulting  from  the  displacement  field  ue  are  well-defined 
2 

in,  e.g.,  L  (Tq),  then  we  have  the  following  characterization  of  the 
solution  u£  of  (4.5.17): 


<«ij<Sc» 


f±  in  0, 


(4.5.21) 


in  the  sense  of  distributions, 

=  E  on  rF  •  °eN  =  tN  +  on  rC  > 


and 


-  3*e 

?eT  "  utN  —  ^~eT^  on  rC  * 
di?eT 


(4.5.22) 


(4.5.23) 


Here  o£l  =*  ®ij(ue)nj  ,  oeN  -  oe  •  N  ,  and  oeT  =  o £  -  oeNN. 

Consider  the  equilibrium  of  the  tangential  stress  on  the  contact 
surface  T^.  If  the  relation  (4.5.11)  is  introduced  into  (4.5.23),  we 
have  the  friction  stresses 


pt.  ( - -+  £u£T)  if  |lleTl  >  C 


•'  lu 


eT1 


2cT 


~  HgT 


if  lueTl  <  e 


(A. 5.24) 


on  the  surface  . 

It  is  clear  that  the  inequality 

if  |ufTl  >e 

if  |ueTi  <  e 

holds.  This  means  that  if  £  -*■  0,  the  friction  stress  never  exceeds 
the  value  -ptN  which  indicates  the  Coulomb  law  of  friction. 
Furthermore,  (4.5.24)  yields 

ueT  =  e  h^N2eT  for  '-eT1  <  e  * 

Thus,  by  passing  to  the  limit  e  >  0,  the  stick  portion  of  the  contact 
surface  is  identified  with  the  set  Se  =  {x  €  :  |ucj(x)|  <  c}. 

A  simple  spring  model  of  this  particular  friction  mechanism  is 
depicted  in  Fig.  4.9.  Note  that  the  initial  tangential  stiffness  at 
the  contact  surface  is  1/e  and  that  some  tangential  displacement  6  is 
reached  before  the  spring  of  stiffness  £  ,  representing  a  hardening 


I  T 1  <  5 


bt^(l  T  ^In^'pl) 


v  mtN 


effect,  is  activated. 
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4.5.2  Finite  Element  Approximations. 


We  shall 


now 


examine 


i 

i 

i 

•' . 

. ■ 

►  . 

s 


approximations  of  the  regularized  problem  (4.5.17)  by  finite  element 
methods.  Applying  the  notations  and  conventions  of  the  previous 
section  (recall  (4.4.28)  -  (4.4.37)),  the  regularized  form  (4.5.17)  is 
approximated  by 

u£  €  Vh  :  a(u^,  vh)  +  I(Dje(Ug) *v![i)  =  f2(vh)» 

Vvh€vJ,  (4.5.25) 

where  the  term  I(Dje(u^)«Vj)  is  the  quadrature  rule  for  the  friction 
term  <Dje(uj}),  v^>  similar  to  the  penalty  terra  in  (4.4.31).  It  is 
clear  that  the  finite  element  equation  resulting  from  (4.5.25)  is 
nonlinear  because  of  the  gradient  Dj^u^)  .  Thus,  to  solve  (4.5.25) 
we  must  generally  resort  to  some  successive  iterative  scheme.  Because 
of  the  nature  of  the  nonlinearity  in  Dje(iJe)>  it  is  natural  to 
consider  a  sequence  of  linearized  problems  which  employ  the  recurrence 
formula, 


h 


wVei^T  + 


i-U» 


-eT 


,  i-1 


-)if 


.  i-1  h  , 
1  ~eT ' 


>  c 


djJ^hS)  - 


(4.5.26) 


^HeT 


if  li_1H^rl  <  e 


eT1 


Then  the  solution  u]}  in  (4.5.25)  is  achieved  as  the  limit  of  the 
sequence  {*Ug}  8UC^  that 


b 
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iUg  €  vh  :  aC1^,  vh)  +  I(Djg(iUg),Vj) 

■  f2(vh).  Vvh  €  vjj, 


(4.5.27) 


for  i  ■  1,2 . m,  which  evolves  from  a  proper  initial  guess  °Ug. 

The  convergence  of  the  iterative  scheme  (4.5.27)  is  more  delicate 
than  the  case  for  the  penalty  formulation  (4.4.33)  since  the  term 

Utjj1  *u£.j.  /  | i_ ^HcT I  destroys  the  monotonicity  of  the  form  (4,5.27), 

unless  the  parameter  t  is  sufficiently  large.  Since  j£  is  convex, 
the  operator  constructed  by  the  quasilinear  form  a(»,»)  +  <Dje (»),•>  is 
strongly  monotone  on  V  x  Vq.  Thus,  it  is  possible  to  establish  the 
existence  of  a  convergent  subsequence.  However,  this  does  not  imply 
that  the  form  a(*,«)  +  <Dj*(»)»*>  is  monotone.  This  possibly  lead  to  a 
sort  of  oscillation  about  the  solution  as  i  ».  Such  oscillations 

might  be  reduced  or  eliminated  by  taking  a  larger  value  of  e(or 

fi)  that  results  in  the  reduction  of  the  effect  of  the  term 
ptNA  u£j/ I  ~eTl  during  the  iteration  process. 


4.5.3  Convergence  of  the  Finite  Element  Method.  Suppose  that  the  unique 
solution  uj}  of  (4.5.25)  is  obtained  by  the  iterative  scheme  (4.5.27). 

It  is,  however,  noted  that  while  the  existence  of  a  unique  solution  to 

(4.5.25)  is  guaranteed,  it  might  not  be  realized  by  this  iterative 

scheme  (4.5.27).  We  shall  verify  the  convergence  of  the  regularization 

finite  element  approximation  by  taking  the  following  steps:  an  estimate 

of  ,uc~9eli  is  first  obtained,  where  is  the  solution  of  the 

problem 


flg  €  vh  :  a(flj?,vh)  +  <Dje(flg),vh>  - 


(4.5.28) 


Vvh  €  . 


Because  of  the  application  of  the  numerical  integration  to  evaluate  the 
quantity 


_  h 


the  quadrature  error  must  be  considered  at  first.  Then  applying  the 
result  in  Theorem  4.8  yields  the  estimate  Iflg-u^lj  in  terms  of  e 
and  £  ,  where  u*1  is  the  solution  of  the  problem 

uh  €  Vh  :  a(uh,vh)  +  j(vh)  -  j(uh)  >  f2(vh-uh>» 

Vvh  €  Vh  .  (4.5.29) 

The  final  step  is  the  estimate  of  lu^-ulj  of  the  solution  u  of 
(4.5.2). 

To  accomplish  the  first  step,  let  us  note  that  the  form  (4.5.25) 
is  equivalent  to 

uj?  €  V*1  :  a(u*?,vh)  +  *(ytN$e(v^))  -  I(ptN4>e(UgT)) 

>  Vvh  €  Vh  .  (4.5.30) 

Similarly,  (4.5.28)  is  equivalent  to 


€  vh  :  a(flj,vh)  +  je(vh)-je(a£)  >  f2(vh-oJ) t 
Vvh  €  Vh  . 


(4.5.31) 
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Lemma  4.6.  Suppose  that  the  quadrature  rule  I  satisfies  the 
condition 

|EI<fhSh)l  <  C6h2<a-»  '£h's-3/2,rc,8hls-l/2,4,rc  (4.5.32) 

for  piecewise  polynomials  contained  in  the  trace  of  the  functions  In 
Vh  on  the  boundary  I*.  Then  there  is  a  positive  constant  C 
independent  of  the  mesh  size  h  and  the  parameters  £  and  e  such 
that 

<  (1+fi)  ^"He's*  ,9j,s»ll^N,s-3/2>hS~1  (4.5.33) 

for  s  >  1  . 

Proof.  From  (4.5.30)  and  (4.5.31),  we  have 

a<He“9e»  <  iV^NMiJeT^  I 

■b|  E  j(  ptj^({ie(fleip)  )  |  . 

Without  loss  of  the  rate  in  the  estimate,  we  can  assume  that  ptjj  is 
also  a  piecewise  polynomial  (see  the  procedure  to  (4.4.62)).  Applying 
the  assumption  (4.5.32),  it  follows  from  (4.4.5)  that 

<  ClptNls_3/2(l+fi'2elg_1/2,4>rc 
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Thus  (4.5.33)  follows.  0 
The  second  estimate, 

I  Og “J}h 1 1  <  C(/e  + /?)  (4.5.34) 

follows  easily  from  the  same  arguments  used  in  the  proof  of  Theorem 
4.8. 

Lemma  4.7.  Suppose  that  the  interpolation  estimate  (4.4.49)  holds 
for  the  function  on  ft  and  on  T  .  Suppose  that 

u  €  Hs(ft),  u|r  €  H8-1/2(rc),  oT(u)  €  Hs-3/2(rc), 

c 

and  ptN  €  Hs"3/2(rc).  (4.5.35) 

Then  for  the  solutions  uh  of  (4.5,29)  and  u  of  (4.5.2  ),  we  have 
the  estimate 

luh-uli  <  C(lul )g)hS“*,  (4.5.36) 

where  C  is  a  constant  Independent  of  the  mesh  size. 

Proof.  Using  (4.5.2  )  and  (4.5.29)  we  obtain  the  estimate 

a(u-u^,u-u^)  <  a(u-u*l,u-v*1)  +  a(u,v-u^) 

+  J(v)-j(uh)  -  f 2(v-uh)  +  a(u,vh-u) 

+  j(vh)-j(u)  -  €  Vh,  X  t  V. 

Taking  v  *  u*1  yields 
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a(u-u^,u-u^)  <  a(u-u*1,u-v^')  +  a(u, v^-u) 

+  j(vh)-j(u)-f2(vh-u). 

The  characterization  of  the  solution  u  of  (4.5.2)  after  applying  the 
generalized  Green  theorem  lead) to  the  inequality 


a(u-uh,u-uh)  <  a(u-u^,u-v^)  +  /p  a^.(u) (v^-u^)ds 

c 

+  /p  —  I2t^^s 


Under  the  assumption  (4.5.35),  we  have 

ralu-u^lip  <  Mllu-u^llj  Hu-v^Ij 

+  “2T(H>"s-3/2,rc  1 Ht~~t " -s+3/ 2 , r  c 
+  ®^N*s-3/2,rc  *HT“'XT,-s+3/2,rc 

Here  we  have  used  the  fact  that  Iv^MujI  <  | uT I  •  Applying  the 
interpolation  estimates  (4.4.49),  we  obtain 

mlu-uhlj  <  CMBu-uhl1hs“1lulg 

+^,2T^H^,s-3/2,rc  +  ,yIN,s-3/2,rc)ch2(s  1),H,s-i/2,rc 

The  estimate  (4.5.35) #is  then  obtained  by  using  Young's  inequality.  □ 
Combining  the  above  results  we  arrive  at  the  estimate  of  the 
approximation  of  the  regularization  finite  element  method  (4.5.25) 

Theorem  4.9.  Under  the  assumptions  (4.4.5),  (4.4.6),  (4.4.49), 
(4.5.32),  and  (4.5.35),  we  have 
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luj-u^  <  C^e  +  ft  +  hS_1>liJls  (4.5.37) 

where  C  is  a  proper  positive  constant  independent  of  h,  t,  an  e.  □ 

If  the  four  node  quadrilateral  element  with  the  trapezoid  rule  for 
the  quadrature  "1"  is  applied,  an  analysis  similar  to  that  leading  to 
Lemma  4.4  shows  that  the  estimate  (4.5.32)  reduces  to 

|E1(ytJj|u|}Tl2)|  <  ch2|»it}»1/2>rcluJTl23/2,4,rc  (4.5.38) 

where  ytj  is  the  interpolation  of  the  ytjj  on  Tq.  Thus,  if  s=*2,  the 
error  estimate  becomes 

luj}-uli  <  C(/F  +  /t  +  h)lul2.  (4.5.39) 

4.5.4  Numerical  Examples.  We  fehall  describe  three  example  problems  for 
friction  contact  problems.  The  first  is  a  footing  problem  in 
foundation  engineering  depicted  in  Figure  4.10. 

Suppose  that  the  foundation  is  composed  of  an  isotropic  elastic 
material  with  Young's  modulus  E=*5xl0  gK/m  and  Poisson's  ratio  0.3 
for  the  shaded  area,  the  remaining  portion  begin  a  second  elastic 
material  which  has  a  Young's  modulu  of  E=2xl0°gK/m  and  Poisson's 
ratio  of  0.3  for  the  other  part.  The  size  of  the  foundation  is  8m  x 
16m,  and  the  friction  coefficient  y  between  the  footing  and  the 
foundation  is  assumed  to  be  yB0.3. 

A  finite  element  model  is  obtained  by  64  four  -  node  elements  as 
shown  in  Figure  4.10  for  the  case  in  which  the  footing  is  pushed 
and  rotated  so  that  the  bottom  line  of  the  footing  lies  along  the 
plane  y  ■  -0.1  +  0.001(x  -  6).  That  is,  the  left  edge  of  the  footing 


is  indented  10  cm  into  the  foundation,  and  the  tangent  is  0.001.  For 
the  choice  of  e  =  10~^/E  and  £=0,  we  can  obtain  the  result  shown  in 
Fig.  4.11  in  which  three  nodal  points  from  the  right  edge  of  the 

footing  are  sliding.  Because  of  the  singularities  on  both  edges,  the 

contact  pressure  is  very  high  on  those  points.  It  is  noted  that  no 
oscillatory  tangential  stress  along  the  bottom  of  the  footing  is 

observed  because  of  the  relatively  course  finite  element  mesh  used. 

The  second  example  is  a  Hertzian  problem  with  friction.  A 
circular  rigid  cylinder  is  idented  into  the  deformable  body  whose 

Young's  modulus  E  is  1000  gK/cm^  and  Poisson's  ratio  v  is  0.3.  Under 
the  assumption  of  plane  strain,  the  problem  is  solved  by  36  four-node 
elements  with  the  depth  of  indentation  d=0.4  cm.  The  radius  of  the 
cylinder  is  taken  to  be  R=4  cm,  and  the  size  of  the  deformable  body  is 
8  cm  x  4  cm.  Using  the  symmetry  of  the  problem,  only  half  portion  of 
the  body  is  analyzed.  Computed  results  for  the  case  y=0.3,  e~2.6  x 

10'6,  and  £=0,  are  shown  in  Fig.  4.12.  The  deformed  configuration  is 
also  given  in  Figure  4.12.  Two  nodes  under  the  circular  rigid  punch  are 
sliding,  and  the  computed  contact  pressure  and  the  friction  stress  are 
plotted  in  Figure  4.13. 

The  last  example  is  an  axisymmetric  problem  for  an  annular  punch. 

The  computed  results  are  compared  with  the  analytical  solution  of 

Shibuya,  Koizumi,  and  Nakahara  [1930]  for  the  case  of  full  adhesion.  The 
deformed  configuration  and  the  contact  stresses  are  shown  in  Fig.  4.J4 
and  4.15,  for  the  case  that  E*1000  g  Kg/cm^,  v*0.3,  and  y*0.3.  The 

size  of  the  domain  Q*2  cm  x  1  cm  by  applying  the  axisymmetry.  The 

domain  is  discretized  by  200  four-node  elements.  We  here  need  not  to 


Distance  along 

60 ttorrv  of  footing 


Figure  4.11  Normal  Stress  on  the  Contact  Surface  (Pressure 
Field  on  the  Bottom  of  the  Footing) 


use  the  non-zero  regularity  parameter  £  .  If  e=10”^/E  is  taken,  the 
iterative  scheme  is  stable,  and  converges  monotonically . 

A. 6.  CONCLUDING  REMARKS 

Starting  from  the  contact  conditions  for  both  large  and  small  deformations, 
the  mathematical  description  (A. 2. 32)  of  the  friction  contact  problem  is  derived 
using  the  Coulomb  friction  law.  Under  the  assumptions  of  small  deformation  and 
quasi-static  motion,  the  set  of  dynamical  equations  (A. 2. 32)  is  reduced  to  the 
incremental  form  (A. 2. 35)  in  the  pseudo-time  t  and  its  increment  At.  The 
variational  statement  corresponding  to  the  incremental  form  (A, 2. 35)  is  derived 
by  using  a  quasi-variational  inequality  (A. 3. 5),  the  existence  of  solution  of 
which  is  proved  only  for  the  case  of  small  friction  coefficient  p. 

In  order  to  obtain  an  iterative  scheme  to  solve  the  quasi-variational 
inequality  that  represents  the  incremental  form  of  friction  contact  problems, 
two  special  cases  are  studied  in  detail.  The  first  case  is  the  one  with  the 
prescribed  tangential  stress  that  is  possibly  caused  by  friction.  Nonlinearity 
then  arises  from  the  unilateral  contact  in  the  normal  direction,  and  is  resolved 
by  applying  the  exterior  penalty  method  (A. A. 8).  Conditions  for  the  convergence 
of  (A. A. 8)  as  the  penalty  parameter  e  tends  to  zero  to  the  corresponding 
Lagrange  multiplier  formulation  (A. A. 26)  are  given  in  Theorem  A. 3.  The  finite 
element  approximation  of  the  exterior  penalty  formulation  is  obtained  using 
four-node  quadrilateral  isoparametric  elements,  and  its  convergence  analysis  in 
Theorem  A. 5  is  established  for  both  parameters  h  and  e,  where  h  is  the 
representative  mesh  size  of  the  finite  element  model. 

The  second  special  case  corresponds  to  problems  with  prescribed  normal 
stresses  on  the  known  contact  surface.  In  this  case,  if  the  normal  stress  is 
smooth  enough,  the  variational  formulation  is  given  by  a  variational  inequality 
of  the  first  kind  as  in  (A. 5. 2).  Because  of  the  non-differentiability  of  the 
functional  j,  that  represents  the  work  done  by  friction  forces,  a  regularization 
scheme  must  be  introduced  to  obtain  a  differentiable  functional  j£.  The 
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existence  of  a  unique  solution  to  the  regularized  problem  (A. 5. 20)  is  shown,  and 
the  convergence  of  the  regularization  is  proved  in  Theorem  8.  It  is  noted  that 
the  convergence  of  the  tangential  stress  o£j  has  not  been  analyzed  in  this 
study.  The  regularized  problem  (A. 5.20)  is  approximated  by  four-node  elements. 
Convergence  of  the  finite  element  approximation  (A. 5. 30)  is  established  in 
Theorem  A. 9,  under  the  assumption  that  the  iterative  scheme  (A. 5. 27)  gives  a 
convergent  result. 


5.  NON-CLASSICAL  FRICTION  LAWS 


The  use  of  the  classical  Coulomb  law  of  friction  in  the  formulation  of 
contact  problems  in  elasticity  leads  to  both  physical  and  mathematical 
difficulties;  the  former  arising  from  the  fact  that  this  law  provides  a  poor 
model  of  frictional  stresses  at  points  on  metallic  surfaces  in  contact  and 
latter  due  to  the  fact  that  the  existence  of  solutions  of  the  governing 
equations  can  be  proved  only  for  very  special  situations.  In  the  present 
article,  non-classical  friction  laws  are  proposed  in  an  attempt  to  overcome 
both  of  these  difficulties.  We  consider  a  class  of  contact  problems  in¬ 
volving  the  equilibrium  of  linearly  elastic  bodies  in  contact  on  surfaces 
on  which  nonlocal  and  nonlinear  friction  laws  are  assumed  to  hold.  The 
physics  of  friction  between  metallic  bodies  in  contact  is  discussed  and 
arguments  in  support  of  the  theory  are  presented.  Variational  principles 
for  boundary  value  problems  in  elasticity  in  which  such  nonlinear  nonlocal 
laws  hold  are  then  developed.  A  brief  discussion  of  the  questions  of 
existence  and  uniqueness  of  solutions  to  the  nonlocal  and  nonlinear  problems 
is  given. 

5.1  Introduction 

In  1781,  the  French  engineer  C.  A.  Coulomb  published  his  "Th£orie 
des  Machines  Simples"  in  which  he  presented  his  celebrated  law  of  friction. 
This  work  earned  him  a  double  prize  from  the  French  Royal  Academy 
of  Sciences  in  1785.  The  classical  Coulomb  law  of  static  dry  friction, 
of  course,  asserts  that  Aelative  A tiding  between  too  bodies  i n  contact 
along  plane  suAfiaces  will  occua  w hen  the  net  sheaA  fioA.ce  panallel  to 
the  plane  Aeaches  a  critical  value  pAopoAtional  to  the  net  nonmal  fioAce 


psieAting  the  tuoo  bodieA  together..  The  constant  o  f  p\opontionaJLLty  iA 
catted  the  coefficient  of  faction. 

It  can  be  argued  that  as  a  basis  for  contact  problems  in  the  theory 
of  elasticity.  Coulomb's  law  is  not  acceptable  from  either  a  physical 
or  a  mathematical  point  of  view.  From  the  purely  physical  side,  it  has 
been  recognized  for  many  years  that  Coulomb's  law  is  capable  of  describ¬ 
ing  only  friction  effects  between  effectively  rigid  bodies  and  gross 
sliding  of  one  body  relative  to  another.  Indeed,  it  is  clear  that  Coulomb 
himself  never  inteded  that  his  law  be  applied  pointwise  in  boundary- 
value  problems  in  elasticity;  the  foundarions  of  continuum  mechanics, 
particularly  the  concept  of  stress  and  the  equations  of  linear  elasto- 
statics,  were  only  fully  developed  many  decades  after  Coulomb  proposed 
his  law,  and  the  first  successful  formulation  of  a  contact  problem  in 
elasticity  came  over  a  full  century  after  Coulomb’s  work.  From  the 
mathematical  point  of  view,  it  is  known  (see  DUVAUT  [1980J  and  also 
DUVAUT  and  LIONS  [1976])  that  if  Coulomb's  law  is  applied  pointwise  in 
contact  problems  involving  linearly  elastic  bodies,  then  the  contact 

stress  o  developed  normal  to  the  contact  surface  is  ill-defined, 

✓  J 

Except  for  some  very  special  cases  (e.g.,  NECAS,  JARUSEK  and  HASLINGER 
[1980])  the  fundamental  question  of  existence  of  solutions  of  the  friction 
problem  is  open  (for  other  related  open  questions,  see  DUVAUT  and 
LIONS  [1976]). 

There  are  several  aspects  of  actual  friction  phenomena  between 
metallic  bodies  that  suggest  alternative  friction  laws  which  represent 
a  marked  departure  from  the  classical  formulations,  First,  we  mention 
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the  obvious  nonlocal  character  of  the  mechanism  by  which  normal  stresses 
are  distributed  on  machined  contact  surfaces.  These  stresses  are  trans¬ 
mitted  over  junctions  formed  by  deformed  asperities  and  are  not  concen¬ 
trated  at  isolated  points  on  the  contact  surface.  Second,  we  note' that 
upon  the  application  of  loads,  experiments  show  that  there  always  exists 
a  small  tangential  displacement  of  points  on  the  contact  surfaces  due 
to  the  elastic  and  elasto-plastic  deformation  of  these  junctions; 

"sliding"  occurs  when  these  junctions  are  actually  fractured.  Since 
these  junctions  can  be  recovered  upon  a  quasi- static  reversal  of  loads, 
the  actual  "adhesion-sliding"  friction  mechanism  is  highly  nonlinear 
and  depends  upon  the  elasto-plastic  properties  of  the  metal  oxide  and 
contaminant  film  on  the  contact  surfaces. 

Independent  of  the  nonlinear  character  of  local  friction  phenomena, 
there  are  also  mathematical  reasons  to  expect  that  a  nonlocal  friction 
law  might  lead  to  a  more  tractable  theory.  Recently,  DUVAUT  [1980]  pub¬ 
lished  a  brief  note  in  which  he  observed  that  the  source  of  difficulties 
in  establishing  an  existence  theory  for  Signorini’s  problem  with  Coulomb 
friction  was  the  lack  of  smoothness  of  the  normal  contact  pressure  0n  . 
Therefore,  by  considering,  instead,  a  proper  mollification  of  the  normal 
stresses  on  the  contact  boundary  he  was  able  to  state  that  the  principal 
obstacles  in  the  way  of  deriving  an  existence  and  uniqueness  theory  for 
contact  problems  with  friction  could  be  overcome.  One  interpretation  of 
such  a  smoothing  of  the  contact  pressure  is  as  result  of  nonlocal  effects 
arising  from  micromechanical  phenomena  taking  place  on  the  contact  regions. 

In  the  present  study,  we  propose  nonlinear,  nonlocal  friction  laws 
for  contact  problems  involving  linearly  elastic  bodies  and  we  present 
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variational  principles  for  contact  problems  in  elastostnt ics  in  which 
these  laws  hold.  Roughly  speaking,  a  nonlocal  friction  law  proposes 
that  impending  motion  out  a  point  oft  contact  bztiozzn  tioo  dzhormablz  bodies 
mM  occ.uk  u‘hzn  the  ihzar  6tru&  at  that  point  reaches  <t  value  propor¬ 
tional  to  a  lozightzd  measure  oh  the  nomal  6trzMzt>  in  a  neighborhood 
oh  the  point .  The  character  of  the  effective  local  neighborhood  and 
the  manner  in  which  neighborhood  stresses  contribute  to  the  sliding  con¬ 
dition  depends  upon  features  of  the  microstructure  of  the  materials  involved. 

While  such  nonlocal  laws  do  lead  to  a  mathematically  tractable 
theory,  they  still  do  not  capture  the  effects  of  the  tangential  elastic- 
plastic  deformations  of  the  contact  junctions  mentioned  earlier.  To 
accomodate  such  effects,  we  present  a  further  amendment  which  provides 
for  small  but  nonzero  elastic  tangential  displacements  at  the  contact 
surface  for  tangential  stresses  below  a  certain  critical  level.  For 
shear  stresses  at  or  near  this  critical  level,  substantially  larger 
motions  can  occur  which  effectively  represent  large  tangential  motions 
such  as  sliding.  This  critical  value  may  be  proportional  to  a  weighted 

measure  of  the  normal  stresses  in  a  neighborhood  of  the  point  on  the 
contact  surface. 

An  interesting  feature  of  our  results  is  that  these  non— conventional 
friction  laws  are  given  in  terms  of  three  positive  material  parameters: 

V,  p,  and  c.  The  parameter  v  is  the  coefficient  of  friction,  although 
its  actual  interpretation  is  somewhat  more  complex  than  that  of  classical 
mechanics.  The  parameter  p  quantifies  the  nonlocal  character  of  the 
response  -  for  p  =  0  a  fully  local  law  is  obtained.  Finally,  c  is  a 
measure  of  the  tangential  stiffness  of  the  elastic-plastic  junctions  on 


the  contact  surface;  the  case  e  =  0  corresponds  to  a  fully  rigid  response 
full  adhesion  or  full  sliding  of  contact  surfaces.  Thus,  by  allowing  p 
and  €  to  approach  zero,  we  can  recover  the  classical,  local,  pointwise 
formulation  of  contact  problems  in  elast istatics  based  on  Coulomb's  law. 

Following  this  introduction,  we  give  a  brief  account  of  the  physics 
of  friction  as  well  as  a  justification  of  specific  friction  models. 

Several  variational  principles  for  boundary  value  problems  in  elasticity 
in  which  nonlocal  and  nonlinear  laws  are  assumed  to  hold  are  derived  in 
Section  3.  The  results  of  some  theoretical  studies  of  these  principles 

are  summarized  in  Section  4.  Our  results  include  conditions  sufficient 
to  guarantee  the  existence  and  uniqueness  of  solutions  to  the  variational 

problems  as  well  as  results  which  establish  the  asymptotic  behavior  of 
solutions  as  parameters  p  and  £  tend  to  zero. 


Ill 


5.2  A  Basis  for  Nonlocal  and  Nonlinear  Friction  Lays. 

5.2.1  Micromechanics  of  Friction.  Friction  phenomena  have  been  the  subject 
of  considerable  experimental  and  theoretical  research  over  -he  last  30 
years  and  its  study  is  a  popular  and  important  aspect  of  modern  mechanical 
engineering  design.  A  standard  reference  is  the  monumental  two-volume 
treatise  by  BOWDEN  and  TABOR  [1950  and  1964]  ;  more  concise  accounts  can 
be  found  in  standard  texts  on  the  subject  (e.g.  RABINOWICZ  [1965]). 

To  understand  friction,  one  must  first  appreciate  the  role  of  the 
microstructure  of  the  materials  involved.  Consider  an  experiment  in  which 
two  metallic  bodies  are  placed  in  contact  along  two  apparently  machined 
flat  surfaces.  At  microscopic  levels,  specifically  at  magnifications  of 

lOOOx  to  5000x,  machined  metal  surfaces  are  seen  to  be  not  smooth  homo¬ 
geneous  planes,  but  rough  contours  with  numerous  irregularities  which  are 
large  compared  with  molecular  dimensions.  We  refer  to  these  deviations 
from  the  plane  as  OApe/tctceA. 

When  we  press  together  two  surfaces,  actual  contact  initially  occurs 
only  at  the  peaks  or  summits  of  the  asperities.  Large  areas  of  the  sur¬ 
faces  are  separated  by  a  distance  which  is  large  compared  with  the  range 
of  molecular  action,  so  that  these  gaps  in  the  surfaces  are  completely 
separated  and  have  no  interaction  with  one  another.  The  load  is,  there¬ 
fore,  initially  supported  at  the  tips  of  the  asperities;  the  area  of 
contact  is  extremely  small,  and  the  pressure  at  the  points  of  contact, 
even  for  lightly  loaded  surfaces,  is  high.  Plastic  deformation  of  the 
tip  of  the  asperities  occurs  at  small  loads  while  the  bulk  of  the  under¬ 
lying  metal  deforms  elastically.  As  the  normal  load  is  further  increased, 


the  asperities  deform  and  fracture  with  the  result  that  the  local  load 
is  distributed  over  an  area  surrounding  each  deformed  asperity.  At  this 
stage,  each  asperity  has  been  flattened  and  the  local  contact  forces 
are  distributed  over  a  neighborhood  of  the  asperity.  The  ACaZ  aAta  0^ 
contact  (as  opposed  to  the  apjMACnt  area)  is  therefore  the  sum  of 

the  areas  of  all  the  surface  irregularities  which  are  touching  .and  which 
support  the  load. 

It  is  often  assumed  that  the  local  plastic  yield  pressure  is 

nearly  constant  and  is  comparable  to  the  indentation  hardness  of  the 
metal.  Under  these  circumstances,  the  real  area  of  contact  for  any  one 
asperity  bearing  a  load  is  =  N^/p^;  for  the  assembly  of  the 

asperities,  the  real  contact  area  is 


(5.2.1) 


where  N  is  the  total  normal  force  pressing  the  surfaces  together. 

The  real  area  of  contact  is,  thus,  proportional  to  the  load  and  indepen¬ 
dent  of  the  size  of  the  surfaces.  Over  these  regions  where  intimate 
contact  occurs,  strong  adhesion  and  welding  of  the  metal  surfaces  takes 
place  and  the  specimens  become,  in  effect,  a  continuous  solid.  We  refer 
to  these  regions  as  junction 6. 

Under  most  working  conditions,  metal  surfaces  are  covered  by  a 


thin  film  of  oxide,  water  vapor  and  other  absorbed  impurities.  The  shear 
strength  of  these  junctions  can  be  strongly  dependent  upon  the  shear 
strength  of  these  surface  films.  In  particular,  it  is  the  shear  strength 


of  this  layer  of  oxide  and  impurities  that  determines  the  coefficient  of 
friction  and  not,  in  general,  the  shear  strength  of  the  parent  metals. 

The  application  of  a  tangential  force  T  creates  a  tendency  for 
the  two  bodies  to  slide  relative  to  each  other.  The  contact  pressure 
must  then  decrease  (since  it  was  near  or  equal  the  plastic  yield  stress). 

Some  microscopic  motion  will  occur  (junction  growth)  even  when  T  is 
small.  If  T  is  steadily  increased,  a  value  of  sufficient  magnitude  to 
fracture  the  contaminant  films  is  eventually  reached  and  gross  sliding  occurs 
It  is  customary  to  set  the  ratio  of  the  magnitude  of  T  at  which  sliding- 
occurs  to  the  net  normal  force  N  equal  to  the  coefficient  of  friction 
V  .  If  s  is  the  average  shear  strength  of  the  interface,  it  follows 
that  approximately 

T  -  A  s  (5.2.2) 

r 

i.e.,  T  is  independent  of  the  apparent  area  of  contact  (since  is). 

Then,  substituting  for  »  we  obtain 


i.e.,  T  is  directly  proportional  to  the  load  (or  v  is  independent  of 
the  load). 

As  an  idealization  of  the  contact  surfaces,  one  may  assume  that  the 
asperities  are  superposed  upon  the  surface  of  spherical  protuberances  with 
a  larger  radius  of  curvature.  It  is  then  possible  to  consider  that 
although  each  individual  asperity  at  the  interface  will  deform  plastically 


the  deformation  of  each  spherical  protuberance  will  be  elastic.  This 
idealization  has  two  purposes:  first,  the  area  of  real  contact  still 
remains  approximately  proportional  to  the  load,  although  the  overall 
deformation  is  elastic  (see  ARCHARD  [1957]);  second,  it  allows  us  to 
treat  each  region  of  contact  as  being  roughly  circular  (see  microphoto¬ 
graphs  in  BOWDEN  and  TABOR  [1964,  p.  71])  and  to  regard  the  contact  pres¬ 
sure  as  being  essentially  symmetric,  attaining  its  maximum  magnitude  at 
the  center  of  the  circle  of  contact,  in  a  manner  consistent  with  the 
well  known  analysis  of  MINDLIN  [1949]. 

We  emphasize  that  the  junctions  through  which  loads  are  trans¬ 
mitted  from  one  body  to  another  are  not  rigid;  indeed,  they  are  com¬ 
posed  of  a  deformable  composite  of  metal,  metal  oxide,  and  surface 
contaminant  that,  for  our  purposes,  can  be  assumed  to  be  elasto-plastic 
or  nonlinearly  elastic.  Several  researchers  have  actually  measured  the 
tangential  micro-displacements  that  occur,  in  friction  experiments  on 
metals,  prior  to  gross  sliding  of  the  surfaces,  and  we  mention  as  examples 
the  works  of  JOHNSON  [1955],  BOWDEN  and  TABOR  [1964]  and  RABINOWICZ  [1965] 
Figure  1  reproduces  a  typical  results  of  static  displacement  tests  of 
JOHNSON  [1955]  which  involved  the  contact  of  hard  steel  balls  with  the 
flat  end  of  a  hard  steel  roller.  Mirco-displacements  are  produced  by  an 
applied  shear  force  varying  progressively  from  zero  to  the  value  necessary 
to  produce  slip. 


5.2.2  A  Nonlocal  Friction  Law.  In  order  to  develop  a  basis  for  a  nonlocal 
friction  law  we  consider  here  the  two  simple  physical  models  shown  in 
Figs.  2  and  3.  Fig.  2,  a  thin  weightless  strip  A  of  length  21 


Disp'jcement  (micro -inches) 


Diameter  of  the  ball  * 
Normal  load  s 


0.  375  in. 
10.5  lbs. 


5.1. 


Measured  shear-tangential  displacement 
variations  after  the  experiments  of 
Johnson  119551. 


-t 
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X 
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a)  A  thin  strip  A  pressed  onto  an  elastic  block  B  with 
surface  asperities;  b)  the  shear  stress  distribution 
at  the  point  of  impending  motion 


Figure  5.3. 
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is  pressed  against  a  fixed  elastic  block  B  by  a  concentrated  normal 
force  N  applied  at  its  midpoint  P  ;  then  a  force  T  is  applied  to 
the  strip  and  is  increased  slowly  in  magnitude  until  motion  (sliding) 
of  the  strip  relative  to  the  block  occurs.  For  this  idealized  situation, 
designed  to  emphasize  the  local  character  of  the  classical  Coulomb  law 
of  friction,  we  are  interested  in  calculating  the  distribution  of  fric¬ 
tional  (shear)  stresses  between  the  strip  and  the  block  assuming  that 
Coulomb's  law  holds  and  that  the  coefficient  of  friction  V  is  given. 

Ke  consider  the  origin  of  a  coordinate  axis  x  along  the  length  of 
the  strip  to  coincide  with  the  midpoint  P  of  the  strip.  Then,  the 
only  point  on  the  contact  surfaces  between  bodies  A  and  B  at  which 
a  resistive  force  can  be  developed  is  at  an  isolated  point  beneath  P. 
Thus,  impending  sliding  is  reached  when  the  resisting  shear  o^,  is 
formally  given  by 

a^(x)  =  vN6(x)  ;  -Jt  <  x  <  l  (5.2.4) 


where  6  is  the  Dirac  delta  corresponding  to  the  point  source  at  the 
origin.  Of  course,  is  merely  the  symbolic  representation  of  the 

distribution 


<cT,$>  -  VN  <6,$>  =  VN<|><0) 


(5.2.5) 


for  all  test  functions  in  the  cliss  PC-£.£)  of  infinitely  differen-r 
tiable  functions  with  support  in  the  interval  (-1,1)  f  where  <<»^  denotes 
duality  pairing  on  distributions  and  test  functions  (i.e.,  the  action  of 


a  distribution  q  on  a  test  function  $  is  denoted  q(4>)  =  )• 

Alternatively,  6  can  be  interpreted  as  the  limit  of  a  6-scquence, 

{(i)  }  ,  U)  in  P(-£,£) : 

p  a<P  p 

i 

<•(0)  =  «(♦)  ■  H»  “  ♦  for  all  *  in  <5.2.6) 

P>0 

Then  we  have,  instead  of  (5.2.5) 

<0  ,♦>  =  v  lira  <N  (0  J>  for  all  <J)  in  P(-£,£)  (5.2.7) 

A  pK)  w 

Me  see  that  the  classical  pointwise  version  of  Coulomb's  law  must  be 
interpreted  in  the  sense  of  distributions  for  this  situation.  As  a 
typical  6-sequence,  we  mention: 

c  exp[p2/(x2  -  p2)],  |x|  <  p 

a)  (x)  =  <  (5.2.8) 

P 

°  ,  (x|  >  p 

A  more  realistic  model  of  friction  from  the  physical  point  of  view 
is  obtained  if  we  take  into  account  the  microscopic  aspects  of  the 
physics  of  friction  described  earlier.  Specifically,  the  contact  sur¬ 
face  of  body  B  will  present  asperities  deviating  from  a  smooth  plane. 

As  the  normal  force  N  is  gradually  applied,  these  asperities  are  grad¬ 
ually  deformed  and  broken  down  until  equilibrium  of  normal  forces  is 
reached.  The  normal  force  reaching  body  B  through  the  strip  A  must 
then  be  distributed  over  the  contact  area  of  the  deformed  asperity  as 
indicated  in  Fig.  3,  We  shall  now  assume  that  the  asperity's  finite 
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transmission  area 
5.2.8  keeping  p  = 
deformed  asperity, 
of  (5.2.7), 


is  accounted  for  by  using  the  6-sequence 
Pq,  Pq  being  the  radius  of  the  contact 
Since  N  =  N(x)  is  now  a  function,  we  h 


{tu  ) 

P 

area 
ave , 


of 

of  the 
instead 


O  (x)  =  VN(y)*0J  (x  -  y) 
P0 


(5-  2.9) 


where  *  denotes  the  convolution  of  the  two  functions.  Thus,  we  have 

arrived  at  a  friction  law  in  which  impending  motion  occurs  at  a  point  x 

on  the  contact  surface  when  the  shear  stress  at  that  point  reaches  a 

value  proportional  to  the  weighed  average  of  the  normal  stress  in  a 

neighborhood  of  the  point.  If  u  is  used  to  characterize  this  weight- 

P0 

ing  function,  then  the  neighborhood  is  a  circular  disc  of  radius  p^ 
centered  at  x  ,  the  maximum  weight  is  given  to  the  stress  intensity  at 
the  center  of  the  disc  (the  contact  area  of  the  deformed  asperity)  and 
exponentially  decreasing  weights  are  assigned  to  stress  intensities  as 
one  moves  from  the  center  of  the  neighborhood  outward  to  the  periphery 
of  the  disc. 

We  can  now  generalize  these  results  to  the  three-dimensional  case: 

let  u^.  denote  the  relative  tangential  component  of  displacement  of  a 

point  x  =  (x^,X2»x^)  on  the  contact  surface  between  two  deformable 

bodies  and  let  o  (u)  and  0_(u)  denote  the  normal  and  tangential 
n  ~  ~  I  ~ 

stresses  on  the  contact  surface  corresponding  to  the  displacement  field 


u 


Then 
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(u)  |  <  vS  (O  (u)) 
*  n  ~ 


impl ics 


'T 


0 


|o.f(u)|  =  vSpo<0n^>) 


implies  that  there  exists 

X  >  0  such  that  u_  »  -Xo_ 
—  ~T  -T 


\ 

>  (5.2.10) 

J 


where 


is  an  operator  mollifying  the  normal  stress  distribution;  e.g. 


S  (o(u))(x) 

po  n  *  - 


y  1)  (-on(u(y)))dy 


(5.2.11) 


where  x  and  y  are  points  on  the  contact  surface  T^,. 

Ve  mention  that  nonlocal  theories  for  other  classes  of  problems  in 
solid  and  fluid  mechanics  have  been  put  forth  by  ERINGEN  and  EDELEN  [1972]; 
a  detailed  account  of  this  work  can  be  found  in  ERIN’GEN  [1976], 

5.2.3  Model  of  Nonlinear  friction.  Both  Coulomb's  law  and  the  nonlocal  law 
given  in  5.2.10  depict  perfect  rigid-adhesion-sliding  conditions  on  the 
contact  surface:  they  assert  that  there  is  absolutely  no  motion  of  points 
of  one  body  relative  to  those  of  another  if  the  tangential  stress  on  the 
contact  surface  remains  below  some  critical  value  T  ;  but  when  this  limit 
is  reached,  unbounded  motions  can  occur,  the  ensuing  tangential  displace¬ 
ment  being  directed  opposite  to  the  tangential  stress  vector.  It  was 
pointed  out  in  the  Introduction  and  in  Section 5.2  that  in  physical  experi¬ 
ments  on  contact,  tangential  displacements  are  produced  by  any  nonzero  tan¬ 
gential  force  developed  on  the  contact  surface  since  elastic-plastic  defor¬ 
mations  of  the  junctions  will  always  accompany  the  application  of  tangential 
forces  (recall  Fig.  1). 

To  model  this  phenomenon,  we  shall  consider  a  family  of  nonlinear 
friction  laws  of  the  form 


whore  the  function 


♦eC> 


satisfies  the  following  conditions: 


1)  <p  is  a  continuous,  monotone,  real-valued  function  of 
the  non-negative  real  numbers  r  ,  depending  on  a  para¬ 
meter  £  >  0,  such  that 

0  <  4>e(r)  £  1  ,  r  >  0 

ii)  lira  <J>  (r)  =  1  for  all  r  >  0 
e  0  e 

iii)  lira  <}>  (r)  =  1  for  all  e  >  0 

r  ■+  oo  e 


\ 


>(5.2.13) 


J 


In  5.2.12,  T  is  a  non-negative  function  of  the  displacement  vector  u 

representing  the  critical  value  that  the  tangential  stress  cannot  exceed, 

i.e.,  x  =  via  (u) |  for  the  local  case  and  T  =  VS  (a  (u))  for  the  non- 
1  n  -  p  n  ~ 

local  case.  It  will  also  be  of  interest,  as  it  will  be  seen 
later,  to  consider  the  case  in  which  T  is  a  given  (known)  function  of 
the  position  vector  x  ,  thus  no  longer  dependent  upon  the  displacement  u. 
As  specific  examples,  we  select  for  <p^  the  following  two  functions: 


A)  <J>e(|uT|) 


(5.  2,14) 


8)  4>£(!ut|) 


tanh 


(5.  2.15) 


It  is  readily  seen  that  both  of  these  functions  satisfy  5.2.13.  The  first 
example  represents  an  "elastic-perf ectly  plastic"  type  response  in  which 
slipping  can  occur  only  after  a  tangential  displacement  |u^,|  >  e  .  On 


the  other  hand,  the  second  example  describes  a  situation  in  which  the  cri¬ 


tical  stress  is  approximated  asymptotically  as  |u^,|  <*>  .  Both  curves 

are  depicted  in  Fig.  4.  We  notice  that  these  are  not  the  only  possible 
choices  and  that  several  others  can  also  be  considered  (e.g.  arctan,  etc.). 
Also,  we  observe  that  the  slope  of  <f>c  at  the  origin,  i.e.,  the  derivative 
of  at  zero  displacement  equals  1/e.  Thus  E  provides  a  measure  of 

the  rigidity  or  stiffness  of  the  elastic-plastic  or  nonlinear  elastic  junc¬ 
tions. 

Finally,  we  wish  to  comment  on  the  combination  of  a  nonlocal  law  with 
a  nonlinear  law.  A  nonlocal-nonlinear  law  of  friction  will  be  of  the 
form (5. 2 . 12)with  T  replaced  by  VS^  (o^(u)),  i.e.. 


a  (u) 

I 


~VSp^an^~^ 


(5.2.16) 


where  is  of  the  form  given  in  (5.2.11).  If  we  allow  p  •*  0  maintaining 

e  >  0  fixed  we  obtain  a  local  nonlinear  law,  i.e.. 


rt 

~T 


(u) 


-V  |on(u)|  *eC|uT|) 


(5.2.17) 


On  the  other  hand,  if  we  allow  e  -*■  0  for  a  given  positive  p  we  recover 
the  perfect  rigid-adhesion-sliding  nonlocal  law  given  in  (5.2.10).  Finally 
if  both  e  -*■  0  and  p  ->  0  the  static  Coulomb's  law  (local)  for  unilateral 
contact  in  the  contact  surface  is  obtained: 

|oT(u)  I  <  V  |onCu)|  implies  uT  =  0 

|o^,(u)  j  =  V  Jo^(u)  |  implies  the  existence  of  X  >  0 

such  that  u„  =* 
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5.3  Variational  Principles  for  Nonlocal  Nonlinear  Friction 


5.3.1  Signorini  Problems  with  Nonlocal  Nonlinear  Friction.  Wg  consider 
here  the  Signorini  problem  of  contact  of  a  linearly  elastic  body  with  a 
rigid  foundation  on  which  nonlocal  nonlinear  friction  laws  hold. 

We  begin  our  analysis  by  considering  a  linearly  elastic  body  the 
particles  of  which  occupy  a  smooth  bounded  domain  ft  in  ,  N  =  2,3, 
with  open  interior  ft  .  The  boundary  T  of  the  body  is  assumed  to  con¬ 
sist  of  three  disjoint  parts,  T  ,  T  and  T  ,  where  T  and  T  are 
the  portions  of  the  boundary  on  which  the  displacements  and  forces  (trac¬ 
tions)  are  prescribed  respectively  and  is  the  candidate  contact  area; 

i.e.  T  is  a  portion  of  the  boundary  which  contains  the  material  surface 
which  comes  in  unilateral  contact  with  a  rigid  foundation  F  upon  the  ap¬ 
plication  of  loads  (see  Fig.  5).  The  external  forces  on  the  body  consist 
of  a  prescribed  body  force  field  of  intensity  f  per  unit  volume  and  of 
surface  tractions  of  intensity  t  per  unit  surface  area. 

We  shall  assume  that  is  perfectly  fixed,  so  that 

u  =  0  on  Tjj 

u  being  the  displacement  field.  On  fp  we  will  have 


ai] (:)oi  ■  'j 


where  a  (u)  is  the  stress  produced  by  u  and  n.  are  the  components 
ij  -  -  i 

of  the  unit  outward  normal  n  to  T.  Here  and  throughout  our  presenta¬ 
tion  Cartesian  index  notation  and  the  summation  convention  are  employed. 
Since  the  body  is  assumed  to  be  linearly  elastic,  Hooke’s  law  holds  so  that 


where  E 


are  the  usual  elastic  constants  of  the  material 


ijk£ 
3u^/3x^  and 


•Eijk£  =  Ejik£  =  Ek£ij  Eij£k 
1  <  i,  j  ,  k,  £  <  N 


Of  course,  particles  within  are  assumed  to  be  in  static  equili¬ 
brium  so  that 


°ij  fa)  » j  +  f  ±  =  0  in  ^ 

The  unilateral  motion  of  particles  of  the  body  on  the  material  sur¬ 
face  rc  is  constrained  by  the  presence  of  a  rigid  foundation  which  is 
a  given  distance  s  from  the  body  prior  to  the  application  of  loads.  Ma 
thematically,  this  constraint  is  represented  by  the  requirement  that  the 
normal  displacement  of  boundary  points  cannot  exceed  s  : 


u  •  n  <  s 


on  r. 


<5.3.1: 


If  u  •  n  *  s  then  contact  is  established,  while  u  •  n  <  s  indicates  the 
existence  of  a  gap  between  the  support  and  the  body.  Thus  5,3.1  repre¬ 
sents  a  non-penetration  condition.  If  contact  is  not  made  (i.e.,  if 


u  •  n  <  s)  then  the  normal  contact  pressure  O^Cu)  3  0  »  where 


Alternatively,  if  u  •  n  =  s  at  a  point  on  P  ,  then  o  must  be  non- 

~  C  n 

posit ive: 

u  •  n  =  s  implies  0^(u)  £  ® 

Thus,  the  unilateral  contact  conditions  on  the  contact  surface  are: 

u  *  n  -  s  <  0  ,  O  (u)  <  0 
—  n  - 

0  (n)(u  •  n  -  s)  -  0 
n  -  -v. 

Condition  0^(u)  (u  *  n  -  s)  -  0  signifies  that  the  pressure  can  only  be 
nor -zero  when  there  is  contact. 

There  remains  the  characterization  of  friction  on  the  contact  surface 
T  .  In  agreement  with  the  nonlinear  and  nonlocal  friction  laws  described 
in  the  previous  Section,  we  will  have 


\ 

If  u  •  n  <  s  ,  0^,(u)  =  0 

If  u  •  n  =  s  ,  /  on  (5.  3.2) 

gT(u)  -  -vsp(an(u))  |~y  y 

where  S  is  of  the.  for:n5. 3. 1J  and  <{,  is  a  nonlinear  function  of  the 
P  E 

type  in  5.2.13  or  hhlh  It  is  important  to  note  that  always 

loT(u) |  <  vSp(an(n)) 

since  (j>  (  |u.  |)  <  1  for  every  ;  >  0  and  any  u  (x)  ,  x  in  T  .  By 
computing  the  inner  product  of  both  sides  of  the  second  equation  0,3,2) 
with  u  ,  we  can  write,  equivalently. 
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c>T(u)  *  -T  +  vSp^°n^^  hTl  =  0 


Summarizing,  the  Signorlni  problem  with  nonlocal-nonlinear  friction 
consists  of  seeking  a  displacement  field  u  ,  which  satisfies  the  following 
system  of  equations  and  inequalities: 


1)  Equilibrium  equations 

(Eljkl  "k.n'-j  +  fi  *  0  ln  n 


(5.3.3) 


2)  Boundary  conditions 

a)  Prescribed  displacements 


ui  =  0  on  rD 


(5.3.4) 


b)  Prescribed  traction 


Eljk«Uk.l"j  '  C1  °"  PF 


(5.3.5) 


c)  Unilateral  constraint 


u  •  n  <  s  ,  a  (u)  <  0  , 
~  -  —  n  -  — 

a  (u)  (u  •  n  -  s)  =  0 
n  -  -  - 


on  r. 


(5.3.6) 


d)  Friction  condition 


aT(u)  =  -vSp(an(u))  <j>c(!uTJ)  -j^y  on  ^ 


(5.3.7) 


where 


S  (a  (u))(x)  =  [  u  ()x  -  y  1 )  C-cr  (u(y)))dj 
p  n  ~  j  p  -  n  ~ 


and  u  and  $  are  given,  respectively,  by  5.2.8  and  5.2.14 
P  ^ 

or  5.2.15. 
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5.3.2  Variational  Principles  for  Non  1  oca  1  Non  1  i  mar  Friction.  Re  1  a  t  ion  sli  i  p 
with  the  Classical  Problems.  We  now  introduce  variational  principles  for 
the  nonlocal-nonlinear  Signorini  problem  and  establish  the  relationship  be¬ 
tween  the  variatonal  and  the  classical  formulations. 

Using  the  notation  introduced  previously,  we  define 

V  =  space  of  admissible  displacements.  A  displacement  vector  v 
will  belong  to  V  if  and  only  if 

1)  v  =  0  on  Fp 

2)  v  produces  finite  (normalized)  strain  energy  in  the  sense 


that  the  norm 


(5. 3.8) 


M'v  ’  Vn'i,!  vij  ^ 


is  finite,  where  dx  -  dx^dx^-.-dx^  . 


K  =  subset  of  V  consisting  of  all  admissible  displacements  v  in 

V  for  which  v  •  n  _<  s  at  all  points  on  the  contact  surface  F  . 

a(u,v)  =  virtual  work  produced  by  the  action  of  stresses  a  (u)  on  strains 

ij 

caused  by  the  displacement  v 


’  Jn  Eiju  V.  V,  d* 


=  virtual  work  done  by  the  frictional  forces  on  the  displacement  v 


=  j  vsp(an(u))  U'E(lyT|)ds 


Here  ds  is  an  element  of  surface  area  on  T  and 
is  rbe  primitive  of  :  ~  ^  •  For  example. 
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Vlvjl) 


|vTl  -  c/2  if  |vT|  >  c 
|vT|2/2e  if  |vT|  <  e 


(5.  3.9) 


Find  an  admissible  displacement  vector  u  '  n  the  set  K 
such  that 


a(u,v-u)  +  j  (u,v)  -  j  (u,u)  _>  f(v-u) 


(5.3.11) 


for  all  admissible  displacements  v  in  K  . 

Inequality  5. 3. 11  is  a  statement  of  the  principle  of  virtual  work 
for  an  elastic  body  restrained  by  frictional  forces  of  the  type  in  5.3.2 
Note  that  this  characterizes  equilbriutn  configurations  by  an  inequality  rather 
than  an  equality  because  of  the  presence  of  the  unilateral  contact  constraint 

u  •  n  <  s  on  T_  .  We  also  notice  that  the  actual  contact  surface  de- 

<*•  ^  v 

pends  upon  the  solution  u  and  is,  therefore,  not  known  in  advance. 

Our  first  major  result  is  stated  in  the  following  proposition,  the 
proof  of  which  is  given  in  the  Appendix: 
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Proposition  3.1:  Let  u  be  a  sufficiently  smooth  solution  of  the 
Signorini  problem  with  nonlocal  nonlinear  friction  5.3.3  -  5.3.7.  Then 
u  is  also  a  solution  of  the  variational  inequality  5.3.11.  Conversely, 
if  u  is  a  solution  of  5.3.11 f  then  u  also  satisfies  the  system  5.3.3  - 
5.3.7  if  these  relations  are  interpreted  in  a  weak  or  distributional  sense. □ 


5.3.3  Other  Related  Friction  Problems.  Several  friction  problems  can  be 
obtained  as  special  cases  of  the  nonlocal  nonlinear  problem  discussed  in 
the  previous  Sections  by  allowing  p  *■  0  or  e  >  0  or  both  p,£  ->  0  or  by 
restricting  the  dependence  of  the  friction  functional  to  its  second  vari¬ 
able.  Results  similar  to  those  established  in  Proposition  3.1  can  also  be 

stated.  Moreover,  if  vS  (a  (u))  is  prescribed  on  T  ,  and,  hence,  is  in- 

p  n  —  C 

dependent  of  the  displacement  u  ,  the  static  friction  problem  thus  ob¬ 
tained  becomes  equivalent  to  a  constrained  minimization  problem  involving 
a  functional  representing  the  associated  potential  energy.  Thus,  we  have 
the  following  cases: 


Case  I  (p  =  0) .  The  friction  law,  which  is  now  of  a  local  type,  is 
given  by  5.2.17.  Thus  the  system  5.  3.3  -  5.3.6  together  with  condition 
5.2.17  will  produce  a  Signorini  problem  with  local-nonlinear  friction 
which  can  be  shown  to  be  equivalent  to  the  variational  principle  5.3.11 
with  the  functional  j  (*, *)  replaced  by 

P  ,  C. 


Vm5 


rc 
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Case  II  (c  =  0) .  For  this  case,  the  friction  law  is  given  by  5.2.10 
If  we  add  this  condition  to  the  system  5.3.3  -  5.3.6  we  will  obtain  a 
Signor  ini  problem  with  nonlocal  friction  which  may  be  seen  to  be  equiva¬ 
lent  to  the  variational  inequality  5.3.11  if  we  replace  j  (•  • )  in  5.3.11 
by  the  functional 

jp(u,v)  =*  |  vSp(On(u))  |vr[ds 

rc 


The  proof  of  this  equivalence  differs  in  some  aspects  from  the  one  given 
in  the  Appendix;  it  can  be  found  in  the  unpublished  report  by  ODEN  and 
PIRES  [1981]. 

Case  III  (p  =  e  =  0) .  When  both  p  =  0  and  e  =  0,  the  corresponding 
friction  law  is  given  by  the  conditions  5.2. 18.  The  system  5.3.3  -  5.3.6 
and  5.2.18  will  then  correspond  to  the  Signorini  problem  with  Coulomb 
friction.  DUVAUT  and  LIONS  [1976]  derived  a  variatoual  principle  charac¬ 
terizing  this  problem  which  is  given  by  5.3.11  if  we  replace  the  func¬ 
tional  j  (•**)  by  the  functional 
P»  t 

jCu.v)  -  J  v  l°n^l  l^x^ds 


Finally  we  mention  as  a  last  special  case,  an  auxiliary  problem 
that  proves  to  be  useful  in  the  next  Section  when  we  establish  the  con¬ 
ditions  for  the  existence  of  solutions  to  the  nonlocal-nonlinear  friction 
problem.  This  auxiliary  problem  involves  a  friction  law  for  which  the 
critical  or  limiting  value  of  the  tangential  stress  is  prescribed  rather 
than  being  determined  by  the  equilibrium  displacement  field  u  .  We  will 


then  have  in  the  nonlinear  case  (e  >  0)  ,  a  law  of  friction  of  type  5.2.12 


where  now  x  is  given  on  the  contact  surface  T  . 

v 


In  the  case  e  =  0 
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this  friction  law  assumes  the  form: 

|o,p(u)  [  <  T  implies  =  0 

lo^Cu)  |  =  T  implies  the  existence  of  X 

such  that  =  -XcJ. 

~T  ~T 

Thus,  depending  on  e  being  strictly  positive  or  zero,  we  have 

Case  IV  (t  fixed,  £  >_0)_.  The  corresponding  Signorini  problem  now 
consists  of  5. 3 . 3  -  5.3.6  together  with  the  friction  condition  5.2,12  (t 
given  on  r^,) .  The  equivalent  variational  formulation  of  this  problem  is  : 

Find  a  displacement  field  u  in  K  such  that 

+  j0  ^  ~  J0  ^  f (v-u)  (5.3.13) 

e  e 

for  all  v  in  K 

where 

jn  -  I  T  *  (|v  j)ds  (5.3.14) 

o£  Jr<_ 

Also,  it  is  not  difficult  to  show  that  in  this  case (5.3. 13)  is  equivalent 
to  the  constrained  minimization  problem  of  finding  u  in  K  such  that 

I  (u)  <  I  (v)  (5.3.15) 

£  —  C  ~ 

for  all  v  in  K  ,  where  the  energy  functional  I^C*)  is  defined  by 

1  <v)  -  ~  a(v,v)  -  f(v)  +  jQ  (v) 

e 


(5.  3.16) 
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Case  V  (t  fixed,  e  =  0).  The  law  of  friction  is  now  of  the  form  (5.3.12) 
which,  when  added  to  the  system  5.3.3  -5.3.6  ,  gives  the  corresponding 

Signorlni  problem.  We  emphasize  here  that  condition  5.3.12  represents 
a  law  of  friction  different  from  Coulomb's  law.  It  is  easy  to  establish 
the  equivalence  between  5.3.3  -  5.3.6  together  with5.3i2  and  5.3.13  if 


we  replace  in  5.3.13  jg  (*)  by  the  functional 

e 


i0<v>  3  |  Tl^T!ds 


(5.3.17) 


The  energy  functional  for  this  case  is  defined  by 


I(v)  =  -  a(v,v)  -  f(v)  +  Jq(v)  , 


v  in  K 


(5. .3. 18) 


Then  problem 


with  jg  (•)  replaced  by  3q(*)  is  equivalent  to  the  prob- 
e 


lera  of  seeking  minimizers  of  the  functional  in  5.3,18  which  satisfies  the 
unilateral  constraints. 


Remark  3.1:  It  is  interesting  to  note  that  when  e  =  0  the  functional 
jg(’)  defined  in  5.3.16  is  non-dif f erentiable  while  the  nonlinear  func¬ 
tional  jg  (•)  given  by  5.3.14  is  differentiable  on  all  of  V  .  Its 
e 

derivative  Djg  (•)  at  the  point  u  ,  in  the  direction  v  is  given  by 

e 


’’Jo/iP-  i  ■  |r  T  ds 


(5.3.19) 
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It  is  then  possible  to  show  that  the  variational  inequality  5.3.13  is  equi¬ 
valent  to  the  variational  inequality 

a(u,v-u)  +  Dj  (u)  •  (v-u)  >  f(v-u)  for  all  v  in  K  (5-3.20) 

-  oc  -  -  -  -  -  ~ 

The  differentiability  of  the  functional  j  (•)  has  important  impl icat ions 

E 

when  we  wish  to  consider  finite  element  approximations  of  the  friction 
problems  discussed  so  far.  In  fact,  the  direct  approximation  of  the  varia¬ 
tional  problem  considered  in  Case  V  (£  =  0)  by  finite  elements  through  the 
minimization  of  the  functional  defined  in  5.3.18  leads  to  a  discrete  sys¬ 
tem  for  which  the  most  popular  methods  for  solving  nonlinear  variational 
inequalities  do  not  apply,  owing  to  the  non-differentiability  of  Jq(*) 

5.3.4  Estimate  of  the  Difference  Between  the  Solutions  of  the  Friction 

Problems  with  e  >  0  (Nonlinear)  and  S  =  0.  We  wish  to  record  here  an 

estimate  of  the  difference  between  the  solutions  of  problem  5.3.11  and  of 

the  problem  considered  in  Case  II  (e  *  0)  of  the  previous  Section,  as  a 

function  of  E  .  The  first  result  is  concerned  with  the  approximation  of 

the  functional  j  (*,*)  by  the  functional  j  (•»•). 

P  P  *  E 

Proposition  3.2:  For  a  given  element  u  in  K  for  which  an(u) 
is  well  defined  there  exists  a  constant  c  >  0  independent  of  E  ,  such 
that 

|jp.e(u,v)  -  jpC«,y)  |  £c  e  (5.3.21) 

for  all  v  in  K  .  Thus,  _  approximates  j  arbitrarily 

closely  as  E  -*■  0  . 0 
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This  result  constitutes  the  basis  for  the  proof  of  the  following 
estimate: 

Proposition  3.3:  Let  u£  denote  a  solution  of  5.3.11  for  fixed 
C  >  0  and  let  u  be  a  solution  of  the  corresponding  variational  in¬ 
equality  obtained  by  setting  e  =  0  in  5.3,11.  Then,  for  a  sufficiently 
small  coefficient  of  friction  V,  there  exists  a  constant  k  >  0  ,  in¬ 
dependent  of  G,  such  that 

llu^  -  u  ||  _<  k  Jz  (5.3.22) 

where  ||  •  ||  is  the  norm  given  in  5. 3.8-0 

Results  similar  to  5.3.21  and  5.3.22  can  be  easily  derived  for  the 
cases  p  =  0  and  T  fixed  (given)  on  Tc  . 

5.4  Existence  and  Uniqueness  of  Solutions  to  the 
Nonlocal  Nonlinear  Frict ion  .Problem 

We  shall  now  establish  conditions  sufficient  to  guarantee  the  exis¬ 
tence  of  solutions  to  the  nonlocal  nonlinear  problem  5.3.11  as  well  as  ad¬ 
ditional  requirements  which  provide  for  uniqueness  of  solutions. 

We  begin  by  considering  an  important  preliminary  result  concerning  the 
auxiliary  problem  5. 3. 13  introduced  in  the  previous  Section, 

Proposition  4.1:  Given  T  _>  0  on  ,  T  smooth  enough, 

(i)  there  exists  a  unique  solution  u  to  5.3.13; 

(il)  the  correspondence  that  gives  for  each  T  the  corresponding 
solution  u  of  5. 3. 13  defines  a  continuous,  nonlinear  map  B,  B(t)  =*  u  ; 
(iii)  the  normal  contact  stress  produced  by  the  displacement  u  , 


o  (u)  =  E  ,  ,  u  .  n  n  is  well  defined,  and  is  continuous  as  a  function 
n  -  ij  k£  k,  £  i  j 

of  u  .  0 

We  next  state  our  second  major  result  which  concerns  the  existence  of 
solutions  to  the  nonlocal  nonlinear  friction  problem: 

Proposition  4.2:  For  the  smoothing  operator  defined  in  5.2.11 
which  transforms  the  normal  contact  stresses  0^  into  regularized  ones, 
there  exists  at  least  one  solution  u  in  K  of  the  nonlocal  nonlinear 
friction  problem  5.3.11  for  each  choice  of  smooth  enough  data  f  and  t  .0 

One  method  of  proof  of  this  proposition  was  suggested  by  DUVAUT  [1980] 
A  complete  proof  for  the  nonlocal  friction  problem  (c  =  0)  is  given  in 
ODEN  and  PIRES  [1981].  Only  the  general  structure  of  the  proof  is  of  in¬ 
terest  here  since  it  suggests  a  means  for  actually  calculating  solutions 
of  the  Signorini  problem.  The  key  steps  in  the  proof  are  outlined  as 
follows  :x 

1.  Pick  an  arbitrary  smooth  enough  T  >  0  on  . 

2.  By  Proposition  4.1,  for  each  such  T  ,  there  exists  a  unique 
solution  to 

a(u,v-u)  +  J  T  (*e(|vT|)  -  ^c(|uT|))ds 

(5.4.1) 

>_  f  (v-u) 

for  every  v  in  K  . 

3.  Let  u  =  B(t)  ,  where  B  is  defined  in  Proposition  4.1. 

4.  Compute  a^fu)  . 

5.  Calculate  VS  (o  (B(t)))  and  check  if  it  is  equal  to  T  .  If  so, 


we  can  write 
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a(u,v-u)  +  |  vSp(0n(u))  (4>£(|vt|)  -  \pc(  ju.fj))ds  >  f  (v-u) 

rc 

for  every  v  in  K  .  If  not,  we  go  back  to  step  2  and  replace 

in  5.4.1  t  by  the  new  value  VS„(o  (B(x)>).  An  iterative 

P  n 

scheme  is  therefore  obtained  by  repeating  in  this  way,  steps  2 
through  5. 

6.  Obviously,  step  5  describes  a  fixed  point  problem  for  the  operator 

T  =  VS_  o  a  o  B 
P  n 

We  must  therefore  show  that  there  exists  at  least  one  element 
X*  such  that  T(t*)  =  x*  .  Then  u*  =  B(x*)  will  be  a  solution 
of  the  contact  problem  with  nonlocal  nonlinear  friction. 

For  small  V  ,  the  composition  T  defined  above,  becomes  a  contrac¬ 
tion  and  the  fixed  point  is  unique.  Hence,  we  can  state  the  uniqueness 
result: 

Proposition  4.3:  If  the  coefficient  of  friction  is  sufficiently 
small,  the  nonlocal  nonlinear  friction  problem  5.3.11  possesses  a  unique 
solution.  □ 

Future  Work.  Numerical  solutions  of  the  nonlocal-nonlinear  friction 
problems  considered  in  this  paper  are  currently  under  investigation. 

These  include  studies  of  behavior  of  the  solution  for  various  values 
of  the  major  parameters;  the  coefficient  of  friction  V,  the  nonlocal 
contact  parameter  p,  and  the  tangential  stiffness  of  the  junctions  e  . 

This  work  is  to  be  the  subject  of  a  forthcoming  paper. 
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APPENDIX 


PROOF  OF  PROPOSITION  1 


Let  u  be  a  sufficiently  smooth  solution  of  the  Signorini  problem 
with  nonlocal  nonlinear  friction  (3.3)  -  (3.7)  .  Then,  the  following 
inequality  holds  for  every  admissible  displacement  v  : 

<?T(u)  *  (vT  -  uT)  +  vS  (°n(u))(^e(|vT|)  -  ¥e(|ut|)  >  0  on  Tc  (A.l) 


In  fact,  since  the  function  ){)(•)  is  convex  and  differentiable. 


*;duTl>  *  (^t  ■  ~t}  -  ~ 


-T  ~T  —  Te  ‘~T 1 


Ml«T|)  TTT  ’  (~T  "  -  ^e^~T^  ■  ^duTl> 

|  U  -  i  ... 


Hence,  since  u  is  such  that  (3.7)  holds, 


aT(u)  •  (vT  -  uT)  =  -vsp(an(u))  <f>£(|uT|)  *  (^T  "  «T> 


-  “  r  *edyTl)  -  ^c(|uT|)] 


oT(u)  •  (vT  -  uT)  +  (|»e(|vT|)  -  ^e(  l«T  |)  1  0 

Since  we  assume  u  to  be  sufficiently  regular  (e.g.  u  (  H  (ft)  0  K  ), 


i  -  2  2  N 

(ft)  H  K  =*  intersection  of  the  Sobolev  space  H  (ft)  =  (H  (ft))  = 

H2 (ft)  x  H2(ft) x  •  •  •  xH2(ft)  with  the  set  K.  For  the  definition  of  the  space 
H2(ft)  see  ADAMS,  R.  A.,  Sobolev  Spaces,  Academic  Press,  New  York,  1975. 


the  following  Green  s  formula  holds  for  every  v  in  K  : 

a(u,v-u)  =  o  (u)(v  -  u.  . )dx 

'  '  "  jq  xj  *•  *-*j 

•  'in°ij<-).j<vi  -  ui,dit 

+  o..(u)n.(v.  -  u  )ds 
J  r  1 J  ~  J  1  i 

But  o.j (u)  ^  =  -f^  in  fi  by  (3.3)  and 

lr°ijC~)nj(Vi  "  Ui)dS  =  i  “  Ui)ds  + 

rF 

+L  °ij  cs)nj  <vi  -  v is 

c 

for  any  v  in  K  since  u  =  0  on  T  by  (3.4)  and  0,.(u)n,  =  t. 

~-D  ij~ji 

T  by  (3.5).  Therefore 
r 

a(u,v-u)  -  f(v-u)  =  |r  0ij(u)nj(vi  -  u±)ds 

C 

=  (  [o  (u)  •  (v  -  u  )  +  0 (u)(v  -  u  )]ds 

Jr  ~  ~T  ~T  n  -  n  n 

C 

Here  vn  “  v  *  n  an<*  similarly  for  .  Hence  for  every  v  in  K 
a(u,v-u)  +  |  vSp(an(u))C^£(lyT|)  -  ^e(|uT|))ds  -  f (v-u)  = 

rc 

“  j  *  (~T  “  +  vSp<°n(u))('J'e(|vTl)  “ 

rc 

-^£(  |uT|)) +0  (u)  (v  -  u  )]ds 


Bi 


O  (u)(v  -  u  )  =  O  (u)[v  -  s  -  (u  -s)] 

n~n  n  n-n  n 

=  0  (u)  (v  -  s) 
n  ~  n 

>  0 

by  (3.6)  and  the  definition  of  the  set  K  .  Hence  applying  inequality 
(A.l)  we  finally  obtain 

a(u,v-u)  +  j  (u,v)  -  j  (u,u)  _>  f(v-u) 

K )  ^  ““  “**  ** 

for  every  v  in  K  ,  i.e.,  variational  inequality  (3.11)  and  the  first 
part  of  the  proposition  is  proved. 

Conversely,  let  u  in  K  be  a  solution  of  the  variational  inequality 
(3.11).  Following  the  proof  given  by  ODEN  and  PIRF.S  [1981,  p.  20-22]  for 
the  case  G  =  0  and  omitting  algebraic  manipulations,  we  are  led  to  con¬ 
ditions  (3.3)  -  (3.6)  if  we  interpret  them  in  the  sense  of  distributions. 
Then,  variational  inequality  (3.11)  will  be  reduced  to 

lT  lonM  *  (vT  -  uT)  +  vSp(an(u))(^(|yTj)  - 

-VlyTl»Jds  -  0 

for  every  v  in  K  .  If  we  let  vT  to  be  of  the  form 

Jr  '  !t  +  0(?t  -  Jr* 

where  0  belongs  to  the  open  interval  (0,1)  and  w^  is  the  tangential 
component  on  T  of  an  arbitrary  disp'  .cement  w  in  K  ,  we  obtain 

v  — 
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|  *  (~T  "  V  +  VSp(an(~))(','e(l”T  +  0(^t“~T)I) 

rc 

-^e( |uT|))]ds  >  0 


Dividing  through  by  8,  taking  the  limit  as  6  -*■  0  and  noticing  that 
is  differentiable  gives 


[°tM  *  (wT 


-  «T)  +  vSp(an(u))  <J>E(|uT|) 


(w?  -  «T)]ds  _>  0 


for  every  w  in  K  .  Finally  taking  first  w_  =  -u  and  then  w  =  2u 
produces 


°T(“)  *  ~T  +  VSp^°n^~^  1~T I?  I  =  0  on  rC 

which  we  have  seen  in  Section  3.1  to  be  equivalent  to  the  law  of  friction 
(3.7).  Q 


6.  SUMMARY  OF  STABILITY  RESULTS  FOR  REDUCED- INTEGRATION- 


PENALTY  METHODS 

A  discussion  of  the  °ccuracy  and  numerical  stability  of  several 
reduced-integration-penalty  methods  for  the  analysis  of  Stokesian 
flow  in  two  dimensions  is  presented.  A  summary  of  results  on  analytical 
studies  of  the  I.BB  condition  is  recorded.  Recommendations  on  which 
elements  provide  good  accuracy  .and  stability  for  use  in  computational 
fluid  dynamics  are  given. 


6.1  Introduction 

In  this  communication,  some  of  the  numerical  and  theoretical 
results  we  have  obtained  over  the  last  several  years  on  reduced- 
integration-penalty  (RIP)  methods  shall  be  summarized.  Complete  proofs 
and  more  detailed  discussions  can  be  found  in  ref erences  ODEN  et  al  [1980, 
1981,  1982]. 


The  basic  problem  to  be  considered  here  is  Stokes'  problem  for 
steady  confined  flow  of  a  viscous  incompressible  fluid,  which  can  be 
characterized  bv  the  iollowing  variational  boundary-value  problem: 

Find  a  velocity  field  u  t  V  and  a  hydrostatic  pressure 
field  p  (  Q  such  that 


a(u,v)  -  ( p , d iv  v )  =  f (v) 
(q,div  u)  =  0 


\  (6.1.1) 


for  all  admissible  velocities  v  in  V  and  all  admis- 
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Here, 


V  =  space  of  admissible  velocities 

=  (v  =  (vj»v2^  I  V£  £  hJ(£1)  ,  fi  a  regular  open  bounded  domain  in  K 
a(u,v)  =  the  virtual  work,  a  continuous  bilinear  form  on  V 


=  v 


Q 


grad  u  :  grad  v  dx 


3Ui  9v^ 

ja  ^  dx;  dx  =  dxidx2 

ij=i  j  j 


-vfo  S 


V  being  the  viscosity  (V  >  0) 


(p,q)  = 


n 


pqdx  =  L  (Q) -inner  product 


Q  =  space  of  hydrostatic  pressures  =  L  (£1) 

=  space  of  Lagrange  multipliers  corresponding  with  the  constraint 
"div  u  =  0  in  fi  " 
f (v)  =  virtual  work  of  body  forces  f  = 


Q 


f  •  v  dx  ,  v  an  arbitrary  "virtual"  velocity  in  V 


We  endow  V  with  the  energy  norm. 


II  Y  Hv 


^—7  9v^  9v^ 

n  8x  3x,  dx 

i,j=l  J  J 


and  Q  with  the  usual  L  -norm. 


When  the  bilinear  form  a(*,*)  satisfies  the  conditions 


a(u,v)  <  M  ||u  |!v  |( v  |iv;  a(v,v)  >  n  ||v  jj“  (fi.1.2) 

for  arbitrary  velocities  u,v  in  V  (with  M  and  a  positive  constants) 
and  if  there  exists  a  constant  S  >  0  such  that 


6  111  llQ/E  i.  S“P 


|  (q,div  v)  [ 


(6.1.3) 


for  all  q  in  Q  ,  then  we  are  guaranteed  the  existence  of  a  unique 

2 

solution  u  in  V  and  p  in  L  (Q)/]R  to  problem  (6.1.1)  [here 
2 

L  (Q)/!R  =  Q/1R  is  the  quotient  space  of  classes  of  square-integrable 
functions  "modulo  constants";  i.e.  the  hydrostatic  pressure  can  be 
determined  only  to  within  an  arbitrary  real  constant]. 

Under  suitable  conditions  on  the  data,  solutions  of  (6.1.1)  are  also 
solutions  of  the  classical  Stokes'  problem. 


-Au  +  Ap  =  f 
u  =  0 
div  u  =  0 


in 

in 

in 


(6.1.4) 


It  has  been  observed  by  many  authors  (see  in  particular  REDDY  [1978,1979], 
HUGHES  [1977],  MAI.KUS  [1977],  and  the  works  of  ODEN  and  associates  [1979-1980]) 
that  considerable  computational  effort  might  be  saved  by  employing  an  ex¬ 
terior  penalty  approximation  of  the  constraint 
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div  u  =  0 

In  fact,  many  ideas  of  convex  optimization  surface  when  we  consider  the 
fact  that  problem  (4)  is  equivalent  to  the  problem  of  minimizing  the 
functional 

J  :  V  -*•  R  ;  J(v)  =  I  a(y»v)  ~  £(v) 
on  the  linear  subspace,. 

K  =  (v  in  V  |  div  v  =  0  (in  Q)} 

One  can  relax  such  constraints  by  appending  to  J  a  convex,  differentiable 
penalty  term  P(v)  =  (2e)  ^(div  v,  div  v)  so  as  to  produce  the  penalized 
functional,  for  e  >  0  ,  given  by 

Je  :  V  -  R  ;  Je(y)  =  J (v)  +  ~  |! div  y  \\*  (6.1.5) 

Minimizers  u  of  J  are  characterized  by, 

~e  e 

u£  £  V  :  a(ue>v)  +  £  ^(div  u£,  div  v)  =  f(v) 

for  all  v  in  V  (6.1.6) 

It  is  informative  to  note  that  an  alternative  way  at  arriving  at 
the  same  formulation  (6.1.6)  is  to  use  the  so-called  perturbed  Lagrangian 
method.  The  Lagrangian  associated  with  J  and  the  incompressibility 


constraint  is 
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L  :  V  x  Q  -*■  1R  ;  L(v,q)  =  J(v)  -  (q,  div  v) 
The  perturbed  Lagrangian  is  defined  by 


LG(v,q)  =  L(v,q)  -  (q,q) 


and  its  saddle  points  (u^.p^)  are  characterized  by  the  system 


a(ue»v)  -  (p  div  v)  =  f(v)  for  all  v  in  V 


(ep£  +  div  u^,q)  =  0  for  all  c  in  Q 


Then  one  immediately  has 


(6.1.7) 


p  =  -£  ^div  u  in  Q 

e  ~e 


(6.1.8) 


and,  hence,  the  first  member  of  (6.1.7)  reduces  to 


a(u£,v)  +  e  ^(div  u£,  div  v)  =  f (v)  for  all  v  in  V 


which  is  precisely  (6.1.6).  Hence  (6.1.7)  and  (6.1.6)  are  equivalent 
formations,  but  (6.1.8)  provides  a  method  for  also  calculating  approxi¬ 
mations  of  the  hydrostatic  pressure  from  the  penalty  approximation  of 
the  velocity  field. 

6.2  Full-Integration  Penalty  Methods 
To  construct  a  finite  element  approximation  of  problem  (1)  [or, 
equivalently,  of  (4)],  we  proceed  in  the  usual  fashion  by  replacing  Q 
by  a  mesh  consisting  of  a  collection  of  E  finite  elements. 
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Continuous  piecewise  polynomial  approximations  of  the  velocities  over 
are  denoted  v^  .  If  is  polygonal,  we  can  usually  construct 
so  that  the  approximate  velocities  be  in  a  finite-dimensional  subspace 

V*1  of  the  space  of  admissible  velocities  V  .  We  generally  denote  by  h 
the  mesh  parameter 


h  =  max  h  ;  h  =  dia(£2  ) 
l<e<E  e  e  e 


where  fi  is  a  finite  element  in  ,  and  by  regular  refinements  of  the 
e  h 

mesh  we  generate  a  family  {V*1}  of  subspaces,  the  union  of  which  is 

h>U 

everywhere  dense  in  V  . 

It  is,  of  course,  also  possible  to  introduce  a  space  QhO  of 
approximate  hydrostatic  pressures  defined  over  the  mesh  £1^  *  these 

being  piecewise  polynomials  not  necessarily  continuous  across  interelement 
boundaries.  But  the  spaces  Q*1  are  not  generally  explicit  in  a  penalty 
approximation  of  (6.1.1).  We  shall  show  below  that,  in  fact,  these  spaces 
are  always  inherently  defined  by  the  manner  in  which  one  approximates 
the  penalty  functional. 

The  first  method  that  one  uses  to  approximate  (6.1.1)  is  as  follows: 


e  /»  h 

For  given  e  >  0  ,  find  u  V  such  that 


rh 

a^Uh’Vh^  +  e_1(div  uh’  div  =  f^~h 


: 


(6.2.1) 


for  all  v.  in  v 

~h  S 


Most  engineers  seem  to  think  method  (6.2.1)  "will  not  work"  and  that 


it  leads  to  a  "locked  solution"  (u^  *  0  as  c  +  0)  ,  but  this  is  not 
the  case.  The  problem  is  that  if  the  terms  in  (6.2.1)  are  integrated  exactly, 
then  the  stability  of  the  method  is  conditional,  depending  on  the  re¬ 
lationship  between  £  and  h  .  For  instance,  Falk  [1975]  and  Falk  and 
King  [1975]  have  considered  finite  element  approximations  of  the  type 
(5.2.1)  with  a  penalty  parameter  of  the  form 


where  y^  and  o  are  positive  constants,  with  a  to  be  chosen  so  that 

an  optimal  rate  of  convergence  is  obtained.  It  is  interesting  to  note 

that  no  such  0  exists  and  that  post  processing  using  an  extrapolation 

technique  is  needed  to  achieve  the  optimal  rate.  However,  without 

2 

extrapolation,  the  best  choice  of  0  for  their  method  is  a  =  -(s-1)  , 

£_2  2 

where  s  =  min(B,k)  and  the  body  force  data  f  is  in  (H  (ft))  , 

B  >  2  ,  and  k  is  the  degree  of  the  complete  polynomial  approximation 
of  the  velocities.  The  point  is  that  these  methods  can  converge,  but 
only  if  £  is  taken  as  a  special  function  of  h  and  even  then  at  a 
suboptimal  rate. 

The  problem  with  method  (6.2.1)  is  that,  from  a  practical  point  of  view, 
it  is  not  satisfactory  because,  for  a  reasonably  fine  mesh,  £  must  be 
taken  so  large  to  produce  a  non-degenerate  solution  that  the  incompress¬ 
ibility  constraint  is  not  adequately  satisfied.  The  mesh  sizes  needed 
to  make  this  method  attractive  are  so  small  that  the  computational  effort 
needed  to  extract  a  solution  is  prohibitive. 
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6-3  Reduced  Integration 


To  overcome  the  difficulties  mentioned  above,  it  has  become  common 
practice  to  use  a  numerical  quadrature  scheme  for  evaluating  the  penalty 
integral  [  div  u  div  v  dx  which  is  of  lower  order  than  that 
required  to  integrate  this  term  exactly.  Let 


I(fg)  = 


I  Ie(fg) 
e=l 


(fg  £  C°  (Q)) 


(6.3.1) 


I  (fg)  =  £  W®  f(£*)  gtf*) 

j=l  J  "J 


denote  a  numerical  quadrature  rule  for  integrating  the  product  of 

0 

piecewise  continuous  functions  fg  over  the  mesh,  where  >  0  are 
the  quadrature  weights  for  element  e  and  are  the  quadrature 

A 

points  in  element  e  .  Suppose  that  a  number  G  of  such  points  must 

be  used  to  integrate  the  functions  div  u,  div  v,  (u  ,v.  £  V*1) 

-n  - n  -n 

A 

exactly  and  that  G  <  G  .  Then  a  r educed- integration-penalty  ap¬ 
proximation  of  problem  (6.1.1)  consists  of  solving,  instead  of  (6.2.1), 
the  following  discrete  problem: 


£ 

For  e  >  0  ,  find  u,  in  V  such  that 

~  n 


a(u^,vh)  +  e  1I(div  u^  div  v^)  =  f (v^) 


(6.3.2) 


for  all  v  in  V 
~h 


Corresponding  to  each  choice  of  V  and  !(*)  there  is  uniquely 
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defined  by  (6,3.2)  a  finite-dimensional  space  0  Q  of  hydrostatic 

pressures  q  .  The  unique  element  p,  in  Q*1  defined  by 
n  n 


£  £ 

I((£p  +  div  u,  ) q  )  =  0  for  all  q,  in  Q, 

n  ~h  n  n  h 


(6.3.3) 


or,  equivalently. 


£ /^e\  -1 ,  .  £  ,,-e, 

ph(5j)  '  ~c  dlv“h<5J) 


(6.3.4) 


is  the  corresponding  approximation  of  the  pressure  p  .  In  fact,  (6.3.2) 

and  (6.3.3)  correspond  to  (6.1.7)  with  (’,*)  replaced  by  I(‘).  For  instance, 

if  V*1  consists  of  piecewise  polynomials  on  rectangular  elements  and 

!(•)  is  2  x  2  Gaussian  quadrature,  Q*1  then  is  spanned  by  piecewise 

bilinear  functions,  discontinuous  across  interelement  boundaries,  with 

nodes  at  the  four  Gaussian  quadrature  points.  If  I(*)  is  one  point 

integration,  Q  is  then  spanned  by  piecewise  constraints,  etc. 

The  key  issue  is  whether  or  not  solutions  of  (6.3.2)  exist  and,  if  so, 

how  they  behave  as  h  tends  to  zero.  Let  Div,  and  V  be  discrete 

n  h 

divergence  and  gradient  operators  defined  by 


Div,  :  Vh  ->  Qh  ;  V,  :  Qh’  -*■  Vh’ 
n  n 


<<!h  •  Dlvh  v =  I(,,h dlv  V  *  -[Vh  >  V 


(6.3.5) 


ti  h 

for  all  in  V  and  q,  in  Q  Then,  in  view  of  the  results  listed 

~h  ^h 

in  Section  1,  a  unique  solution  to  (6.3.2),  (6.3.3)  exists  whenever  (6.1.2) 

holds  for  all  u,  ,v,  in  V*1  (which  is  always  true  if  (6.1.2)  holds  and 
~n  -n 
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and  VhC  V)  and  if,  in  analogy  with  (6.1.3),  the  following  discrete 
Babuska-Brezzi  condition  holds: 


There  exists  a  constant  >  0  such  that  for 
all  hydrostatic  pressures  q^  in  , 


Ph  ^h  Q/ker  7  ±  SUp 
n 


I(qhdiv  vh) 


lh  11 1 


(6.3.6) 


Here 


ker  Vh  =  {q^  in  Q  |  l(q^div  v^)  =  0 


for  all  v  in  V  } 
~n 


(6.3.7) 


Error  estimates  can  be  obtained  in  certain  caser .  Suppose  that, 
as  £  0  ,  (6.3.2)  and  (6.3.3)  lead  to  the  mixed  finite  element  problem 


-  1<phdlv  V  * f  'V 

I (q^div  u^)  *  0 


(6.3.8) 


for  arbitrary  vh  and  qh  .  Let  EI  denote  a  generic  integration 
error  defined  by 


Setting  v 


Ei  ■  EiViv  V  ■  V  dlv  V  -  I(t,hdiv  V 

v,  in  (6.1.1)  and  substracting  (6.3.7)  gives 
~h 


(6.3.9) 


0 


(6.3.10) 


a(u  "  V  ~h)  "  (p  "  V  diV  ~h  "  EI  = 


Using  (6.1.1)  and  (6.3.10),  we  have 


2 

a  11“  '  “h  11  v  -  a(~  “  V  ~  “  ~h) 


=  -  v  "  "  V  +  a(“  ■  V  Yh  '  “h} 


<  M  ||  u  -  uh  |iv  ||u  -  vh  Hv  +  (p  -  ph,div(vh  - 


<  fllu 


“h  11 V  +  2a  II  “h  ~  ~h  II V 


V  +  EI 


+  lip  -  Ph  llqt  11“  -  “h  1IV  +  ’  Yh  II  v1  +  Ei 


Thus, 


U  -  uh  II  <  c(  inf  ||u  -  vh  |jv  +  Up  -  Ph  II 0)  +  Ex 

v,  in  V 
~n 


(6.3.11) 


P  -  Philo  i  ‘"£  h  f  Up  -  IjJq  +ilph  -  II 


,h  in  Q 


(6.3.12) 


i  inf  h  IIP  -  PhliQ  +HPh  *  %  Ho/ker  A 

q,  m  Q  x  1 


i>h- vwvh-ir:,E 


I  a (u  -  uh»yh>  +  Ejt 


u  -  u. 


Ei  (ph-dlv;h) 

n-hHv 


etc.  Again  we  note  that  the  convergence  and  stability  of  the  method 

are  strongly  tied  to  the  constant  8^  ,  the  characterization  of 

ker  V  ,  and  the  interpolation  properties  of  Q*1  and  V 
h 


6.4  Summary  of  Some  Stability  Results 


A  mathematical  analysis  of  the  discrete  Babuska-Brezzi  condition 

(6.3.6)  has  been  made  by  Oden  and  Kikuchi  [1982],  Oden,  Kikuchi  and 

Song  [1981]  and  Oden,  Jacquotte  [1982]  for  several  finite  elements  for  a 

two-dimensional  Stokes'  problem  on  a  uniform  mesh.  We  shall  summarize 

chese  results  here  which  pertain  to  the  behavior  of  the  "LBB-constant 

8  and  the  stability  of  the  pressure  calculations.  We  use  the 
h 

notations 

=  space  of  complete  piecewise  polynomials  of  degree  k 
over  an  element 


=  space  of  tensor  products  of  complete  polynomials  of 
degree  k 


18  =  the  eight-node  isoparametric  element 

Results  are  summarized  in  Table  1.  In  this  table,  figures  1,  2, 

and  7  "lock”  at  small  values  of  the  penalty  parameter  c  .  This  means 

that  for  a  given  mesh  size  h  ,  e  cannot  be  taken  arbitrarily  small, 

as  noted  earlier.  Of  course,  for  an  acceptable  c  for  reasonable 

mesh  sizes,  e  is  so  large  that  the  constraint  of  incompressibility 

is  not  adequately  satisfied.  Hence  these  elements  should  generally  be 

avoided.  Elements  2,  4,  5,  8,  11,  and  14  are  unstable  since  B,  =  0(h)  . 

h 

Remarkably,  these  instabilities  frequently  are  not  observed  on  uniform 
meshes  when  the  solution  is  very  smooth.  Mild  irregularities  in  the 
solution  or  small  perturbations  in  the  mesh  may,  however,  produce  violent 
oscillations  in  computed  pressures  the  magnitudes  of  which  increase 
without  bound  as  h  tends  to  zero.  In  many  cases,  however,  these 
oscillations  disappear  upon  "filtering"  the  pressure  solutions  (i.e. 
upon  averaging  the  pressures  over  one  or  more  elements) .  In  the  case  of 
elements  2  and  14  it  has  been  proved  mathematically  (by  N.  Kikuchi  and 
the  author)  that  certain  filtering  schemes  will  produce  a  stable  and 
convergent  method.  However,  it  is  not  known  if  filtering  can  be  used 
to  stabilize  and  salvage  the  remaining  unstable  elements. 

Elements  6  and  10  lead  to  stable  and  convergent  schemes  and  are 
quite  robust  in  the  sense  that  they  are  insensitive  to  singularities  in 
the  solution.  However,  they  are  not  too  accurate  and  converge  at  a 
suboptimal  rate. 

Elements  5  and  9  are  calculated  using  the  perturbed  Lagrangian 
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Ideas  discussed  in  Section  1:  a  piecewise  linear  approximation  of  the 
regularized  pressure  p£  is  computed  over  each  element.  Then  (6.1.7) 
leads  to  a  discrete  system  for  each  element  of  the  form 


£?e 


f  +  a 


eMp  +  Bu£  =  0 


(6.4.1) 


where  K  is  the  element  stiffness  matrix,  B  the  non-rectangular 
constraint  matrix,  f  the  load  vector,  a  the  "connecting"  vector 
(which  sums  to  zero  upon  connecting  elements  together  to  form  the  mesh), 
and  M  is  the  Gram  matrix  corresponding  to  the  linear  shape  functions 
for  p_  .  Thus, 

p£  =  -e“1M“1BTuE  (6.4.2) 

£ 

i.e.,  the  fact  that  p,  is  discontinuous  across  interelement  boundaries 

h 

makes  it  possible  to  eliminate  the  pressure  at^  the  element  level  by 
(6.4.2).  Then  th  penalty  approximation  over  an  element  is  characterized 
by 

(K  +  £_1  BTM  B)u£  *  f  +  0  (6.4.3) 


Thus,  no  reduced  integration  rule  is  actually  used  in  constructing 
elements  5  and  9. 

Element  9  is  clearly  the  superior  of  any  listed:  it  is  unconditionally 
stable,  it  provides  both  velocity  and  pressure  approximations  which 
converge  at  the  optimal  rate,  and 
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ker  V,  =  ker  V 
h 

Element  13  is  somewhat  of  a  novelty.  While  element  5  yields 
unstable  pressure  approximations,  ODEN  and  JACQUOTTE  [4]  have  shown 
that  a  composite  of  three  Q?/P^  elements  (no.  9)  and  one  I8/P^ 
element  (no.  5)  is  stable. 

The  behavior  of  elements  11  and  12,  marked  with  an  asterik,  is 

only  conjectured  here  and  has  not  been  rigorously  proven. 

Extensions  of  these  results  to  three-dimensional  elements  are 
straightforward . 


7.  SUGGESTED  AREAS  FOR  FUTURE.  RESEARCH 


During  the  course  of  this  project,  two  areas  have  emerged  as  representing 
important  research  subjects  needing  significant  amount  of  additional  study. 
These  are:  1)  the  study  of  nonlinear  and  non-classical  friction  laws  in 
contact  phenomena  in  solids  and  structures  and  2)  a  continued  study  of  models 
and  numerical  techniques  for  the  analysis  of  finite  elastoplastic  deformations 
of  the  type  characteristic  of  metal  forming  processes. 

The  idea  that  new  descriptions  of  friction  are  necessary  to  adequately 
describe  the  phenomena  of  contact,  impact,  and  wear  of  deformable  bodies 
will  have  a  significant  impact  on  broad  areas  of  applied  mechanics  in  engin¬ 
eering.  It  will  mean  that  a  variety  of  new  models  and  results  will  need  to 
be  developed  to  adequately  describe  such  phenomena  as  load  reversal  on  contact 
surfaces,  heat  generation,  abrasion  and  wear,  impact,  dynamical  friction  ef¬ 
fects,  even  fracture  initiation  and  growth.  The  multitude  of  this  phenomena 
in  which  ver;  crude  friction  models  have  been  used  in  the  past,  must  ulti¬ 
mately  be  re-examined  in  some  detail.  This  will  represent  a  research  effort 
of  very  large  proportion,  but  should  ultimately  have  a  significant  pay-off 
in  terms  of  the  liability  of  mathematical  models  and  numerical  techniques  for 
simulating  the  nonlinear  behavior  of  complex  structures. 

The  reliability  of  most  of  the  numerical  simulators  of  metal  forming 
processes  is  very  much  in  doubt.  Recent  results  seem  to  indicate  that  most 
of  the  popular  methods  may  be  marginally  stable  and  produce  stress  approxi¬ 
mations  which  are  very  sensitive  to  perturbations  in  the  mesh  material  pro¬ 
perties.  However,  if  true,  this  would  be  a  very  undesirable  situation,  since 
these  factors  play  a  fundamental  role  in  the  prediction  of  residual  stresses 
in  machine  parts  and  structures.  Therefore,  a  careful  mathematical  and  nu¬ 
merical  study  of  mathematical  models  and  numerical  methods  for  handling  these 
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classes  of  problems  seems  to  be  in  order  and  to  represent  an  important 
are  for  future  research. 


REFERENCES 


ARCHARD,  .1 . F .  [1957],  "Elastic  Deformation  and  the  Laws  of  Friction,"  ?ioc . 

Roy.  Soc.  A243,  pp.  190-205. 

AUBIN ,  J.P.  [1979],  Applied  Functional  Analysis,  John  Wiley  &  Sons,  New  York. 

BABUSKA,  I.  and  AZIZ,  A.K.  [1972],  "Survey  Lectures  on  the  Mathematical 

Foundations  of  the  Finite  Element  Methods,"  The  Mathematical  Foundations 
of  the  Finite  Element  Method  with  Applications  to  Partial  Differential 
Equations,  Ed.  by  A.K.  Aziz,  Academic  Press,  New  York,  pp.  1-359. 

BOWDEN,  F.P.  and  TABOR,  D.  [1950],  The  Friction  and  Lubrication  of  Solids, 
Clarendon  Press,  Oxford. 

BOWDEN,  F.P.  and  TABOR,  D.  [1964],  The  Friction  and  Lubrication  of  Solids, 

Part  II,  Clarendon  Press,  Oxford. 

BREZIS,  H.  [1968],  "Equations  et  Inequations  Non  Lineaires  dans  les  Espaces 

Vectoriels  en  Dualite,"  Ann.  Inst.  Fourier,  Grenoble,  Vol.  18,  pp.  115-175 

DUVAUT,  G.  [1980],  "Problemes  Mathematiques  de  la  Mecanique  -  Equilibre  d'un 
solide  elastique  avec  contact  unilateral  et  frottement  de  Coulomb," 

C.R.  Acad.  Sc.,  Paris,  t.  290,  Serie  A,  pp.  263-265. 

DUVAUT,  G.  and  LIONS,  J.L.  [1976],  Inequalities  in  Mechanics  and  Physics, 
Springer-Verlag,  Berlin,  Heidelberg,  New  York. 

ERINGEN,  A.C.  (ed.)  [1976],  Continuum  Physics,  Vol.  IV  -  Polar  and  Nonlocal 
Field  Theories,  Part  II  by  D.G.B.  Edelen,  "Nonlocal  Field  Theories," 

Acad.  Press,  New  York. 

ERINGEN,  A.C.  and  EDELEN,  D.G.B.  [1972],  "On  Nonlocal  Elasticity,"  Int.  Journal 
of  Engineering  Science,  Vol.  10,  p.  233. 

FALK,  S.R.  [1974],  "Error  Estimates  for  the  Approximation  of  a  Class  of 

Variational  Inequalities,"  Mathematics  Computer,  Vol.  28,  pp.  963-971. 

FALK,  R.S.  [1975],  "An  Analysis  of  the  Penalty  Method  and  Extrapolation  for 
the  Stationary  Stokes  Problem,"  Advances  in  Computer  Methods  for 
Partial  Differential  Equations,  Edited  by  R.  Vichnevetsky ,  AICA 
Publication,  pp.  66-69. 

FALK,  R.S.  and  KING,  J.T.  [1975],  "A  Penalty  and  Extrapolation  Method  for 
the  Stationary  Stokes  Equation,"  SIAM  Journal  of  Numerical  Analysis, 

Vol.  13,  No.  5,  pp.  814-829.  '  . 

HLAVACEK,  I.  and  NECAS,  J.  [1970],  "On  Inequalities  of  Konn’s  Type," 

Arch.  Rational  Mech,  Anal.,  Vol.  36,  pp.  304-334. 

HUGHES,  T.J.R.  [1977],  "Equivalence  of  Finite  Elements  for  Nearly  Incompressible 
Elasticity,"  Journal  of  Applied  Mechanics. 

JOHNSON,  K.L.  [1955],  "Surface  Interaction  Between  Elastically  Loaded  Bodies 
Under  Tangential  Forces,"  Proc.  Roy.  Soc.  A,  Vol.  230,  pp.  531-548. 


166 


KIKUCH1 ,  N.  [  1932  ],  On  Regular  izat  ion  of  the  Coulomb  Friction  haw,  Report 
No.  3,  Department  of  Mechanical  Engineering  and  Applied  Mechanics, 
University  of  Michigan,  Ann  Arbor. 

K1KUCHI ,  N.  and  ODEN,  J.T.  [1983],  Contact  Problems  in  Elasticity,  SIAM, 
Philadelphia . 

LIONS,  J.L.  and  STAMPACHIA,  G.  [1967],  "Variational  Inequalities,"  Comm. 

Pure  Appl.  Math.,  Vol.  20,  pp.  129-148. 

MALKUS ,  D.S.  [1977],  "Finite  Element  Analysis  of  Incompri ssi ble  Solids," 
Dissertat ion ,  Department  of  Mathematics,  Boston  University. 

MINDLIN,  R.D.  [1949],  "Compliance  of  Elastic  Bodies  in  Contact,"  J_.  Appl . 

Mech. ,  Vol.  16,  pp .  259-268. 

NECAS,  J.,  JARU!>EK,  J.  and  HASLTNCER,  J.  [1980],  "On  the  Solution  of  the 

Variational  Inequality  to  the  Signorini  Problem  with  Small  Friction," 
Boll.  U.  M.  I.  (5),  17-B ,  pp.  796-811. 

ODEN,  J.T.  [1982],  "Finite  Element  Methods  for  Constrained  Problems  in 

Elasticity  and  Fluid  Mechanics,"  Proceedings  of  Symposium  on  Finite 
Element  Methods,  Hefei,  China,  Science  Press,  People's  Republic  of 
China,  (to  appear). 

ODEN,  J.T.  [1980]  "RIP  Methods  for  Stokesian  Flows,"  TICOM  Report  80-11, 

Austin,  [an  extended  version  is  to  appear  in  Finite  Elements  in  Fluids, 
Vol.  IV,  John  Wiley  &  Sons,  Ltd.,  1982]. 

ODEN,  J.T,  and  DEMKOWICZ,  L.  [1981],  On  Some  Existence  and  Uniqueness  Results 
in  Contact  Problems  with  Nonlocal  Friction,  TICOM  Report,  81-10,  The 
University  of  Texas  at  Austin. 

ODEN,  J.T.  and  JACQU0TTE,  0.  [1982]  ,  "A  Stable  Second-Order  Accurate,  Finite 
Element  Scheme  for  the  Analysis  of  Two-Dimensional  Incompressible 
Viscous  Flows,"  Proceedings  of  the  Fourth  International  Symposium  on 
Finite  Elements  in  Flow  Problems,  Tokyo ,  Japan. 

ODEN,  J.T.  and  KIKUCHI,  N.  [1980,1982],  "Finite  Element  Methods  for  Constrained 
Problems  in  Elasticity,"  TICOM  Report  81-10,  Austin,  (also  appeared  in 
International  Journal  of  Numerical  Methods  in  Engineering,  1982,  Vol.  18, 
pp.  701-725.  '  . 

ODEN,  J.T.  and  PIRES ,  E.B .  [1981].  Contact  Problems  in  Elastostatics  with 
Nonlocal  Friction  Laws,  TICOM  Report  81-12,  The  University  of  Texas 
at  Austin. 


ODEN,  J.T. ,  KIKUCHI,  N.  and  SONG,  Y.J.  [1981],  "Penalty-Finite  Element 
Methods  for  the  Analysis  of  Stokesian  Flows,"  TICOM  Report  81-11, 

Austin,  (to  appear  in  Computer  Methods  in  Applied  Mechanics  and  Engineering). 

RABINOWITZ ,  E.  [1965],  Friction  and  Wear  of  Materials,  John  Wiley  &  Sons,  Ltd., 

New  York. 


167 


REDDY,  J.N.  [1978],  "On  the  Accuracy  and  Existence  of  Solutions  to  Primitive 
Variable  Models  of  Viscous  Incompressible  Fluids,"  International 
Journal  of  Engineering  Science,  Vol.  16,  p.  921. 

REDDY,  J.N.  [1979],  "On  the  Finite  Element  Method  with  Penalty  for  Incompressible 
Fluid  Flow  Problems,"  The  Mathematics  of  Finite  Elements  and  Applications 
I IT ,  J.R.  Whiteman  (ed.),  Academic  Press,  London,  pp.  227-235. 

SHIBUYA,  T.,  KOIZUMI,  T.  and  NAKAHARA,  I.  [1980],  "An  Elastic  Contact  Problem 
for  a  Half-Space  Indented  by  a  Flat  Annular  Rigid  Stamp  in  the  Presence 
of  Adhesion,"  Proceedings  of  the  25th  Japan  National  Congress  for  Applied 
Meehan ics ,  Vol.  25,  University  of  Tokyo  Press,  Tokyo,  pp.  497-504. 


